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The System Lime-Alumina-Water at 1°C 
Elmer T. Carlson 


A study was made of phase equilibria in the system CaQO-Al,0;-H,0 at 1° C. 


The 


stable phases are believed to be gibbsite, tetracalcium aluminate hydrate, and Ca(OH)», 
but no solubility curve was established for gibbsite because of extreme slowness of reaction. 
he compound 3CaO-Al,O;-6H,0 appears to be metastable, being very slowly transformed 


to one or the other of the hexagonal hydrates. 
lished for CaQ-Al,O;-10H,O and 2CaQO.-A],O3-8H,0. 
(a metastable invariant point) is about 0.10 g of Al,O; and 0.38 g of CaO per liter. 


Metastable equilibrium curves were estab- 
he concentration at their intersection 
The 


dicalecium compound can readily be dried to a 6H,O stage, giving a characteristic X-ray 


pattern, 


During the long storage period required to approach equilibrium, most of the 


reaction mixtures picked up a little CQ), apparently by absorption through the polyethylene 
containers, with formation of a carboaluminate, probably 3CaQ-Al,O3-CaCO;-11H,0. 


1. Introduction 


The system lime-alumina-water, largely because 
of its relation to the setting and hardening of hy- 
draulic cements, is a field in which there has been a 
great deal of research over a period of many years. 
The literature on the subject has been thoroughly 
reviewed by Steinour (1951) [1],' by Bogue (1955) 
2], and by Lea (1956) [8]. Extensive studies of 
phase equilibria in the system at room temperature 
have been reported by Bessey (25° C) [4], by Wells, 
Clarke, and MeMurdie (21° C) [5], and by D’Ans 
and Eick (20° C) [6]. Equilibrium studies at tem- 
peratures from 50° to 250° C have been reported by 
Peppler and Wells [7] and at 100° to 1,000° C by 
Majumdar and Roy [8]. 

The phase diagrams given In the reports of the 
three investigations at room temperature are in fair 
agreement, in general; but all three are admittedly 
incomplete or uncertain in some respects, and on a 
few points there is direct conflict of opinion. Among 
the questions left unsetiled are the following: (1) 
The existence of a hexagonal tricalcium aluminate 
hydrate, usually formulated 3CaO-Al,O,-12H,O. By 
many authorities such a compound is considered 
“well-established”, and D’Ans and Eick [6] agreed 
with this view, though they did not attempt to show 
an equilibrium curve for it in their phase diagram. 
On the other hand, Wells, Clarke, and MeMurdie [5] 
concluded that the supposed compound was actually 
amixture of the dicaleium and tetracalcium alumi- 
nate hydrates. (2) The existence of more than one 
form of teiracalecium aluminate hydrate. This ides 
was firsi pui forth by Assarsson [9], who concluded, 
on the basis of slight differences in index of refraction 
and X-ray pattern, that two very similar forms, 
which he termed @ and 8, existed at room tempera- 
ture. D’Ans and Eick show widely different meta- 
stability areas for the two forms in their phase dia- 
gram, whereas Wells, Clarke, and MeMurdie report 
only one form. (3) The area of stability, or meta- 
stability, of monocalcium aluminate hydrate. Al- 
though Assarsson [9] reported the formation of such 
4compound many years ago, its existence has been 
tonfirmed only recently, and none of the published 
phase diagrams shows any stability area for it. 
tenes 
' Figures in brackets indicate the literature references at the end of this paper. 
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The present study was undertaken in the hope of 
shedding some further light on the above questions, 
among others. The temperature, 1° C, was selected 
because Assarsson [9] had indicated that the forma- 
tion of the monocalcium aluminate hydrate was 
favored by low temperature. Insofar as the infor- 
mation gained may be related to the practical prob- 
lems of hydraulic cements, there is the further con- 
sideration that temperatures near the freezing point 
occasionally are encountered during the placing and 
hardening of such materials. 


2. Experimental Work 
2.1. Procedures and Materials 


An ordinary household refrigerator, equipped with 
a special control device and a blower, served as a cold 
cabinet for storage of the reaction mixtures. The 
temperature was kept at 1° +1° C, the relatively 
wide range being chiefly due to variations from point 
to point within the chamber rather than to the cycling 
interval. This range of temperature was acceptable 
because the solubilities change very little with 
temperature. 

The primary starting materials were reagent-grade 
calcium carbonate, crystalline hydrated alumina 
(gibbsite), and distilled water. For most of the ex- 
periments the water was freshly boiled to expel car- 
bon dioxide. Secondary starting materials, namely, 
calcium aluminates, were prepared by heating 
mixtures of calctum carbonate and gibbsite in the 
desired proportions. For the greater part of the 
subsequent work, a monocalcium aluminate prepara- 
tion was desired that would yield a calcium aluminate 
solution of high concentration. This was obtained 
by heating the mix at 1,300° C. since preliminary 
experiments showed that the use of higher tempera- 
ture, to effect complete combination, resulted in a 
comparatively unreactive product. In the prepara- 
tion of calcium hydroxide solutions, calcium car- 
bonate was converted to the oxide, slaked, shaken 
with water, and the mixture allowed to stand until 
it had settled clear. 

When pulverized monocalcium aluminate is shaken 
with water it dissolves in amounts up to 2 g per liter 
or more, forming a solution slightly richer in lime than 








the parent solid, because of hydrolysis [10]. The 
solution thus formed is highly supersaturated with 
respect to certain hydrated compounds, but the initial 
rate of precipitation is slow, so that if the mixture is 
filtered after being shaken 30 to 60 minutes, the re- 
sulting filtrate will remain clear for several minutes 
or even hours. The same result can be attained 
somewhat more conveniently by substituting an 
an aluminous cement for the monocalcium aluminate. 
In the present study, however, the idea of using 
aluminous cement was abandoned when it was found 
that the cement at hand produced solutions contain- 
ing an appreciable amount of sulfur as sulfide. 

Supersaturated calcium aluminate solutions were 
prepared as described above, cooled in the refriger- 
ator, and then mixed with varying proportions of 
previously cooled calcium hydroxide solution in 
l-liter polyethylene bottles. These were stored in 
the cold-cabinet, shaken at intervals, and examined 
from time to time. 

By similar methods the various hydrated calcium 
aluminates were prepared, filtered off, and used in 
subsequent experiments on solubility in water and in 
various concentrations of calcium hydroxide. Thus 
equilibrium was approached on the one hand from 
supersaturation, and on the other, from under- 
saturation. Additional solubility experiments were 
performed with 3CaO-Al,O;-6H,O and 4Ca0-3 Al,O,- 
3H,O0. These compounds were prepared hydro- 
thermally, the first by heating a mixture of Ca(OH), 
and gibbsite at 200° in the presence of water, the 
second similarly from Ca(OH), and CaQO-Al,O; at 
250°. 

Two methods of sampling the solutions for analysis 
were employed. If the solution was clear, the sample 
was taken directly with a pipet, usually 50-ml. If 
any turbidity was noted, the liquid was first filtered 
with the aid of a suction flask and a fritted glass 
crucible, fitted with a stopper and a tube that ex- 
tended into the reaction vessel. Contamination with 
carbon dioxide thus was held to a minimum. (Fil- 
tration of lime solutions through fritted glass intro- 
duces a small error, apparently because of adsorption 
of lime on the glass. This error was disregarded, as in 
no case did it appreciably affect the position of the 
equilibrium curves. ) 

Aluminum and calcium oxides were determined by 
conventional methods. Because of the low concen- 
trations, especially in the case of Al,O,, the errors 
were relatively high in the analyses of 50-ml samples; 
however, the final analysis in each case was done on a 
sample of 100 to 400 ml. 

When portions of the solid were filtered off for 
analysis or examination they were not washed with 
water. Dissolution of the precipitate was thus 
avoided, and the minute amount of dissolved matter 
adhering to the precipitate was considered negligible. 
Usually the precipitates were dried in a desiccator 
containing calcium chloride and a CO, absorbent. 
More rapid drying, when desired, was accomplished 
by washing with alcohol and ether. 

Precise determination of the optical properties of 
the hydrated calcium aluminate usually is difficult 
because of the minute size of the crystals, especially 





in the case of the plate-form phases. The crystals 
formed at 1° C were much smaller than those ob- 
tained under similar conditions at 25° C, so that 
identification by microscopic means often was impos- 
sible. X-ray diffraction, however, proved to be g 


very useful method for identification, even when 


several phases were present. The apparatus used 
was a Geiger-counter diffractometer, with copper 
Ka radiation. The platy phases ordinarily were 
strongly oriented, and gave very strong reflections 
at the lower angles. With certain precipitates cop. 
taining 2CaQO-Al,O;-8H,O it was found essential to 
obtain the X-ray pattern while the material was stil] 
damp, 
drying. 


2.2. Solid Phases in the System CaO-Al,O,-H,O 


a. Hydrated Alumina 


Various investigators [5, 6, 7, 10] have found gibb- 


site, Al,O;-3H,O, to be the stable form of hydrated 


alumina in this system at temperatures from 20° to 
150° C. In the present study its stability at 1° C 
is assumed rather than proved. Actually, as sug. 
gested by D’Ans and Eick [6], the true stable phase 
may be anhydrous Al,O; rather than gibbsite, but 
if so, the transformation does not occur even in 
geologic time. In the present study, gibbsite was 
identified in the precipitates from certain reaction 
mixtures, for example, No. 1 in table 3. In other 
cases, hydrated alumina in some form was believed 
present, even though the X-ray patterns gave no 
evidence of it. Occasionally, after filtration of a 
mixture, the filtrate would contain a little very fine 


suspended material, not readily soluble in dilute HCl. | 


This was assumed to be alumina since all the other 
expected solid phases are readily soluble in acid, 
Again, the examination of one exceedingly finely 
divided precipitate (No. 1, table 2) under the electron 
microscope revealed aggregates of fibers believed to 
be alumina gel. 


b. Calcium Hydroxide 


This is a stable phase in the system at room tem- 
perature [5]. The solubility of coarsely crystalline 
Ca(OH), in pure water was reported by Bassett [11] 
to be 0.130 g of CaO per 100-g¢ solution at 0° C, 
and 0.128 at 5°. He also reported values of 0.148 
and 0.151 for finely crystalline Ca(OH), at 0° C. 
The marked effect of particle size on the solubility of 
Ca(OH), at somewhat higher temperatures has been 
critically studied by Hedin [12]. In the present 
study no attempt was made to fix the solubility 
precisely, since the particle size was not subject to 
control. The value 1.30 g of CaO per liter at 1° has 
been accepted on the basis of Bassett’s data. 


c. 3CaO-Al,O;-6H,O 


The isometric crystalline compound of composition 
3CaO-Al,O;-6H,.O has been shown to be a stable 
phase in the system CaQO-Al,O;-H,O at temperatures 
from 20° to 200° C [4, 5,6, 7,8]. As will be brought 
out later, it probably is not a stable phase at}1° ©. 
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d. CaO.Al,O;-10H,O 


Hydrated monocalcium aluminate was described, 
and its composition discussed, in an earlier paper 


(13]. It exists as a metastable phase in the system 
at 1°. It was obtained as crystals ranging from 


about 254 down to submicroscopic size in different 


preparations. The larger crystals were well-devel- 
oped hexagonal prisms with flattened terminal 
pyramids. 

e. 2CaO.Al,O;-8H,O 
Dicaleium aluminate hydrate, crystallizing in 


hexagonal plates, is a well-established metastable 
phase in the system at room temperature, and, with 
decreasing stability, even as high as 90° C [4, 5, 6, 7]. 
Whether the crystals really are hexagonal or of lower 
symmetry is not definitely known [1]. Despite this 
uncertainty, the phase is generally referred to as 
hexagonal, and this usage will be followed herein. 
MecMurdie [5] indexed the X-ray lines on the basis 
of hexagonal symmetry, and his indices of 001 and 
002 for the first two reflections will be adopted in 
the following discussion. Although the formula 
commonly adopted shows 8 molecules of water, the 
water of hydration is known to vary somewhat with 
the degree of drying. According to Salmoni [14] 
the fully hydrated compound may contain 11 mole- 
cules of water, but this can be reduced to 7 by 
drying the compound over CaCl, and to 5 by heating 
at 100°, without essential change in the crystal 
structure. Further drying is said to alter the 
structure irreversibly. 

Although the dehydration of dicalcium aluminate 
hydrate was not systematically investigated in the 
present work, certain data were obtained that appear 
to necessitate a slight modification of the picture 
presented above. Drying of the precipitate from 
the moist condition to the octahydrate stage (over 
a saturated solution of NH,Cl) resulted in only a 
very minor change in crystal-cell size. Thus the 
largest d-spacing (001) dropped from 10.9 or 10.8 
A down to 10.6. The same effect has been noted by 
others [5]. There was evidence that this continues 
on further drying, down to doo,;=10.4 A (approxi- 
mately), but beyond that point an abrupt change 
was observed. A specimen dried to constant weight 
over CaCl,, and having an approximate composition 
2CaO-Al,O;-6H,O, gave an X-ray pattern indicating 
spacings distinctly different from those of the octa- 
hydrate. The largest spacing was 8.8 A. The only 
change in appearance under the microscope was a 
sight inerease in refractive indices. Numerous 
preparations showed both hydrates coexisting in 
varying proportions. 

An X-ray pattern of the lower hydrate is given in 
table 1, along with one of 2CaO-Al,O;-8H.O. The 
octahydrate pattern is included merely for compari- 
son; the relative intensities given indicate a high 
degree of orientation of the mounted sample. No 
doubt this factor likewise affected the relative 
intensities in the second pattern, though to a smaller 
extent. Certain similarities between the two pat- 
terns are evident. On the basis of MeMurdie’s 


indexing, the pattern of 2CaQ-Al,0O,-8H,O [5], 
together with the fact that the crystals do not lose 
their platy habit on dehydration to the 6H,0 stage, 
it appears likely that the lattice remains conceal 
unchanged except for shortening along the c-axis. 
This suggests that the water thus removed may 
have been present in loosely bound layers within the 
structure. The change from 8 to 6H,O is completely 
reversible and occurs rather rapidly. Thus certain 
samples actually went from the higher to the lower 
hydrate while on the X-ray diffractometer, during a 
period of less than 1 hour. 

The amount of water in the lower hydrate is 
known only approximately. The material used for 
the X-ray pattern (table 1) had the composition 
2.10CaO-Al,0;-6.14H,O. Another preparation had 
the composition 1.95CaQ-Al,O,-5.22H,O, but this 
was shown to contain a small percentage of the cubic 
3CaO-Al,0,6-H,O0 (which would also necessitate the 
presence of some form of Al,O;); the latter figure is 
therefore less reliable. The value, 6H,O, appears 
most probable from the data at hand. 














TaBLe 1. X-ray diffraction lines of 2CaO-Al,03-8H,0 and 
2CaO.-Al,O3-6H,O 
2Ca0-AloO3-8H20 2Ca0-AloO3-6H20 
Spacing Relative Spacing | Relative 
intensity intensity 
A A 
10. 6 100 8.8 100 
5. 27 8O 4. 80 13 
3.49 20 4. 36 65 
2. 88 14 3. 80 5 
2. 67 4 3.60 2 
2. 62 7 2. 86 22 
2. 53 5 2.72 ll 
2.49 7 2. 57 8 
2.44 5 2. 47 12 
2. 39 3 2. 39 8 
2. 22 3 2. 32 8 
2.12 2 2. 29 3 
2.03 2 2. 23 9 
1, 954 2 2.14 2 
1. 935 3 2. 04 ll 
1. 740 2 1, 802 4 
1. 662 7 1. 742 4 
1. 528 2 1. 653 8 
1. 440 3 1. 623 7 
1, 424 1 1, 441 3 
1. 352 4 
1. 292 4 
1, 289 3 
1. 260 2 








f. 4CaO.Al,O;-13H,O 


Tetracalcium aluminate hydrate crystallizes in 
hexagonal plates, almost indistinguishable under the 
microscope from those of dicalcium aluminate 
hydrate. Like the latter compound its existence is 
well established, as a result of numerous investiga- 
tions [1, 4, 5, 6,9]. The comments given above con- 
cerning the symmetry of the dicalcium aluminate 
hydrate apply equally to the tetracalcium. Just as 
with the dicalcium compound the amount of water 
of crystallization is in some doubt. Most of the 
compositions reported in the literature range from 
12 to 14H,0.[5]. The formula with 13H,O has 
been generally adopted, but certain work of the 





Building Research Station in England as reported 
by Lea [3] indicates that the fully hydrated com- 
yound may contain 19H,O. Some investigators, 
fed by Assarsson [9], believe there are two forms, 
designated a and 8, with the same degree of hydra- 
tion but slightly different degrees of stability in the 
system; also perhaps a third form with one 
molecule of water. 

In the present study only one form of tetracaleium 
aluminate hydrate was produced, judging from the 
X-ray diffraction patterns. (Certain apparently 
contradictory evidence will be discussed below.) No 
special study was made of the amount of water of 
hydration. <A slight shift in the X-ray diffraction 
lines with increased drying was noted. 


less 


g. 3CaO.Al,O,;-CaCO;-11H,O 


This compound obviously falls outside the system 
CaQ-Al,O,;-H,O; nevertheless, it may appropriately 
be considered here. During the study of X-ray 
diffraction patterns of numerous preparations in the 
high-lime region of the system, a series of lines was 
frequently noted very close to those of 4CaQ-Al,O,- 
13H,O. The maximum spacing averaged about 7.7 
A, as compared with 8.2 for the tetracalcium com- 
pound. At first the new phase was tentatively con- 
sidered to be the 8-form of this compound, but when 
it later appeared in the low-lime region of the system 
as well, another explanation was sought. To ascer- 
tain whether it might be a carboaluminate such as 
that reported by Bessey [4], an attempt was mad> 
to prepare such a compound. From a mixture of 
calcium aluminate, calcium hydroxide, and ammo- 
nium carbonate solutions a precipitate was obtained, 
consisting of minute platy crystals having indices 
and birefringence somewhat higher than those of 
4CaO.-Al,O;-13H,O. The X-ray pattern agreed per- 
fectly with the series of lines in question. Although 
the amount of material was insufficient for a precise 
analysis, the composition was reasonably close to the 
formula 3CaQO-Al,O,-11H,O suggested by Bessey. 
Further confirmation was deemed unnecessary, 
because just at that time a publication by Turriziani 
and Schippa [15] appeared describing the preparation 
of a carboaluminate of composition 3CaQ-Al,O,- 
CaCQO,-10.SH,0 and reporting an X-ray pattern ir 
good agreement with those under consideration. 

In only one of the preparations (No. 13, table 2) 
were the carboaluminate crystals of sufficient size to 
permit a determination of their indices of refraction, 
and even in this it was difficult because of the mask- 
ing effect of the other phases. The high and low 
indices were 1.554 and 1.527, respectively. 

The question as to the source of the CO, naturally 
arises. Although there was undoubtedly a little 
entry of atmospheric CO, when the bottles were 
opened for removal of samples, it seems likely that a 
more important factor was the passage of CO, 
through the walls of the polyethylene bottles. Poly- 
ethylene is known to be somewhat permeable to 
CO, [16], although the effect apparently was rather 
In 


slow under the conditions of the present study. 





most cases the amount was barely 
shown by the X-ray patterns. In 


enough to be 
DY a few it Was 
much higher, possibly indicating that the bottle egy 
did not form a perfect seal. 


h. Other Aluminates 


Certain other hydrated calcium aluminates tha} 
have been reported by various investigators may also 
be mentioned. No evidence of a hexagonal trical. 
cium aluminate hydrate was found in this study 
Compositions having approximately this ratio jp. 
variably gave X-ray patterns indicating mixtures of 
di- and tetracalcium aluminate hydrates. This nega. 
tive evidence supports the view of Wells, Clarke, and 
MeMurdie [5] that the alleged tricalcium compound 
does not occur. 

Another tricalcium aluminate hydrate, in acicular 
form, and assigned the formula 3CaQ-Al,O,-18 to 
21H,O, has been reported by several investigators 
[2]. It was not found in the present study. How- 
ever, the author has previously prepared a compound 
of this general description and approximate composi- 
tion. It was found only in aqueous extracts of an 


TA 


} 


aluminous cement, and then only as a transitory ' 


phase. The solutions contained a trace of sulfide. 
which may or may not be significant. 

The more basic aluminates that have been re- 
ported were not considered in this study, since they 
probably do not possess a field even of metastability 
in the system. There remains the compound 4Ca0. 
3A1,0;-3H,0, which is formed under hydrothermal 
conditions but is unstable below 200° C. Although 
it has no place in the system at 1° C, a pair of ex- 
periments were conducted to study the course of its 
reactions. 


2.3. Precipitation Experiments 


A series of 10 supersaturated calcium aluminate 
solutions were prepared as described above, by dilut- 
yng a master solution with various amounts of satur- 
ated lime water. They were then stored in the cold- 
cabinet. These are designated by numbers 1 to 10 
in table 2. From time to time, samples of the solu- 
tions were taken for analysis, and small portions of 
the precipitates for examination. Data obtained at 
intermediate periods have been omitted from the 
table, except W hen needed to show the course of the 
reaction. The figures given for initial concentrations 
(columns 2 and 3) were calculated from the measured 
volumes of the analyzed master solutions. Concen- 
trations at later periods (columns 5 and 6) were ob- 
tained by analysis, and the composition of the pre- 
cipitate (column 7) was calculated from the difference. 
In the last column, the solid phase present in greatest 
amount is given first. 

The change in concentration during precipitation 
may be visualized better by reference to figure 1, in 
which the concentration data for solutions 1 to 10, 
table 2, are plotted on rectangular coordinates. _Ini- 
tial concentrations are represented by the row of; 
circles at the right. The lines extending downward 
to the left represent the changing concentrations. 
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TaBLE 2. Data showing changes in concentration and solid phases present during precipitation from calcium aluminate solutions 





at 1°C 
Initial concentration Concentration after (Molar ratio 
reaction period CaO0:AleOs,, 
Experiment Reaction period - in precip- Solid phases present * 
itate (cal 
AlO CaO ALO CaO culated) 
g/liter g/liter g/liter g/liter 
tt = . 1. 266 0). 828 5 ssnapessumeneiaa 
1. 278 0). 833 ff as a 0. 144 . 190 1. 03 Ca0O-AlyOs-10H20; probably a little AlsOs hyd.; trace of 
1..--- ‘ ™ cx 2CaO0-AhOs-8H20. 
13 mo 112 . 152 1.06 Same. 
{1 day i 436, . 393 1. 24 siete _ ria cee 
2 me 1. 1% 879 11 days : . 162 . 316 1.04 CaO-AlyO3-10H20; trace of 2CaO-AlgOs-8H20. 
(75 days . 104 . 300 1.01 Same, plus trace of 3CaO-AlyO3-CaCO3-11H20. 
| 
1 day , . 572 412 2.07 ee 
7 days ; 422 . 402 1. 58 omnia wie ee ee | 
1.022 gor 11 days. . 332 . 330 1. 57 2CaO-AlyOs-8H;:0; traces of CaO-AlpO;10H,0, 4Ca0O- 
- ‘ AloO3-13H 0, and 3CaO-AlO3-CaC O3-11H20. 
13 mo 112 . B88 1. 07 2CaO-AlgO3-8H20; CaO-AlpOs-10H20; 3CaO-AlyOs.- 
CaCO; 11H,0. 
poe ind) 404 2. 06 _ wae se a 
0. 89 170 11 days . 314 32u 2.01 2CaO-AloOs-8H20; trace of 3CaO-AlOs-CaCO;3-11H20. 
1 mo . 316 333 2.00 _ ce 7 Le Se 
15 mo ‘ . 120 3u4 1.3 2CaO-AlyOs-8H2O; CaO-AlpOs-10H20; 3Ca0-AlLOs-- 
CaCO;y11H,0. 
{1 day . 230 412 2.04 ‘ enoeusese . : 
767 1.014 11 days .174 345 2.05 2CaO-AlpOs-8H)0; trace of 3CaO-AhO3-CaCOs-11H,0. 
9 mo 203 . a 2.09 Same. 
{1 day 168 . 512 2. 11 ee SE ee -- | 
. 639 1. 058 +11 days . 110 442 2.12 2CaO0-AlhOs-8H20; 3CaO-AbOs-CaC Oz 11H,0. | 
19 mo 085 . 384 2.21 Same. 
j} day . 106 . 43 2. 54 : eebneane won 
7 ll 1. 108 411 days . 096 . 556 2.42 2CaO-AleOs-8H20; 4CaO-AlpOs- 13H20. 
11 mo . . 078 . 3 2. 54 Same. 
{1 day . : . 140 670 3. 62 . ‘ a eaecees oo | 
. ISS 1. 154 +11 days seicinaiol OYS . 570 3.73 41CaO-AloOs-13H20; 2CaO-AlO;-8H)20, 
11 mo_. . .072 . 499 3. 83 Same, plus 3CaOQ-AlhO3-CaCO;3-11H;0. | 
1 day . . 090 . 870 3. 52 - aaiiaia — EE Oe ee 
256 1.191 11 days . 026 . 726 3. 67 4CaO- AloOs-13H20; 2CaO-AlO3-8H20. 
{11 mo-. 016 . 695 3. 76 4Ca0-AlO3-13H20; 3CaO-AlO3-CaCO;-11H20. 
ji day. i . 046 1.111 2. 99 a at eee a | 
it 128 1, 246 11 days heal . 004 1.018 3. 34 4CaO-AlgO3-13H20,. | 
{9 mo ; . 004 0. 963 4.15 Same. 
\3 hr. 7 570 . 908 1, 81 2CaO- AloOs-8H20; CaO-AloOs-10H20. | 
1] : 1.310 0. 970 1 day. _ 208 . 322 1.07 CaO-AlO3-10H30. 
7 mo.. ‘ ov2 . 260 1. 06 Same 
. hr . 328 . 366 2.16 2Ca0-AlO;3-8H20; trace of 3CaO-AlO;3-CaCO;y-11H;20. 
12 0. 916 1. 066 2 days , . 312 . 364 2.12 Same. 
“ : 5 mo. . 240 . 332 1.97 PSS - I E 
lio |” ae . 153 . 372 1. 65 2CaQO-AlO;-8H20; CaO-AleOs-10H20. 
1 mo....-- . 449 .317 Se Pee eee ee chiaidiotan 
13 166 0. 397 > one....- me 364 312 0. 27 4 lis bi PSS eae ee = fie 
7 sp esecee . 152 . 229 . 57 Gibbsite; CaO-AlO;-10H20; 3CaO-AlO3-CaCO3-11H20. 
ee . 100 . 209 . 54 Same. 
|! a . 386 . 269 1. 32 Gibbsite. 
4b. . 466 . 327 =a . 320 . 266 0. 76 . OR ee 
{14 mo-- icine . O85 . 200 . 61 Gibbsite; CaO-AlO3-10H20; 3CaO-AlO3-CaCO3-11H20. 
5 310 610 fl mo... — .104 . 410 1.77 Gibbsite; few hexagonal plates. 
at . {11 mo_...--- - . 103 . 399 1.85 Gibbsite; 2CaO-AleO3-8H20; 3CaO-AlgO3-CaCO3-11H20. 
16 155 893 fl mo : — . 038 . 691 3. 12 Gibbsite; hexagonal plates. 
e i . . 11 mo.. . 022 . 629 3. 61 Gibbsite; 4CaO-AlpOs-13H20; 3CaO-AlOs-CaC O3-11H20. 








* As used here, “‘trace’’ is equivalent to ‘‘barely sufficient to be detected on the X-ray diffraction pattern.” » Gibbsite added to solution. 


The calculated molar ratio of CaO to Al,O; (C/A) in 
the precipitate in the early stages is indicated in each 
Case, 

Solution 1 was essentially unchanged after 1 day, 
but subsequently a vary voluminous precipitate ap- 
peared. The X-ray diffraction pattern showed it to 
be chiefly CaO-Al,O,-10H.O, with a detectable amount 
of 2CaO-Al,O,-8H,O. Under the electron microscope 
the predominant phase was observed to consist of 
slender prismatic crystals, less than 0.1 g in thick- 
ness. A few much larger hexagonal crystals, pre- 
sumably of 2CaO-Al,O;-8H,O, were noted. There 


were also a few masses of felted fibers, finer than the 
prismatic crystals. It could not be determined 
whether these were of the same composition or 
whether they were a different phase, possibly amor- 
phous hydrated alumina. Precipitation was essen- 
tially complete after 11 days, and no change in ap- 
pearance was observed after 13 months. 

In the other members of the series, precipitation 
was more rapid, increasingly so as the CaQ/AI,O, 
ratio increased. No. 2 yielded a precipitate con- 
sisting largely of CaO-Al,O,;-10H,O, appearing as ro- 
settes of prismatic crystals much coarser than those 
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Fiaure 1. Change in concentration of calcium aluminate 
solutions resulting from precipitation at 1° C. 


C/A=molar ratio, CaO/AhO;. Numbers refer to experiments in table 1, 


in No. 1. A detectable amount of the dicalcium 
compound was also present. After 75 days a steady 
state apparently had been reached, with no change 
except for the appearance of a detectable amount of 
the carboaluminate. At that time the precipitate 
was filtered off for use in further experiments. 

The change in slope of the concentration curve for 
solution 2 suggested that the initial precipitate may 
have been richer in lime than was indicated by the 
analysis after 1 day. Another solution of compara- 
ble concentration therefore was prepared and ex- 
amined at shorter intervals. ‘This is designated No. 
11 in table 2; it has been omitted from figure 1 to 
avoid confusion. The initial precipitate was ob- 
served to be 2CaQ-Al,O,-8H,O, and the ratio of CaO 
to Al,O, after 3 hr was nearer 2 than 1. However, 
after 1 day, only CaQ-Al,O,;-10H,O was observed. 
Thus the initially formed dicalctum compound must 
have redissolved, in agreement with the observed 
change in slope of the concentration curve. 

Solutions 3, 4, 5, and 6 all precipitated 2CaO. 
Al,O;-8H,O as the primary phase. The calculated 
CaO/Al,O; is consistently slightly higher than 2. 
Any attempted explanation would be speculative, 
but it is clear that the inclusion of a small amount of 
tetracalcium aluminate or carboaluminate would 
tend to raise the ratio. The precipitate in No. 3 
underwent partial conversion to CaQ-Al,O,;-10H,O 
within 7 days, as indicated by the change in CaO 
Al,O; ratio, later confirmed by X-ray diffraction. 
The same reaction occurred in No. 4, but in this case 
it did not start until after 1 month. 

Experiment 12, table 2 (omitted from fig. 
intermediate between Nos. 4 and 


1) was 
5 in initial CaO 





Al,O; ratio. It behaved like No. 4, except that th 
appearance of the monocalcium aluminate was gtijj 
longer delayed. Nos. 5 and 6 showed none of the 
latter compound even after 9 months. 

Solutions 7, 8, and 9 precipitated mixtures of 
2CaO.-Al,O;-8H,O and 4CaO-Al,O;-13H,O, with very 
little change after 11 days except for the appearance 
of small amounts of carboaluminate. The mola 
CaQO/Al,O; ratio in the precipitate was between 2 ang 
4. The initial ratio in the precipitate from solution 
10 was calculated to be 2.99. This is suspected of 
being in error, as the low-alumina concentration jp 
this region greatly magnifies the effect of smalj 
analytical errors. The X-ray patterns of the pre. 
cipitate showed only 4CaO.-Al,O;-13H,0. 

To recapitulate: No. 1 precipitated only Ca0. 
Al,O;-10H,O; Nos. 11, 2, 3, 4, and 12 gave 2Ca@. 
Al,O;-8H,0 first, which was replaced by CaO.Al,0, 
10H,O, more slowly and less completely as the pro- 
portion of CaO in solution increased; Nos. 5 and § 
gave only 2CaQO-Al,O;-8H,O; Nos. 7, 8, and 9 yielded 
mixtures of the dicalcium and tetracalcium aluminate 
hydrates; and No. 10 precipitated only the tetra. 
calcium compound. 

Reaction mixtures 13 to 16 constitute a series ip 
which the initial concentration of both Al,O, and 
CaO was much lower than in the foregoing group 
One gram of solid gibbsite was added in each cass, 
in an attempt to induce precipitation of this stable 
phase instead of the metastable aluminates. As in- 
dicated by the results in table 2, this attempt was 
unsuccessful. Aside from supplying additional data, 
this group is of interest chiefly in that the crystalline 
phases formed very slowly, yielding crystals of ade. 
quate size for optical study. 


2.4. Solution Experiments 


It is evident from the foregoing that the approae 
to equilibrium from supersaturation at 1° C is a ver 
slow process. Attempts were therefore made to ap 
proach it from the opposite direction, that is, by de 
termining the solubilities of the various solid phases 
in different concentrations of calcium hydroxide. In 
each experiment the solid phase was added in moé- 
erate excess to pure water or to lime water of known 
concentration. The procedures for storage, analysis 
and examination were the same as for the precipite- 
tion experiments. 


a. Hydrates Formed by Precipitation 


It will be convenient to consider first, as a group 
the three ternary compounds found in the precipite 
tion experiments; namely, CaO-Al,O,-10H,O, 2Ca0 
Al,O;-8H,O, and 4CaQO-Al,O;-13H,O. Data for thi 
group are given in table 3, concerning which a fev 
explanatory notes are in order. Final concentrations 
were determined by analysis, whereas initial concel 
trations were calculated. In certain cases the solid 
were introduced in the form of a suspension in th 
solution in which they were prepared. This accounts 
for the fact that not all the experiments started from 
zero Al.O,; concentration. 
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TABLE 3. 


Solubilities of CaO-Al,O3-10H,0, 2CaO-Al,0;-8H,0, and 4CaQO-Al,O3-13H,0 in calcium hydroxide solutions 








———— 
<< 





























Initial concentration Dura- Final concentration | 
Experi- Solid phases initially present tion of | ___| Solid phases present at end of experiment 
ment contact 
Al:Os CaO | | AlzOs CaO | 
giliter | g/liter Months g/liter g/liter 
tcl { 0.000 | 0. 000 12 0. 070 0.169 | CaO-Al,O3-10H20; calcite, gibbsite. 
2 __...\}CaO-AlzOs-10H20. 056 | 21 2 | 119 222, | CaO-Al:Os-10H20; carboaluminate. 
eee | | . 056 283 | 2 .118 . 290 CaO-AlsO3-10H20; very little carbo- 
aluminate. 
| | 
4._..-.-| CaO-Al203-10H20+2Ca0-Al203-8H20-..--. -000 | . 000 5 . 130 . 368 Ca0-Al:03-10H20; 2Ca0-Alz:03-8H20. 
a ae f .000 000 | 65 | .61 339 ' 
5..-----| 19Ca-AlzOs 8H:0-.. --} ‘ooo | . wai 6 | * 182 * 324 }2Ca0-AlsOs-8H:0; carboaluminate. 
hail f -0o1 | 35 4 - 081 .419 | 2CaO-AlpO3-8H20; CaO-AlsO3-10H20. 
a _...-[}8Ca0-AlsOs-SHisO ®........-..e-nccrcccccccce 001 | 52 ort . 049 .5387 | 2CaO-AlsOs-8H:0. 
— | 001 69 4 | 024 .630 | 2Ca0-Al,:0s-8H20; 4Ca0-AlsOs-13H20. 
10... -] '{ . 005 514 | 4 . 040 . 550 
. " ‘ -005 | 743 CO | ° . 69% ‘ 
I1.----- 2030. Al0r-8H20-+4Ca0- Al:Os-13H30.... —. = | ss | ‘aa Ga) :|{2Ca0-AlaOs-8H20; 4Ca0-AlsOx-13H20. 
eee - 003 1.143 | 4 | -001 1. 098 
| | 
4. . 000 0.000 | 1 | 020 | O34 | 
. 000 -000 | 1 028 | . 396 
a .000 | 000 | 1 026 | . 442 
— . 000 000 | 1 024 | . 489 
13... . 000 . 000 1 . 026 . 565 
O....--| sad. AlsOp 18HsO ©... .............------- “Oe naa : as 610 | 4Ca0- Al:05-13H20. 
Weeks 
0... .000 4 | 6 009 | .758 
ae . 000 65 | 6 . 008 | . 868 
eee -000 | . 80 | 6 . 001 . 971 
3... . 000 1.00 | 6 002 | 1. 103 
| 








s Contained a minor amount of CaO-Al:03-10H20. 
> Contained Ca(OH)s. 


Great difficulty was experienced in obtaining 
single-phase preparations for starting materials. 
Thus the dicalcium aluminate hydrate used in ex- 
periments 7 to 9 contained some monocalcium alu- 
minate hydrate. This phase persisted in No. 7 
but dissolved in the two more basic solutions. 
Tetracalcium aluminate hydrate usually was ob- 
tained as a mixture with the dicalctum compound. 
(Such a mixture was used in experiments 10 to 13.) 
In order to avoid the less basic aluminate it was 
found necessary to have free calcium hydroxide 
present. The latter contaminant dissolved com- 
pletely during the solution experiments, and in 
experiments 14 to 19 the equilibrium concentration 
of Ca(OH), was varied merely by varying the 
amount of solid originally used. 

In this group of experiments a steady state was 
attained in a week or two, no doubt because of the 
extreme fineness of the crystals. The final concen- 
trations given in table 3 are represented in figure 
2 by filled circles, and the tentative solubility curves 
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tCaQ-Al,O;-13H.0 in calcium hydroxide solu- 


Curves are based on data from table 3 (filled circles). 
Circles) are included for comparison. 


Data from table 2 (open 


In experiments 14 to 19, concentration of Ca(OH): in solution was varied by using increasing amounts of the solid preparation. 


are drawn with reference to these points. The open 
circles represent the final data from the precipitation 
experiments (table 2), and are included here to facili- 
tate comparison of the two sets of data. 

Curve ABC, figure 2, represents approximately 
the solubility values for 2CaOQ-Al,0O,;-8H,O. The 
corresponding concentration values from table 2 
lie somewhat above the curve, but not far from it. 
Curve CD fits equally well the data for solubility 
of 4CaQO.Al,O,-13H.O or for a mixture of the tetra- 
and dicalcium compounds, including the values 
from both table 2 and table 3. The horizontal 
line FC fits the data for 4CaO-Al,0O,-13H,O in the 
region below 0.6 g of CaO per liter. (There is, of 
course, no counterpart of this in table 2.) Line 
EB represents roughly the solubility of CaO-Al,O,- 
10H.O. Here the values from table 2 (open circles) 
all fall below the curve, as does the value for experi- 
ment 1 in table 3. 

There should be an invariant point (metastable) 
at the intersection (B) of the curves for the mono- 
and dicalcium aluminate hydrates. The ee 
obtained in experiments 4 and 7 in table 3, and 3 
and 4 in table 2, ended essentially as mixtures of 
CaO.-Al,O,-10H,O and 2CaQ-Al,0,;-8H,O. Three of 
these had compositions very close to B, probably 
within the limits of accuracy of the procedures, but 
the fourth was appreciably removed. This may 
reasonably be attributed to lack of attainment of 
equilibrium. 

Point C does not represent an invariant point. 
The compound 4Cx0-ALO,.13H,0 appears to be in 
equilibrium all along the line FCD. The fact that 
2C'aQ-Al,O;-8H,0 coexists, seemingly indefinitely, 
with 4CaO-Al,O,-13H,0 between C and D (and in 





fact, beyond D to the limit of Ca(OH), solubility), 
points to extreme slowness of attainment of equilib- 
rium in the area. If this view is correct, there 
should be two solubility curves between C and D, 
one for each hydrate, but they must be practically 
coincident. 


b. 3CaO.A!.O;-6H,O 


From earlier cai [5, 6, 7] it is known 
that the isometric 3CaQ-Al,O,-6H,O exists as a 
stable phase in the system CaO-Al,O;-H,O from 20° 
to about 200° C. The experiments summarized in 
table 4 comprise a series in which the solubility of 
this compound in water and in calcium hydroxide 
at 1° C was sought. The results were highly erratic. 
The most interesting fact, however, is that the 
original solid phase was largely replaced by 2CaO. 
Al,Oz;-8H,O in experiments 3 and 4 and partially 
by 4CaQ-Al,O;-13H,O in Nos. 7, 9, and 11. Inas- 
much as these hexagonal phases were known to be 
absent from the original solid, it must be concluded 
that 3CaQ-Al,O,-6H.O is not stable in these areas. 
The final concentrations here fall close to the solu- 
bility curves for the respective hexagonal phases. 

The final concentration data in table 4 are plotted 
in figure 3, superimposed upon the solubility curves 
of the hexagonal aluminates, borrowed from figure 
2. It is evident that the plotted points, taken as 
a single group, cannot be connected by any reason- 


able curve. The experiments may, however, be 
separated into two groups: (a) those in which 
the residue was partially transformed to di- and 


tetracalcium aluminate hydrates, and (b) those in 
which the 3CaQ-Al,0,-6H,O remained unchanged 
(except for minor formation of carboaluminate). 
Group (a) is represented in figure 3 by filled circles, 
which fall not far from curve ABCD, the solubility 
curve of the di- and tetracalcium aluminates. Group 
(b) is represented by open circles, which appear to 
be scattered at random appreciably below curve 
ABCD. 

It is concluded that 3CaO-Al,0,-6H.O is less stable 
ae 2CaO.-Al,O,-8H,O and 4CaQO-Al,O,-13H.O at 1 
’, but that it goes into solution extremely slowly. 
i or some une xplaine “l reason, certain experuments in 
group (b) failed to reach equilibrium concentrations 
even after as long as 16 months of storage. Others 
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Figure 3. Solubility of 1Q-Al,O;-6H2O in calcium hy- 


droxide solutions at 1° C. 


Open circles represent mixtures in which the solid phase was essentially yp. 
changed; closed circles those in which phase transformations had occurred, 
Solubility curves, taken from figure 2, are included as a frame of reference, 


in group (a) reached concentrations relatively quickly 
(within 1 month for experiments 3 and 4), at whieh 
time they were supersaturated with respect to the 
hexagonal phases, whereupon the latter precipitated 
out. Even here the reaction had not gone to com- 
pletion after prolonged storage. It is conceivable 
that the small amount of carboaluminate appearing 
during the long storage period may have formed 4 
protective film on the crystals. 


c. 4CaO.3Al,0;-3H,O 


Although the 4:3:3 ternary compound is stable 
only above approximately 200° [7], it was con- 
sidered of interest to study its behavior at 1° C. In 
water the compound dissolved slowly, without forma- 
tion of a precipitate, and after 10 months had reached 
the surprisingly bigh concentration of 0.382 g of 
Al,O; and 0.285 g of CaO per liter. By reference to 
figure 2 it appears that such a solution must be highly 
supersaturated with respect to CaO-Al,O;-10H,0. 

The solubility in caletum hydroxide was studied at 
only one concentration, 0.58 g of CaO per liter. The 
reaction proceeded slowly, and involved simultaneous 
dissolution of the 4:3:3 compound and precipitation 
of 2CaQO-Al,O,-8H,O. The final concentration was 
0.092 g of Al,O, and 0.435 g of CaO per liter, which is 
close to curve BC, figure 2. 


Gibbsite (Al.O,-3H,O) 


An attempt was made to determine the solubility 
of gibbsite in a series of lime solutions ranging from 
0.2 g of CaO per liter up to saturation, but the results 
were unsatisfactory. After 4 months at 1° C, the 
AI,O, in solution was barely detectable. In an effort 
to speed up the process, the mixtures were allowed to 


TABLE 4. Solubility of 3CaO-Al,03-6H,O in calcium hydroxide solutions 
Initial con-| Duration Final concentration 
Experiment centration, 0 Solid phases present at end of experiment 
CaO contact 
ALO CaO 
g/liter Months a/liter g/liter 
l 0. 000 16 0.117 0. 235 3CaO- AloOy 6H2O; 3CaO-AlyOsCaC Oy 11H,0 
2 O00 15 134 268 Do 
} 000 9 166) 300 2CaO0-AlpOs-8H20; 3Ca0-AlOs-6H20; 3CaO-AlO;-CaC Ov 1LH,O 
4 O00 Q 164 3AY Do 
5 128 ] 100 338 3CaO-AlpOs-6H20; 3CaO-AloOs-CaC O;3-11H2,0 
6 257 16 028 319 3CaO0-AleOy- 6H2QO; 3CaO-AleOsCaC Oy 1LH20 
7 434 LS 040 515 3CaO0-AlOs6H2O; 4CaO-AloO;-13H,O; very littl 3Ca0-AhO 
‘aC Oy llLH,O 
8 535 14 O10 : 554 3CaO- AloOy 6H2O; 3CaO-AlyOs-CaC Oy 11LH2O 
) 778 14 O18 22 3CaO-AloOs-6H2O; 3CaO-AloOs- 11LH2O; 4CaO0-AlOs 13H2O. 
10 1. 039 14 003 1. 055 3CaO-AlpOs-6H20; 3CaO- AlpOs-CaC O3-11H20 
11 * 1.476 10 OO! 1. 371 Ca(OH)o; 3CaO-AlpOs-6H2O; very little 4CaO-AlyO s-13H,O 


* Contained excess Ca(OH), in suspension 
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stand 6 months at room temperature, after which 
they were again stored at 1° for 10 months. In the 
first solution (0.2 g of CaO per liter), the Al,O; con- 
centration rose from 0.002 to 0.012 g per liter during 
the storage at room temperature, and failed to change 
on subsequent storage at 1° C. This figure approxi- 
mates the solubility in lime solution of the same con- 
centration at 21° C [5] and the true value at 1° would 
be expected to be lower. 

In all the more basic solutions (0.4 g of CaO per 
liter and higher) there was precipitation of 4CaO. 
Al.O,-13H,O, indicating that the gibbsite had dis- 
solved sufficiently to exceed the solubility of this 
aluminate. Final concentrations fell close to curve 
FCD, figure 2. 

Because of the long time involved, further work 
on the solubility of gibbsite in lime solutions at 1° C 
was abandoned. The only conclusions to be drawn 
from the present data are: (1) That gibbsite appears 
stable in lime solutions up to 0.2 g of CaO per liter; 
2) that the equilibrium concentration of Al,O, at 
this point is between 0.002 and 0.012 g per liter; and 
2) that at some higher lime concentration, the solu- 
bility curve of gibbsite intersects that of 4CaQO- 
Al,Og-13H.0. 


2.5. Phase Equilibrium Diagram 


From the foregoing data it is possible to construct 
at least a partial diagram of stable and metastable 
phase equilibria in the system CaQ-Al,O,-H,0O at 1 
C. Such a diagram is presented in figure 4. Cer- 
tain assumptions and qualifications are necessary. 

One assumption is that the carbon dioxide in 
solution and the carboaluminate precipitated were 
without effect. While this cannot be strictly true, 
it is evidently not far from correct. Carbon dioxide 
in very small concentration is known to react with 
lime solutions, forming calcite. In the present work 
the carboaluminate formed in preference to calcite 
in nearly every case, hence the CO, remaining in 
solution must have been even less. Also, the pre- 
cipitated carboaluminate probably is the least 
siluble of the phases present (possibly excepting 
gibbsite), and it would therefore not alter the equi- 
librium concentration. 

The position of the solubility curve of gibbsite, 
GF, is merely an estimate. The position of D, the 
solubility of calcium hydroxide, is based on values in 
the literature [11]. Curve CD, of course, cannot 
meet the horizontal axis at D; rather there must be 
an invariant point nearby, with an Ai,O; concentra- 
tion about 0.001 g per liter. Two curves are 
assumed to exist between this invariant point and C; 
the lower one the stable curve for 4CaO-Al,O,-13H,0O, 
the upper a metastable curve for 2CaO-Al,O,-8H,0. 

Curve ABC is essentially a compromise between 
data obtained respectively from precipitation and 
from solution experiments. The upper end (near 
A) takes into account the inflection points of certain 
lines (3 and 4) in figure 1. Curve EB is somewhat 
uncertain as to position, as there is considerable 
scatter among the points on which it is based. 

With these qualifications and assumptions in 
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mind, the diagram may now be summarized. The 
stable phases are gibbsite, tetracalctum aluminate 
hydrate, and calcium hydroxide. The position of 
point F, the invariant between gibbsite and 4CaO. 
Al,O,-13H,O0, is not known. The invariant between 
4CaO.-Al,O,-13H,0 and Ca(OH)., at D, is estimated 
to be at 0.001 g of Al,O, and 1.30 g of CaO per liter, 
and the solubility of Ca(OH), (from the literature 
{11]) is likewise placed at 1.30 g of CaO per liter. 
The compound 2CaQO-Al,O;-8H,O is in metastabie 
equilibrium with solution along curve ABCD, the 
section CD being practically coincident with the 
curve for 4CaQ-Al,O,-13H,O. Nothing is known as 
to its extent above A. Curve EB represents the 
equilibrium of CaQ-Al,O,-10H,O and solution; its 
leftward extent is unknown. The point B is a meta- 
stable invariant between 2CaQ-Al,O,-8H,O, CaO. 
Al,O,-10H.O, and solution. 

Metastable solutions to the right of AB tend to 
precipitate 2CaQO-Al,O,-8H,O until the concentration 
reaches AB, whereupon CaQ-Al,O;-10H,O is formed, 
and the concentration shifts toward B (see table 2, 
experiments 3, 4, 11). Solutions richer in lime (for 
example, No. 7 in table 2) precipitate both di- and 
tetracalcium aluminate hydrates until the concen- 
tration reaches BC. At this point, presumably, 
4CaQO-Al,O,-13H,O continues to precipitate out, but 
as long as there is any 2CaQ-Al,O,-8H,O left, the 
latter will continue to control the concentration of 
the solution, which thus will still lie on BC. In 
time the dicaletum compound should be used up, and 
the concentration fall to curve FC, but this stage 
was not reached in any of the precipitation experi- 
ments described above. 

The area under GFCD represents stable solutions. 
Solutions above GF are supersaturated with respect 
to gibbsite, but the rate of precipitation is extremely 
slow. Solutions in the general area ABE can pre- 
cipitate CaO-Al,O,-10H,0. 


3. Discussion 


3.1. Relation to Earlier Data 


After reviewing the literature on the system CaQ- 
Al,O,-H,O, one is impressed by the wide variance, 
not so much in the actual data as in the interpreta- 
tion of the data. Some of the results obtained in the 
present study may help to explain some of the contra- 
dietory conclusions. 

Several investigators of the calcium aluminates 
have stressed the necessity for rigid exclusion of 
carbon dioxide; yet it is probable that complete 
absence of CO, has seldom, if ever, been attained. 
Even freshly boiled distilled water retains some CO,, 
and the rubber stoppers commonly used for closures 
of the reaction flasks are by no means impermeable 
to atmospheric CO,. Formation of calcite has, in 
fact, been reported in some cases [10]. This phase, 
if present, is conspicuous under the microscope by 
reason of its high birefringence and distinctive 
crystal shape. It is only natural that the absence 
of calcite might lead the investigator to assume the 
absence of CO, from his reaction vessel. On the 








other hand, if any of the carboaluminate, 3CaQO.- 
Al,O,;-CaCO,-11H,O, were present, it might readily 
be mistaken for 4CaQO-Al,0,;-13H,O. This might 
help to explain the discrepancies in reported re- 
fractive indices (as has been suggested by Bessey [4]), 
as well as the presumed multiplicity of forms of tetra- 
calcium aluminate hydrate. By reference to tables 
2 and 3 it may be noted that the carboaluminate 
was present in mixtures including almost the entire 
range of lime concentration. It may be significant 
that D’Ans and Eick [6] showed a similar range for 
the solubility of their §-tetracalcium aluminate 
hydrate, whereas the a-form was confined to the more 
basic region approximating that indicated by Wells, 
Clarke, and McMurdie [5] for the tetracalcium 
compound. 

Another possible source of confusion is the ap- 
parent instability of 2CaO-Al,O;8H,O on drying. 
As was brought out above, this compound by itself 
underwent a dehydration to 2CaQ-Al,O,;-6H,O, on 
drying over CaCl,, and the product gave an X-ray 
diffraction pattern that might conceivably have been 
mistaken for that of 4CaQ-Al,0,-13H,O. This re- 
action was reversible. A different phenomenon was 
observed with precipitates consisting of mixtures of 
di- and tetracalcium aluminate hydrates. After 
drying over CaCl,, such precipitates often gave 
patterns showing the lines of 4CaQ-AlI,O;-13H,O 
only. Presumably the dicalctum compound had 
undergone decomposition to the tetracalcium com- 
pound and alumina.’ This effect was not reversible. 
For precipitates of this type it was found necessary 
to obtain the X-ray pattern while the material was 
moist. Obviously this change on drying could lead 
to a false conclusion as to the solid phases prese 
in contact with the solution. 


3.2. Relation to the System at 21° C 


If the equilibrium diagram (fig. 4) is studied in 
comparison with the diagram published by Wells, 
Clarke, and MeMurdie [5] for the same system at 
21° C, certain marked differences are noted, as 
well as some similarities. Curve ABCD in figure 4 
has its counterpart in the 21° C diagram, and the 
same area of contusion between di- and _ tetra- 
calcium aluminates exists in the approximate range 
CD. However, the 21° C diagram does not show 
a stable solubility curve for 4CaQ-Al,O,;-13H,O com- 
parable to FC. It should be noted in this connec- 
tion that D’Ans and Eick [6] show a metastable 
solubility curve for the 8 form of this hydrate in 
the general area of FC. The from the 
1° diagram of a stability field for the isometric 
3CaQ-Al,0;-6H,0 is also a difference of major 
importance. The existence of a metastable solu- 
bility curve for monocalcium aluminate hydrate at 
1°, of course, was not unexpected. The other differ- 
ences between the diagrams for 1° and 21° are more 
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? While this paper was in process of publication, two other papers bearing in 
part on this phenomenon appeared in the literature. The authors, Roberts![18] 
and Schippa and Turriziani [19], report similar observations, but offer a different 
explanation According to their views, the X-ray lines here attributed to di- 
calcium aluminate hydrate may actually represent a different form of tetra- 
calcium aluminate hydrate (with 19 H30, according to Roberts). 


26 
~ 05} 
oe 
z 
© 04 
_ 
é 
5 93 
9 
$0.2 SCED* ALS, OM 
- Ca0:-A 2° HC 
~ " Calo 
4CoO-AbOx 13H Ca(OH 
¢ GIBBSITE —_. J 
= r Se 
02 a 4 5 6 8 9 0 3 ‘ 
ocO IN L T N ; 
Fiaure 4. Portion of the system CaQO-Al,03;-H,0 at 1° C, 


ABCD, metastable equilibrium curve for 2CaO-Al.0;3-8H20; EB, metastabie 
solubility curve for CaO-AlO3-10H20; GF, (assumed) stable equilibrign 
curve for gibbsite; FCD, stable equilibrium curve for 4CaO-Al2,03.13H,0, 


surprising, and it appears that further studies at 
one or more intermediate temperatures may be 
needed to define all the stability and metasta bility 
areas. 


3.3. Application to Cement Hydration 


The system CaQ-Al,O;-H,0 is not directly appli- 
cable to portland cement, which is predominantly a 
calcium silicate cement, with enough calcium sulfate 
added to combine with most of the alumina in the 
early stages. On the other hand, it has been shown 
[10] that the hydration of aluminous cements, at 
least to a first approximation, can be explained in 
terms of reactions in the system CaQ-Al,0,-H,0. 
Aluminous cements are known to react with water 
to form highly supersaturated calcium aluminate 
solutions slightly more basic than monocalcium alu- 
minate. The crystalline phase that precipitates first 
may be CaQ-Al,O,;-10H,O or 2CaQ-Al,O,-8H,O, or a 
mixture of these, depending on temperature, basicity, 
and probably other factors. From the data in table 2 
it may be concluded that, at 1° C, CaO-Al,O,-10H,0 
either would precipitate first or would form rapidly 
after an initial precipitation of the dicaleitum com- 
pound. Earlier work [10] indicates that, at room 
temperature, 2CaQ-Al,O;-8H,O forms initially and 
persists for long periods. However, it is clear from 
other data in the same report that the reactions of 
cement in paste form are not necessarily identical 
with those in the presence of a large excess of water. 
Thus CaQ-Al,O,-10H,O, rather than 2CaQ-Al,0; 
8SH,0, was shown by X-ray diffraction to be the 
predominant phase in pastes of aluminous cements 
after hydration at room temperature (24° C). From 
the data in the present report, coupled with that 
given in an earlier paper [13], it is concluded that 
monocalcium aluminate hydrate is also the chief 
crystalline product of hydration of aluminous ce 
ments at 1° C. 

In connection with the hydration of cements, the 
carboaluminate, 3CaQ-Al,O,-CaCO,-11H,O, should 
also be mentioned. From tables 2 and 3 it is apparent 
that this phase was encountered throughout. the 
range of lime concentrations studied, despite efforts 
to exclude CO,. Most other investigators likewise 
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have sought to eliminate the extra component. In 
practice, however, cement is always mixed and placed 
in the presence of atmospheric CO», hence the carbo- 
aluminate may form freely. Evidence of it has been 
found in the X-ray patterns of hardened pastes of 
aluminous cements. In one of the studies previously 
referred to [10] it was reported as “phase x’; the 
position of the diffraction lines now permits its 
identification as the carboaluminate. 

Some of the recent investigations of Farran [17] 
into the nature of the contact film between hardened 
cement and aggregate are also of interest in this 
connection. Farran found that aluminous cement 
appeared to react superficially with calcite aggregate. 
The X-ray pattern of the reaction product included 
a series of lines that agree well with the pattern 
obtained in the present study for the carboalumi- 
nate. Farran was of the opinion that the pattern 
indicated the presence of a solid solution between 
4Ca0-Al,O;-13H,0 and a carboaluminate. Such a 
composition previously had been suggested by Bessey 
4] for the mineral hydrocalumite. In the present 
work, however, and solely on the basis of the X-ray 
patterns, the carboaluminate and 4CaO-Al,0;-13H,O 
appear to exist side by side, w ith no evidence of solid 
solution between the two. 

These questions cannot be resolved without further 
data. There can be little doubt, however, that 
earbon dioxide may play a significant role during, 
and subsequent to, tne hardening of aluminous 
cements, and probably of other hydraulic cements as 


well. 
4. Summary 


A study of the system CaQ-Al,O,-H,O at 1° C has 
been made. The stable phases appear to be gibbsite 
(Al,O.-3H.O), 4CaQ-Al,O;-13H,O, and Ca(OH)». 
Metastable equilibrium curves were determined for 
2CaO0-Al.O;-SH.O and CaQO-Al,O,-10H.O. <A meta- 
stable invariant point for the latter two hydrates 
and solution exists at about 0.10 g of Al,O; and 0.38 
g of CaO per liter. The compounds 2CaQ-Al,O;- 
8SH,O and 4CaQ-Al,O,-13H,O crystallize together 
over a wide range of concentration, but the former is 
presumed to be metastable throughout the range. 
The isometric 3CaQ-Al,O;-6H,O appears to be meta- 
stable at 1° C, being transformed very slowly to one 
or more of the hexagonal hydrates. 

No evidence was found of the existence of a hex- 
agonal tricalcium aluminate hydrate. Only one form 
of the tetracalcium aluminate hydrate was observed. 
A similar phase, which was observed in minor 
amounts in most of the solid residues in contact with 
solutions, especially after prolonged storage, was 
identified as a carboaluminate, probably 3CaO-Al,O;- 
CaCO;-11H,O. The hydrate 2CaO-Al,O,-8H,0, on be- 
ing dried over CaCl., was found to lose 2 molecules of 
water reversibly. The lower hydrate has a character- 
istic X-ray pattern. Coprecipitated mixtures of 
2CaO0-Al,0;-8H,O and 4CaQ.-Al,O,;-13H,0, after dry- 
ing over CaCl,, showed lines only of 4CaQ-Al,O;- 
13H,0 on the X-ray patterns. The dicalcium com- 
pound presumably was decomposed to the tetra- 
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calcium and amorphous alumina. This reaction was 
not reversible. 

Owing to the extreme slowness of dissolution at 
: ° C, the solubility curve of gibbsite was not estab- 
ished. 


During the earlier stages of this study the author 
was ably assisted by G. J. Halley. The X-ray 
diffraction patterns were prepared by H. E. Swanson 
and coworkers of the Microstructure Section. Their 
assistance is hereby gratefully acknowledged. 
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Heat Content of Zirconium and of Five Compositions of 
Zirconium Hydride from O° to 900° C' 


Thomas B. Douglas and Andrew C. Victor 


Using a precise Bunsen ice calorimeter, a “silver core” furnace, and a “drop”? method, 
the heat content (enthalpy) of zirconium metal and five samples of zirconium hydride (24 
to 52 atomic percent hydrogen) was measured over the range 0° to 900° C. Thermal 
hysteresis of the hydrides was investigated in several cases, and corrections were applied 
for the impurities in the samples. At each temperature up to 550° C the heat contents 
of the hydrides relative to 0° C were found to vary linearly with composition, a fact in 
accord with the phase diagram of the system. ; ; 
sample reached its maximum value between 550° and 600° C, a behavior that can be inter- 
preted as due to the endothermic formation of the high-temperature beta phase. 


1. Introduction 


The fact that palladium metal absorbs consider- 
able amounts of hydrogen gas has been known for a 
very long time. Although many metals form similar 
alloylike solid hydrides, those of the family of chem- 
ical elements titanium, zirconium, hafnium, and tho- 
rium have become well known because of their ease 
of formation, wide ranges of composition, and relative 
stability. In the case of zirconium no gaseous hy- 
drides are known, but the metal expands and absorbs 
continuously variable amounts of hydrogen even at 
room temperature, forming brittle, metalliclike ma- 
terials whose compositions approach ZrH,. Only 
those with the largest amounts of hydrogen have 
decomposition pressures that exceed 1 atm below 
1,000° C. 

Many of the physical properties of the zirconium 
hydrides have been investigated intensively, espe- 
cially in recent vears. These properties include their 
erystal structures as determined by X-ray and neu- 
tron diffraction, their paramagnetic susceptibilities, 
and the diffusion rates and equilibrium pressures of 
hydrogen at various temperatures and compositions. 
This paper reports an extensive series of measure- 
ments of the heat content of zirconium and five 
compositions of its hydrides from 0° to 900° C. 


2. Samples and Their Compositions 


The samples of zirconium hydride whose heat con- 
tents were measured are listed in table 1 in the order 
of increasing hydrogen content. Their compositions, 
corrected for chemical impurities in the manner pres- 
ently to be described, are expressed in several forms 
that are commonly used. Zirconium metal is in- 
cluded with the hydrides, as it is a member of the 
series not only in a formal sense but also in terms of 
its properties. The zirconium sample was crystal- 
bar material, and was obtained from the U.S. Bureau 
of Mines, Albany, Oregon. 

In the case of pure zirconium hydrides the three 
ways of expressing the hydrogen content used in 
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Center, under USAF Delivery Order No. AF 33 (616) 56-21 The contract was 
initiated under Project No. 7360, ‘‘Materials Analysis and Evaluation Tech- 
niques,” Task No. 73603, ‘“‘Thermal Measurements,” 


13 


The effective heat capacity of each hydride 


TABLE 1. Samples and their compositions 














Composition 

Sampk Source of sample et 
| Weight % | Atom % Atomic 

H(w) | H(A) ratio H/Zr 

(x) 

1 U. 8. Bureau of Mines 0. 000 0.0 | 0. 000 
2 Wright Air Dev. Ctr . 356 24.5 . 324 
3 Sylvania . 611 35. 7 . 556 
4 do . 79 41.2 . 701 
do 1, 092 50.0 . 999 
in Wright Air Dev, Ctr 1.170 51.7 1. 071 








table 1 are related to one another, and to the hydro- 
gen as parts per million by weight (ppm), by the 
following numerical formulas. (The symbols are 
defined in table 1.) 


ppm = 10,000w= 110.5.A/(1—0.00989.A) = 


11,0502/(1 +-0.011052) (1) 
A=90.50w/(1+0.895w) = 1002/(1 +2) (2) 
r=0.905w/(1—0.01w)=.A/(100—A). (3) 


The methods used in preparing the hydride sam- 
ples were as follows: Samples 2 and 6 were prepared 
by the Power Plant Laboratory of the Wright Air 
Development Center. The zirconium used had been 
prepared by the Carborundum Metal Corporation. 
The hydriding temperature was 1,800°F (982°C), 
and a period of 2 hr was allowed in order to preclude 
temperature gradients along the bar. The furnace 
was shut down thereafter and the bar cooled in the 
furnace for at least 8 hr. The percentages of hydro- 
gen furnished with the samples were determined by 
the manufacturer by three different methods, de- 
scribed as vacuum fusion, weight gain, and volu- 
metric analysis. The three methods agreed very 
closely, and the estimated accuracy was stated to be 
of the order of 0.5 percent of the hydrogen present. 
Samples 2 and 6 were each supplied in two pieces, 
which were measured as a single sample to increase 
the precision, as in each case both pieces were stated 
to have the same composition and were found to have 
room-temperature densities that were identical 
within the uncertainty of measurement. 








Samples 3, 4, and 5 were supplied by Sylvania 
Electric Products, Inc., Atomic Energy Division. 
It was stated that at the time of their preparation 
only a commercial grade of zirconium was available. 
The metal was heated to 920°C, and exposure to 
hydrogen gas was then begun. As the system 
slowly cooled (reaching 800°C in about 45 min), 
the pressure of hydrogen was maintained equal to 
the equilibrium pressure of the solid phase down to 
0.2 atm. No further adjustments of pressure were 
made during the subsequent slow cooling overnight 
to room temperature, because it was believed that 
changes in the composition of the solid during this 
time were insignificant. It was stated that the 
greatest difference in temperature in any one sample 
did not exceed 4 deg C. ‘Tests for homogeneity were 
not made on these particular samples, but were on 
others, which showed no composition differences 
within the accuracy of analyzing for hydrogen. The 
compositions of samples 3, 4, and 5 given in table 1 
are based on analyses at Sylvania Electric Products, 
Inc., using a combustion method. Additional speci- 
mens that were approximate duplicates of these three 
compositions also were supplied. No heat-content 
measurements were made on them, but the parallel 
specimens were ground to a particle size of 60 to 100 
mesh and were analyzed at the National Bureau of 
Standards for carbon and hydrogen by a combustion 
method. The hydrogen contents thus found differed 
from the values reported by Sylvania by amounts 
that, expressed as percentages of the hydrogen 
present, varied from —1.5 to +3 percent. , 

The chemical impurities found by analysis in the 
six samples are given in detail in table 2. These 
values were obtained by qualitative spectrographic 
analysis, except where otherwise indicated, and pre- 


TABLE 2. Impurities in the zirconium hydride samples 
(Percentage by weight 
Samp 
Element 
1 2 3 { f 
ae 0.0145 0.01 20.600 »0. 44 a (). 48 0.01 
Fe 030 l a &S es 53 s 41 l 
Gc b. (2 b. 06 t ; 32 $2 b 06 
N. b . 005 Ol Ol 
oO b 006 017 O17 
Si 02 e 48 a 3 a { 
Ni 002 007 a (4 as 08 » 07 007 
Cr OOo ol 10 se OH 06 ol 
Al O02 Ol 12 a 12 17 ol 
| Mn ool Ol n Ol . Ol Ol Ol 
Ti 0035 5 06 ‘ Ol . 01 
Sa < OO1 OOS . OS . ol " Ol OOs 
Seoeaee OOl 005 ol ol 01 OOS 
Sn < OO0OS5 oo OO Ol Ol OOS 
Be 0001 0001 
Vv —s y 00 OOO] OO01 0001 005 
Ag a a ee OO1 OO! OO! 
| Ca anaes eer oe 4 Ol <. 01 Ol 
_ ane Ool cai ane Ol Ol Ol 
| Zn 005 
i 00005 -_ . 0001 0001 OOO! - 
B m 00005 Ooo1 Ooo! OoOd 
Mo OO Ooo! Ooo! Ooo! 
Co 0005 oool Ooo! Ooo1 
Total. % 0.09 0. 25 2.6 1.9 1.9 0. 25 
’ 


* By quantitative spectrographic analysis. » By chemical analysis 
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sumably no value was determined where none jg 
listed. Samples 3, 4, and 5 themselves were not 
analyzed for carbon, but the value given is probably 
approximately correct, as it is the mean of those 
found (0.29 to 0.35%) for the duplicates of these 
samples. ; 
The method of correcting the hydrogen COMposi- 
tions for the impurities shown by analysis is ¢gp. 
sistent with the method of correcting the thermal 
values described later, and was as follows: The totg] 
mass of those impurities whose small atomic sizes 
indicate clearly that they must be interstitial instead 
of substitutional in the zirconium lattice (carbop 
nitrogen, and oxygen) were subtracted from the total 
sample mass. The remaining impurities were they 
assumed to be replaced atom for atom by zirconium 
before recalculating a corrected sample mass and q 
corrected percentage by weight of hydrogen. 


3. Measurements of Heat Content 


The calorimetric apparatus and method used hay 
been described in critical detail in several reports 
and papers of the Bureau during the past 10 year 
[1].” 

In the heat-content measurements the sample was 
contained in a thin-walled cylinder of the alloy 
80Ni-20Cr sealed by a gold gasket. The hydride 
samples were not clad with any protective coating 
to prevent the escape of hydrogen, for in the worst 
case (at 900° C and 1l-atm pressure of hydrogen) th 
4 em® of gas space in the container will hold onh 
0.08 mg of hydrogen gas, and this is only about 0.04 
percent of the total hydrogen in the sample used 
Small losses in mass attributed to loss of hydroger 
did occur at the highest temperatures (800° C and 
above), as discussed later, but these aggregate losse 
equaled only a small fraction of the hydrogen in th 
solid sample. As the container appeared to remail 
sealed in most cases, these losses of hydrogen pre 
sumably occurred by diffusion through the containe 
walls. 

Briefly, the method of heat measurement was & 
follows: The specimen and container were suspended 
in the furnace for sufficient time to ensure ther 
reaching the known constant temperature of the 
furnace, which was always adjusted to a rounded 
value. They were then allowed to fall into a pre 
cision Bunsen ice calorimeter, which, by absorbing 
mercury equivalent in volume to the shrinkage 
caused by the melting of ice, measured the heat 
evolved by the sample and container in cooling te 
0°C. A similar measurement on the empty con 
tainer at the same furnace temperature provided 
correction for the heat the container itself contributed 
and for the small amount of heat lost during th 
fall into the calorimeter. The furnace temperatur 
was measured by a calibrated platinum-rhodium 
thermocouple, and below 650° C also by a calibrated 
platinum resistance thermometer to increase th 
accuracy and guard against any otherwise unknowl 


? Figures in brackets indicate the literature references at the end of this papt 


he 
tv 
te 
he 
se 
Ol 
al: 
at 
the 
te! 
tic 
zir 
th 
th 
be 


CO) 
pa 
L 


ha 
hig 
ba: 
fro 
rel; 
tur 
san 
anc 





none js 
rere not 
robably 
of those 
Of these 


‘OMposi- 
1s con. 
thermal 
‘he total 
Nic sizes 
| instead 
(carbon. 
the total 
ere then 
irconium 
ss and q 


nt 


sed have 
| reports 
10 years 


nple was 
he alloy 
hydride 
» coating 
he worst 
yen) the 
iold only 
yout 0.04 
ple used 
hv drogen 
° C and 
ate losses 
en in the 
o remail 
gen pre 
contained 


it Was as 
uspended 
ure their 
‘e of the 
rounded 
to a pre- 
absorbing 
shrinkage 
the heat 
‘ooling to 
ipty con- 
provided 
ntributed 
uring the 
n perature 
-rhodium 
-alibrated 
rease the 
unknow! 


i of this pape. 





shift in the calibration of either the thermometer or 
the thermocouple. Before being sealed, the sample 
container was filled with pure helium gas, which also 
flowed through calorimeter and furnace cores during 
the measurements. ‘This minimized the times re- 
quired to reach thermal equilibrium and the oxida- 
tion of the container surfaces. 

The heat leak of the calorimeter is normally small 
and quite constant, averaging only several tenths of 
a calorie per hour. In many of the measurements 
on the zirconium hydrides, particularly those with 
the highest hydrogen contents, the sample con- 
tinued to evolve heat in the calorimeter (up to a 
total of 1 cal extra) for as long as 1% hr beyond the 
90 min normally required for reaching thermal equi- 
librium. This slow evolution of heat was presum- 
ably due to a slow approach to final phase equilib- 
rium, and was taken into account when the precision 
of the measurements seemed to justify doing so. 
One of the authors has used the heat-content data 
reported in this paper to derive, by two independent 
methods, the shape of one of the boundaries of the 
zirconium-hydrogen phase diagram (the hydrogen- 
rich a boundary) from 0° to 550° C. The approxi- 
mate agreement of the results, using the two methods, 
is evidence that in this temperature range approxi- 
mate phase and thermal equilibrium were obtained 
in the present heat-content measurements [2]. 

The basic heat-content data obtained are given in 
detail in tables 3 to 8. Heat is expressed in calories 
(1 eal=4.1840 absolute joules). The second column 
indicates the relative chronological order of the indi- 
vidual runs with the sample, so that the effect of 
intervening exposure of the sample to other tempera- 
tures may be traced. 

Each value recorded in the sixth column for the 
heat found for the empty container is the mean for 
two or more measurements at the particular furnace 
temperature, with the following exception. The 
heat content of the empty container was measured 
several times during the series of measurements, but 
only at 100-deg-C intervals from 100° to 900° C (and 
also at 550° C), and the values for all other temper- 
atures were found from smooth plots of the devia- 
tions from simple monotonic empirical functions of 
temperature. In the present case, such interpola- 
tion led to errors less than those caused by the 
zirconium hydride samples themselves, although in 
the case of samples capable of being measured with 
the highest precision, this procedure would hardly 
be accurate enough. Different values for the empty 
container at the same temperature arose when in 
cases of emergency a new container or container 
part had to be substituted. 

The tabulation ‘corrections for losses of H,’’ that 
had occurred during exposure of the sample to the 
highest temperatures were made on an empirical 
basis, being estimated from the drift of heat values 
from run to run at the same temperature, from the 
relative duration of exposures to the high tempera- 
ture, and from the losses of mass undergone by the 
sample. At two different temperatures, such as 850° 
and 900° C, it would be expected that the rates of drift 
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of the heat value and the sample mass would be ap- 
proximately proportional to the equilibrium pres- 
sures of hydrogen. This relationship was found to 
hold roughly, and led to the belief that the correc- 
tions applied are of at least the right order of magni- 
tude, especially as the drifts were often decidedly 
outside the usual precision of measurement with this 
apparatus. (Presumably no problem of possibly 
sluggish reproportionation between two solid phases 
existed above 800° C.) It will be noted that this 
method of correction leads to exact agreement in the 
next-to-the-last column when only two measurements 
were made at the given temperature. 

Alternatively, it would have been possible, in 
principle, to apply corrections for the losses of hydro- 
gen calculated from a comprehensive correlation of 
the heat-content data with the phase diagram of the 
zirconium-hydrogen system. However, this was not 
practical for two reasons. In the first place, while it 
appeared that two of the samples each lost a total of 
8 mg of hydrogen (most of this at 900° C), some of 
the loss may have been obscured by small gains of 
oxygen simulataneously acquired by the sample, and 
the two corrections are widely different per milligram 
change in sample mass. In the second place, the 
theoretical correlation would often have led to cor- 
rections that were uncertain because they involved 
small differences between relatively large numbers. 
It is interesting to note, however, that such a correla- 
tion was made and that in nearly every case the 
corrections so predicted had not only the same sign 
but also the same order of magnitude as those found 
empirically. 

The hydride compositions were corrected for the 
chemical impurities by subtracting the mass of car- 
bon, nitrogen, and oxygen and then assuming the re- 
maining impurities to be replaced by the same number 
of zirconium atoms. The corresponding corrections 
to the heat values shown in tables 3 to 8 were made 
in two steps. First, the additively calculated heat 
contributions by the interstitial elements carbon, 
nitrogen, and oxygen were subtracted. Then it was 
assumed that the heat capacity of zirconium is the 
same per atom as that of the remaining impurities it 
was assumed to replace. However, in the case of 
carbon and silicon, appreciable deviations from equal- 
ity of atomic heat capacities at the same temperature 
were taken into account by using the actual heat 
capacities of the free elements. 

From the third columns of tables 3 to 8 it will be 
noted that in some runs the sample was first held for 
a while above the temperature of the measurement 
before being cooled for a final period of equilibration 
at the nominal temperature. (This elevated tem- 
perature was actually only 30 deg higher for zirco- 
nium at 870° C, and only 40 deg higher for the hy- 
drides at 500° C.) These runs, in which equilibrium 
was approached from a higher temperature, were 
made for comparison with those approached from a 
lower temperature in order to set limite to the error 
in heat content that could have resulted from lack 
of phase equilibrium, and separate mean values have 
been recorded in the last columns of the table to show 





TABLE 3. 





rime 1 
Furnace Order in if 
temper- which Sample 
sure, f run was mass 
mace With 
(+50 t sample 
fi mi min } cal 
if { 31 20. 0753 258. ¢ 

li v0 31 20. OF 54 JSS. f 
100 } 24 0 3y 20. 0757 28S. 6 

| 25 0 44 20. 07 28k. 

1s ) 38 “). O757 505.3 
19 0 33 20. 0757 504. 3 
A) 0 37 20. 0757 503. 2 
21 0 1) 2). 0757 593. 4 

} 22 0 tt 20. 0757 502. 2 

aw \ 23 ) 3t 20. 0757 593. 8 
) ) 41 20). O757 504.9 

27 0 32 20. O757 592. 7 

| 45 0 42 20. 0757 592.9 

16 Tt 53 4%). O757 593. 6 

7 0 43 20. O757 916.5 

s 0 35 20). O7 57 915.1 
300 9 0 47 0). O757 915.4 

| 10 0 18 ”). OTST 915.0 

ll 0 $1 20. O757 915.0 

12 0 31 20. O757 1240. 6 
13 0 31 0. 0757 1240. 4 
400 47 0 17 20. O757 1239. 2 
| 4s 0 37 0), O757 1240.8 
{ 49 0 51 20. 0757 1241.1 
} v0 30 at). 7 1574.4 
4 0 30 ( 757 1578. 0 
“ 0 $1 Pa 07 7 1576 9 
6 { {3 A), ( 157 
14 0 $1 7 1575 
15 0 $4 20. 07 1575. 2 
600 l { 3 a] r 1923. 6 
l 2 0 ) 20. O757 1923. 1 

28 { t f 7 2200). 
H 24 ( 7 2288. 7 

| 30) { 42 2 7 2288. 7 

| 31 ( 0). OTST 206. 4 
SOO s2 ; 2) Mod 2HH6. 4 

} 37 ( ”). OT57 onM4. & 

3s 0 ti 2%). O757 HH5. 4 

33 0 64 0, O757 2863.8 
1 0 h2 0. O757 284. 3 
2 ; 0 2 0. { 7 HT 6 
S00 1) 0 3 0). OTST ORO. 7 
Q 0 7s 15. O77¢ 2528. 3 
Ho 0 iH 15. O776 4, ee) 
- j 44 5s 2 20. 0757 3123.9 
St = - « 

\ 41 O4 ) at 3123. 2 
a7 i 0 Th 15. O77¢ 27H 
SSO) 55 0 57 15. O776 2781 

4) 0 63 ] O77 2801. 3 
885 v2 v 72 15, O776 2803. § 

| 5S 0 0 15. 077 9709. 5 
Su“y) ‘ 0 61 ] O77¢ ORTR 
RY 57 0 57 15. 0776 2835. 2 

$5 0 fZ a) O757 325. 1 
a 0 i) MM). O57 3263. 3 
4] ( Ph ") OTAT 275K § 
Yoo - « 
_— 42 0 61 20. O757 3258.8 
51 0 #2 15. O77 9854. 2 
53 0 | 1 O776 9853. 5 
s Omitted from the mean becaure of inferior precision 


these differences. It will be noted, as one would ex- 
pect, that in every case the higher heat was found 
when the approach was from a higher temperature. 
Also, longer times at the nominal temperature usually 
narrowed the discrepancy between the two directions 
of approach, though not rapidly enough to be a 
practical means of coming as close to equilibrium as 
the precision of the apparatus enables one to measure. 


Individual measurements of heat content 


Heat observed 
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zirconium metal (sample 1) 


Corrections for com- Corrected H HHo°c, 
position net for sample 


Container | Loss of Hy | Impurities | Single run Mean 
only set 
cal cal cal @ cal ¢ cal g 
150. 0 0. 00 0. 00 6.91 
1M). 0 .00 oO 6. 91 | 
- 6. ¥l 
150.0 Oo 00 6. 90 
1s. 0 Ow OO 6. ¥l 
310. ¢ 00 01 14.17 
310. ¢ UU Ol 14.12 | 
310. ¢ 00 Ol 14. 07 | 
310.6 00 Ol 14. OS 
310. 6 OO Ol 14. 02 
310. ¢ 00 Ol 14.10 f 14. 08 
310. f om Ol 14. 15 
310.6 oo Ol 14. 04 
311. ¢ OO ol 14. 00 
311. ¢ OO Ol 14.04 
476.7 OO 4 8 (21.87 
476.7 OO 4 21, St 
476. 7 oo — 04 21. 81 1 SO 
476. 7 00 4 21.79 
476. 7 00 04 21. 79 
648. 2 oO 03 29, 48 
(48, 2 Ow O38 2Y. 4 | 
651.0 oo 03 29 36 29. 44 
49.2 Oo —. 03 20. 44 
bAY 29. 45 j 
825.0 Ooo 04 37. 29 
825.0 Ooo 4 37. 47 
R25. 0 oo 04 t7. 41 on 
825.0 oo 04 37. 35 D4. Dé 
825. 0 on 04 $7. 34 
R25. ( ow — (4 37. 33 } 

10090 oo (Vv 45.47 } . 
1009. 5 Ot 04 45. 45 a5, At 
1204. 0 Oo OS 4. 02 
1214. 0 Oo OS nd. OS 3. O7 
1204. { ”) OS 53. 95 | 
1408.0 x OS 62. 61 
1408. 0 “ OS #2. 61 | F - 
1408. 3 Oo OS 62. 5l | o-_ . 
1408. 3 00 Os 62. 54 | 
1510.8 . 00 10 67. 3 
1510.8 00 10 67. 32 
1511.1 rT 10 67.17 a 
1511.1 00 10 67.14 _ 
1513 Oo 10 67, 2 
1513. 5 00 10 67. 23 
1552. 7 00 09 78.17 ee 

&. If 
1552. 7 00 09 it" a 
Lat f x oo a 70. lf 
1576. 1 n 09 79. St 79. St 
1586. 5 mM 9 mt) 4S } 
1586, 00 oy SD). 62 } SO. 4U 
1586 oo ou Mi). 4 
1597.0 oo ov RO. Ol SO. U1 
1007 OO 09 81.33 81.33 
1614.3 00 09 $2.14 
1614.3 on OW RD? OF 
1615 } im) oa 1 7 81.9] 
161 4 Ow — Ss] i 
1618.0 oo Ho RT oO 
1618.0 ) 09 RI. &5 } 


When a hydride sample had been measured“at high 
temperatures, such as 900° C, before being remeas- 
ured at some lower temperature in the manner 
described above, the remeasurements were begun in 
the usual wav by approaching the furnace tempera- 
ture from a lower temperature, and the corrections 
applied for previous loss of hydrogen were calculated 
so as to give perfect agreement with the results ob- 
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TABLE 4. Individual measurements of heat conten 
Time interval 
Furnace Order in at 
temper- which Sample bs 
iture, t run was mass 
made With 
(+50 t sample 
( min min g cal 
{ 13 0 67 19. 8280 206. 7 
100... 1 14 0 43 19. 8280 296. 4 
. f 15 0 42 19. 8280 619.7 
0. - | 16 0 R0 19. 8280 618.9 
. 6 0 63 19. S280 963.9 
300. - . 7 0 53 19. 8280 963. 5 
f 4 0 64 19. S280 1337. 6 
40K 5 0 50 19. 8280 1338. 7 
l 0 61 19. 8280 1733. 1 
500 2 0 48 19. S280 1733. 2 
| 3 40 71 19. 8280 1733. 4 
. 23 0 69 19. 1375 1881.0 
jad 24 0 52 19. 1375 1882. 2 
{ 8 0 64 19. 8280 2062. 1 
550_.. 4 0 4 19. S280 2032. 2 
| 22 0 78 19. 1375 2028. 8 
10 0 52 19, S280 2561.8 
11 0 72 19. 8280 2560. 9 
ww 12 40) 74 19. 8280 2565.3 
17 4) v6 19. 8280 2566.8 
f 18 40 ai) 19. S280 2741.6 
ww 1 19 0 87 19. 8280 2741.1 
f 20 0 70 19. 8280 3083. 1 
SOU , OF 0 57 19, 8280 3079.3 
y 25 0 61 19. 1375 3675. 6 
JOU 1 26 0 76 19. 1375 3668. 9 


*s Given no weight. 


tained before the exposure to higher temperatures. 
This method of applying corrections for previous loss 
of hydrogen did not, of course, affect the principal 
purpose of making a new series of measurements, 
which was to ascertain the difference in heat resulting 
from heating and from cooling the sample to the de- 
sired furnace temperature. 


4. Mean Heat Capacities of Zirconium and 
the Zirconium Hydride Samples 


The heat content of zirconium has been measured 
by Mixter and Dana (0° to 100° C) [3], Jaeger and 
Veenstra (21° to 801° C) [4], Coughlin and King 
25° to 1,098° C) [5], Skinner (25° to 1,000° C) [6], 
Redmond and Lones (150° to 1,050° C) [7], and 
Scott (60° to 960° C) [8]. Reviewing the first two 
investigations, Kelley [9] gave equations for the 
heat content and heat capacity over the range 25° 
to 627° C. Seott’s zirconium had a purity com- 
parable to that indicated in table 2, but presumably 
most of the earlier measurements were on samples 
containing considerable hafnium for which no cor- 
rections were applied. According to Dulong and 
Petit’s law, in the case of zirconium each 1 percent 
of hafnium would lower the observed heat content 
and heat capacity per unit mass by about ' percent. 

Pure zirconium metal undergoes a transformation 
of crystalline structure at a temperature that, partly 
due to its sensitivity to small amounts of certain 
impurities, does not seem to be known very pre- 
cisely. Coughlin and King [5] observed a rapid 
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zirconium hydride 0.356% 


Corrections for com- 
position 


Corrected Hi— Ho’ c, 
net for sample 


Mean 


Container | Loss of Hy | Impurities Single run 
only (set) 
cal cal g-! cal g=! cal g-' cal g™! 
1%). 3 0. 00 0. 00 7.39 | 7.38 
150.3 00 .00 7.37 j si 
311.6 00 —. 02 15. 52 j - 
311.6 “00 —"02 3.43 jf 1550 
477.3 .00 —.04 24. 50 } ‘ 
477.3 . 00 —, 04 24. 48 f 24. 49 
650. 6 .00 —.05 34. 59 ] os 
650. 6 00 —.05 34.65 |f 34. 62 
812.6 .00 —.07 46. 36 \ _ of 
812.6 “00 = 6.35 |\jf %% 
$12.6 .00 —.07 46. 37 46. 37 
899.8 +.15 —. 08 51.34 \ ‘ 
899.8 —.15 —.08 51.41 | 51-38 
918.0 —.15 —. 08 57. 62 } 
918.0 —.15 —. O08 ®(56. 11) 57. 86 
918.0 +15 —. 08 58. 11 | 
1010. 5 +-.15 —. 08 78. 16 \ } 
1010.5 +15 ~_ 08 m1. |f 14 | 
1010. 5 +-.15 —. 08 78. 33 ] =o oF 
1010. 5 +15 —. 08 78. 41 { 78. 37 
970.1 oO —.11 89. 24 89. 24 
970.1 00 —. 11 89. 21 89. 21 
1132.9 +. 12 —.12 98. 36 | ‘ 
132.9 +. 22 12 98.26 |f %-3! 
1615.4 +. 35 —.13 107. 87 \ ” = 
1615.4 +. 48 -. 107.65 |f 107. %6 


Apparently the time in the furnace was inadequate, 


rise in heat content over a range of some 80 deg C. 
Since it is known that oxygen or nitrogen, when 
dissolved in the metal, has such an effect, they 
assumed their sample to contain one or both of 
these impurities. The heat-content measurements 
of Redmond and Lones [7] were of low precision, 
and they reported no evidence for a heat-absorbing 
transformation. Though two investigators [6, 8] 
reported 870° C as the transformation temperature, 
two compilations [10, 11] have given 862° and 863° 
C, respectively, and the last temperature has been 
somewhat arbitrarily assumed in this paper. (The 
assumption of another of these temperatures would 
have had very little effect on the value of the heat 
of transformation derived from the data.) 

The present authors were limited by their ap- 
paratus to temperatures not higher than 40 deg C 
above the transformation temperature. As _ the 
preliminary measurements at 870° and 900° C had 
indicated a value for the mean heat capacity between 
these two temperatures of 11 cal/g atom-deg C, a 
value that seems too high for this element, addi- 
tional measurements of heat content were made at 
5-deg intervals over this temperature range. The 
mean relative heat content at each of these temper- 
atures (last column of table 3) is plotted in figure 
1. Considering the very small temperature in- 
tervals and the large magnitudes of all these meas- 
ured heats, the precision is good; but (above the 
transformation temperature) there is a decided 
negative curvature at the lower temperatures, which 
may be due to small amounts of oxygen and nitrogen. 








TaBLe 5. Individual measurements of heat content—zirconium hydride, 0.611% H by weight (sample 3) 


Time inter, He 
Furnace | Order in at 
temper- which _ Sam pie 
ature, ¢ run was mass 
made With 
1+ sample 
Cc 7 cal 
1 ( j 19. 6516 291. 6 
9 ( { 19. 651¢ 308. 5 
100 ; ) 42 19, 6516 103. 7 
{ ( 52 19. 6516 105. 0 
. ( 4 1 O51 (40. 2 
AM 4 37 > ¢ ‘ 640.1 
” 0 i7 19. 6516 1011.4 
m 8 0 7 19. 6516 1011. 2 
) 0 1? 19. 6515 1402. 7 
| 10 ( s 19. 6516 1403. 0 
1 ( I 19, 651¢ 1402. 5 
19 0 10) 19. 6516 1404. 4 
400 20 0 1 19. 6516 1404 
on 0 19, 651¢ 1402. 7 
2t { i 19, 6516 1402. 6 
27 0 ‘ 19. 6516 1402.0 
Ss 0 41 19, Hit 1402. 9 
11 0 54 19. 6516 1831.5 
12 0 1“) 19. 6516 1S). ¢ 
10 0 2 18, SSH 1789. 6 
Wy 4] 0 71 18. RRM 1791.7 
4? 41 ’ 1S, S850 ‘'7YlL. 1 
{ tie | 18. S850 i, 700. 
. 
lf 0 ‘ 19. 6516 2063. 6 
( 17 0 19. 6516 262. 2 
| LS 0 4( 19. 6516 2061.8 
29 { i 19. 6516 2386, 1 
75 ai) 0 7 19. 651 2404.8 
| ; l 19. 651 2318. 4 
14 4] 19. 6516 2H. 2 
l i ~f 10. O51 20. 7 
‘ 1) LS. S850 2578. 4 
600... - ; . acm 
r 0 “I S. SSO) is 
Is 41 4 18, S850 282. 1 
iy 0 7s 18, S850 240. 1 
1 0 tp 19. 6516 9953. ¢ 
{2 0 y. 19. 6516 2055. 0 
‘ 44 0 i) 18, SS5O 2870. 4 
ow 45 0 61 18. R850 INTO. 
6 10 ti 18. S850 2871. ¢ 
17 10 us 1S. SS5O INT 1 
{ 91 0 54 19. 6516 174 
Ww ‘ 29 0 10) 19. 6516 3175. 1 
{ 23 0 10) 19. 6516 3624.9 
SUU | 24 0 1) 19, 6516 3623. 1 
{ 34 0 10 19. 6516 41059. 9 
oad 0 ’ 1Y. 6516 1053.8 
*® Omitted from the mean because of inferior pre 


The straight line was drawn on the assumption 
that the scatter of points at and above 880° C 
random. The mean heat capacity of 8-zirconium 
corresponding to the slope of this line is subject to 
considerable uncertainty because of the smallness 
of the total temperature interval. Extrapolation 
to the transformation temperature leads to a value 
of 9.8 cal g—' for the heat of transformation of the 
metal. The following values (in the same units) 
have previously been given: 4.7 [12], 7.7 [10], 10.1 
[5], 10.9 [8], and 11.4 [6]. The large spread of these 
values is undoubtedly largely due to the usual 
necessity of interpolating the heat content of this 
element through a considerable temperature range 
because the transformation indicated by the thermal 
data was not sharp. 

In table 9 are given the values of mean heat 
capacity of zirconium observed in the present in- 
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Corrected Ti— Fhec, 
net for sample 


Corrections for com- 
position 





Container Loss of H; Impurities | Single run Mean 
only set 
cal cal g cal ¢ cal g cal g 
150.4 0. 00 0.16 a(7. 02 
150. 4 oo 16 7. 63 - ¢ 
150. 4 oo 16 7. 64 | To 
150. 4 00 16 7. 60 j 
10.3 00 ) 16. 42 } . 
10.3 00 36 16. 42 16, 42 
475.3 00 58 6. 70 ~ "1 
4175.3 00 —_ 58 26, 69 - 
648. 0 00 S 7. 55 
H48. 0 Oo S5 7.57 
648 OO &5 7.55 
648. 0 00 s 64 
H48, 0 00 SJ 7. 64 ri 
648 00 & 7.56 
O48. | 00 8 7. 55 
HAS, ( 00 x 2 
648. 0 Oo sf 7. 57 
SOF " 1.12 0.0 } 

R96. 3 oo 1. 12 10, ON | an 
897.0 L O9 ui. 19.04 : 
827.0 09 1.12 0. OF 
827.0 i 09 1.12 0.0 e . 
827.0 (0 1.12 50. 00 0. U 
116.8 00 1. 20 7.07 | 
O16. 8 ”) 1. 29 7. Of 7. 02 
916.8 00 —1.29 6. US | 
962.9 ” —!] 9 70. 90 ] 
u62.4 " l 2 71. R45 , O1 
O62. 9 ) 1. 52 67. 4 | 
1009.8 OO 1.72 S242 } 
10090, 8 Oo 1.72 82.20 | , 
8 f 

1010. 5 £1. 05 1.72 &2 
1010 +1.0 1.72 82. 37 
LOO, 5 +1.0 1.72 85.7 ~ os 

— 85. ¢ 
1010. 5 1.0 —1.72 85. 62 . 
1106.8 0. 00 1. 90 72, O8 | 
1106, 8 Oo 1. 90 72.15 2. 12 
1115.7 +1. 11 —1. 90 212 || — 
1115.7 $1.11 1. 90 92.14 
1115.7 +1. 11 1. 90 2.19 |) 24 
1115.7 1. 11 —1.90 12. 28 } » 
1205. ¢ 0.00 —? 05 QS 13 } 98. 1 
1205. 6 OO 2.05 YS. 17 { 
1409. 1 +-. OF 2. 34 110. 48 { 
1409. 1 +. 16 2. 34 110. 48 a. 48 
1616. 1 +. 38 2. 61 122.13 |) 9 13 
1616. 1 oY 2. 61 122. 1 sss : 


vestigation. These~ values were calculated from 
the differences between successive values in the last 
column of table 3 up to 850° C, and from the line in 
figure 1 above 863° C. Also given for comparison 
are the percentage deviations of the smoothed values 
of other investigators from the values of the second 
column. 

Similarly, the mean observed heat capacities of 
the five samples of zirconium hydride whose heat 
contents were measured were calculated for the 
temperature intervals of measurement from. the 
successive differences in the last columns of tables 
4 to 8. These values are given in table 10. (The 
compositions may be identified from table 1.) In 
the intervals where the apparant heat capacity was 
found to be unusually large, particularly between 
550° and 600° C, extensive phase transformation 
was occurring, and as a consequence the effective 
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TaBLeE 6. Jndividual measurements of heat content 


Time interval 


Furnace Order in at— 
temper- which - Sample 
ature, ¢t run was mass 
made With 
t+50° t sample 
Cc min min g cal 
f 1 0 56 19, 3007 303, 2 
100 | 2 0 39 19. 3007 303. 2 
: 3 0 55 19. 3007 (42. 6 
20) 4 0 40 19. 3007 642.8 
0 56 19, 3007 1019.8 
WU. » 6 0 41 19. 3007 1019. 0 
, 0 58 19. 3007 1414.1 
4100 s 0 10 19. 3007 1413.6 
) 0 wa) 19. 3007 1847.4 
10 0 41 19. 3007 1847.7 
( 0 60 19. 3007 1847.6 
mn 31 0 57 19. 3007 1845. 8 
32 0 60 19. 3007 1845.8 
$3 0 60 19. 3007 1846. 7 
{ 13 0 6 19. 3007 2085.9 
0 ‘aie 14 0 187 19. 3007 2091. 1 
| 15 0 52 19. 3007 2084. 2 
{ ") 0 6 19. 3007 2217.4 
, 21 0 126 19. 3007 2199. 5 
| 2 0 62 19, 3007 2189. 6 
18 0 43 19, 3007 2508. 5 
19 0 SY 19. 3007 2505. 1 
ll 0 56 19. 3007 9722.3 
l 0 9 19. 3007 2723.9 
35 0 42 16. 7834 2507.0 
om 36 0 73 16. 7834 2510. 6 
10 61 16. 7834 2535, 8 
ss 40 95 16. 7834 2532. 4 
{ ) 0 51 19. 3007 3037. 2 
17 0 40) 19. 3007 3035. 2 
3] 0 43 16. 7834 2779. 5 
— 10) 0 6n 16. 7834 2779. 2 
41 11 66 16. 7834 2778.9 
412 41 90 16. 7834 2780. 2 
- 23 0 v1) 19. 3007 3252.5 
- 24 0 41 19. 3007 3251. 6 
| 0 45 19. 3007 3691. 4 
aM) ni 0 42 19, 3007 3688. 2 
27 0 15 19. 3007 3689. 1 
OR 0 5 19, 3007 4127.6 
YOO. 24 0 4) 19, 3007 4124.1 
| 34 0 50 19. 3007 4118.6 


Hieat observed 


zirconium hydride, 0.769% H by weight (sample 4) 


— | 
| 
| 


Corrected Hi— Tec, 
net for sample 


Corrections for com- 
position 


Container | Loss of Hz | Impurities | Single run Mean 
only (set) 
cal cal g=' cal g-' cal g~! cal g 
150. 4 0. 00 0.12 7.79 | 7 20 
150. 4 . 00 —.12 7. 80 j ms 
310.3 . 00 —. 29 16. 92 i] » on 
310.3 00 —. 29 16.94 |f 16. 93 
475.3 00 —. 48 27.73 \ 97.71 
475.3 . 00 —. 48 27. 69 j = 
648. 0 00 —. 71 38. 9S \ ae OF 
648. 0 00 —.71 | 38.96 If 38. 97 
826. 3 00 —. 92 51. 98 
826. 3 00 —.92 | 52. 00 
826.3 +. 05 —. 92 | 52. 05 as ¢ 
826.3 05 —, 92 51. 95 51. 99 
826.3 0 92 51. 95 
826. 3 0 —. 92 52. 00 
917.4 00 —1,06 9. 48 } 
917.4 OO 1. 06 9, 75 9. 54 
917.4 00 —1. 06 9. SY | 
935.9 .00 1.11 65. 29 } 
035. 9 00 111 64. 36 64.5 
935. 9 00 a, 63.85 || 
963. 7 .00 1,28 78. 76 } 78. 6 
963. 7 00 —1. 28 78. 538 j 
1009.8 00 1, 40 87.33 
1009.8 00 —1. 40 87. 41 = op 
1021.7 +17 1. 40 87.27 87. 34 
1021.7 .17 -1. 40 87.48 
1021.7 17 —1. 40 88. 99 | 
1021. 7 17 1. 40 any => 
1107.0 00 —1. 55 98. 46 | 
1107.0 . 00 1. 55 98. 35 
1115.7 83 —1. 55 93.42 if %%4! 
1115.7 83 —1. 55 98. 40 
1115.7 83 —1. 55 98. 38 \ an 
1115.7 +. 83 —1. 55 98.45 |f #2 
1205. 6 00 —1. 67 104. 38 ) es 
1205. 6 00 1.67 | ras jf 08 
1409. 1 1 07 1. 90 116, 42 } 
1409. 1 +. 14 —1. 90 116. 32 116. 39 
1409. 1 21 —1. 90 116.44 | 
1616. 1 +-. 48 —2.12 128. 48 | 
1616. 1 +. 68 2.12 128. 54 128. 48 
1616. 1 +. 93 2.12 128. 47 [ 





heat capacity was not only changing rapidly with 
temperature but the mean values recorded are 
subject to large uncertainties because of the impos- 
sibility of reaching equilibrium in the measurements 
within a reasonable time. On the other hand, in 
those temperature ranges where mean heat capacity 
appears to vary regularly with temperature (par- 
ticularly from 0° to 540° or 550° C), it is justifiable 
to assume the recorded values to be the instantaneous 
heat capacities in the middle of their respective 
temperature ranges and to interpolate for inter- 
mediate temperatures in the usual way. 

The regions of regular variation can be more 
quickly located by an inspection of figure 2, since 
for each composition, each two adjacent points are 
connected by a straight line whose slope is pro- 
portional to the corresponding value in table 10. 
It is also justifiable to interpolate for values for 
intermediate compositions at any temperature up 


to 540° C, as discussed in section 5. 


5. Discussion of the Heat Capacities of the 
Hydride Samples 
The observed discontinuities in the heat-capacity- 


temperature functions of the five samples of zirco- 
nium hydride can be interpreted qualitatively in 


terms of the phase diagram of the system. One 
version of this diagram is shown in figure 3. The 


three solid solutions indicated are designated a, 8, 
and y. <A similar diagram was given by Vaughan 
and Bridge as a result of their high-temperature 
X-ray diffraction studies [13]; but the results of dif- 
ferent investigators, using various methods, do not 
agree well as to the shapes of some of the phase-field 
boundaries, two of which seem to be particularly 
uncertain and are therefore shown as dashed lines. 
The dotted vertical lines and their points indicate 
the compositions and temperatures for which meas- 
ured values of relative heat content are reported in 
this paper. 
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TABLE 7. Individual measurements of heat content zirconium hydride, 1.029% H hy weight sample 7 
in erval leat observed Corrections for com Corrected Hi—Hoec, 
Furnace Order in at position net for sampk 
temper- whicl Samr 
ture, run was 
made With Container Loss of Hy Impuritic Single ru Mean 
t+ i sample only set 
Cc ’ cal cal j cal q cal q-! 
l 0 l 18. 4319 106, 5 150.0 0.00 0.12 Q 15 } . 
100 2 0 44 18, 9319 306, 5 150.0 00 12 8 15 8.15 
3 0 0 18. 9319 655. 7 210.6 Oo 0) 17 93 
4 0 +] IS. US19 656. 3 110.6 Oo +x 17. O68 
. 5 ( $1 18. 0319 657.5 310. ¢ oo bv) 18, 02 | 
otk 14 0 42 18. 9319 679.4 310.6 00 u) 19. 18 [ 18. 00 
15 0 5 18. 9319 658. 0 310.6 Oo 30 18. OF | 
lf ) 3 18. 9819 H58. 0 310. € oOo - 1) 18. OF | 
f 4) a) * 18. O319 1050. 1 476.7 ow 44 ”) St) 
7 0 +] 18. U319 1048. 7 476.7 mw 40 0 72 
= 8 0 15 18. 9319 1048. 4 476.7 Oo 1 29. 71 7 
ww 17 0 4 18, 9319 1046. 3 476.4 OO 10 20 Fi 9, OY 
18 0 ‘4 18. 0319 1045.8 476.4 Oo 14 24 r 
19 0 i. 18. O319 1049. 0 47 4 Oo 14 74 
ll 0 i 18. 93190 1462. 4 (48, 2 oo 7 (2 97 
12 0 4 18. 9319 1461.3 H4a. 2 OO 7 42. 22 
13 0 3 18. 9319 1461. 2 (48. 2 00 7 ‘2. 21 
100 33 0 4 18. 9319 144. 1 48. 2 On 7 $2. 37 42. 30 
34 0 (i 18. 9319 1464.9 (48,2 Oo 73 $2 4] | 
41 0 1s 18, 319 1456. 0 (48, 2 + 3S 73 42. 32 
42 0 61 18. 9319 1455 HAS. 2 +. 38 73 12.30) 
f Q ( 18. 9319 1910.4 St 0 mw ) h. 36 
10 0 { 18. 0319 1911.7 R25 0 ™ 7 6. 43 
' 43 f | 7 1796.3 812. ¢ L 30 97 6.39 P 
An 44 0 7 17. 2397 1706. 3 812.6 + 30 O7 ‘ 
45 4( 2 2307 1798. 7 812. ¢ ‘ 7 ; 
tt i) “ rf SOO. 8 SI. 4) 97 A ; 
v1) ( & ou ) 2164.4 917. 2 iM 1.11 (4 
, 21 0 7 8. 9319 2160. 3 917.2 0 1.11 64. 5 
- | 2 62 18. 9319 2158. 4 917.2 00 Li 64. 4 A. oi 
yA ( 18. USL9 J1LAS. 4 G17 Oo 1.11 4.4 
- f " 0 iO 18. 9319 2384. 2 ra { 0 1.2 S 2 . 
i 7 ‘ 64 8. 9319 2376. 8 3. 4 " 1.3 ‘ lin 
24 { 18 ) 2461.3 1009 " l an a4 
25 ( &. ) 2AH2. ¢ 1009 iM l s 7 7 
- RO. 07 
47 0 l 239 2408. 0 ou? 9 10 1.33 wt) ee 
OUU . is ( 72 ] 234 2408. 3 O92. 9 0 13 wu) ¢ 
4 4) 9 l {07 2416.0 o027 9 ) 1 a] 9 
7 : S1. 1 
”) 40) iv l 3Y7 2415.5 992.9 l ; SI ) 
an yan! 0 Hm 18. O319 ZRTs. 4 1107.9 ”) 1.49 > | op { 
o. 29 0 ’ 18. 9319 IR77.8 1107.9 0 1 49 > i «. Ve 
30 0 4 18. 9319 4179 1204. 0 . 1. 67 iT OS 
$1 ( | & USI19 4182. 4 124. 0 oo Lf 2 8 
32 ) 75 18. 319 sI81.9 1204. 0 TT 1 A7 > QI 102. 77 
700 1 ( ,. 2397 2001. ¢ 1182.9 ; 1. 67 ( 
2 0 PM) 208 1182.0 4 1.¢ ss 
; 435 f 20 S011. 1 1182.9 $ 1. 67 lf ; 
4 43 | 17. 220K 20905. 7 1182.9 4 1.¢ ( ( 10 
. ; 0 rr ) 478.7 1305.3 ”) 1 2 112. a8 13 
wv) 36 0 ( 18. 9319 4479.8 1305. 3 00 1. 82 113.04 113 
37 0 Pa) 8 O810 ah 4 1408. 0 Oo 1.97 123, - 
ow ? . ts 0 S 18. & ) 3740. 3 1408. 0 Oo 1. 97 123. 83 - 
40 ) ‘ a o ) s86, 2 1510.8 -~ 4 2.14 134. 30 
- 0) { 4 15. & } 1086, 4 1510.8 - 4 > 14 4 34 134 
RO ~ 17 206% 2R26 & 1510.8 4 2 14 134. 39 34. 32 
n ( l 20 3830. 5 1510.8 1s 214 134. 32 
, f 7 0 f 2H 144.3 1615.4 .. 26 t 141 
on | 5 0 ~ 17. 2260 4028. 4 1615. 4 8 2.2 j i41 
* Omitted from mean because o | ! Given no weight ADI y thet ‘ \ lequ 
ry - ‘ . ~~ ‘ rn , 
The mean observed heat contents up to 500° C of temperature up to 550° ¢ The sample ZrHo.x 


zirconium and the five hydride samples, taken from 
the last columns of tables 3 to 8 and converted to 
cal mole of ZrH,, are plotted against the over-all 
composition z in figure 4. The linearity of the-values 
for the five hydride samples is in ceneral rood and 
accords with the prediction of the phase diagram 
(fig. 3) that these five samples are mixtures in dif- 


ferent proportions of the same two phases at each 
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obeys this relation only up to 540° C In fact, the 
proximity of the observed points to the fairly steep 
lines at the higher temperatures may be considered 
to afford an approximate check on the general inter- 
consistency of the hvdride compositions found by 
analysis. The coefficients of the equations of the 
straight lines drawn through these points in figure 4 
were determined by the method of least squares and 
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TABLE 8 Individual measurements of heat content zirconium hydride, 1.170% H by weight (sample 6) 
ime interval Ileat observed Corrections for com- Corrected F:— Hooc, 
Furt O at position net for sample 
1M Samplk 
ur 4s 
( With Container | Loss of H Impurities | Single run Mean 
imple only (set 
( cal cal cal g cal g cal g=' cal g 
s Is il 276 122.9 0. Of 0.01 8 21 } g 93 
- , 276, 8 199 9 ) Ol 25 j 
( 18 l j 251.9 oo 03 18. 29 18. 29 
s 044.8 SR5_8 oo 04 30. 10 | 99 90 
x4 IS. 6311 M3 SS5. 8 OO 04 20. SS | = 
l 1321.0 24.2 ”) ( 42.72 42.72 
8 #31 1735. 2 HAG. G 0) Os 7.27 27 
{ 18. 6311 2 Hoe. ¢ Oo OS 57. 49 7.49 
® yw] 74 y ON 66 3 | » 
8 177.4 $ ) s 66. 31 { 66. 42 
s 2451 815.2 ’ 0 87. 71 7.71 
5.¢ 2458. 0 SI 2 ” ” SS. ON SS. OS 
14.8 tor 1] 110. 43 | > 
219 14.2 1207. ¢ il 110. 40 110. ¢ 
9] 527 8 1410.7 1 120.8 } 
2 g9.8 1410.7 TT) 13 130. 98 130. 92 
219 70 17 14 137. 57 — 
219 NN 3 1 14 137. 57 j set 
. 2 196. 2 g ( 1.11 l 144.0 | . 
5 219 172.8 1618. 0 2. 52 ] 144.0 p 144.0 
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are given in table 11. These relations can be used 
as a means of interpolating below 550° C for smooth 
values of the relative heat contents of compositions 
on which no measurements were made. 

The phase diagram shown in figure 3 indicates a 
first-order phase transition at a eutectoid tempera- 
ture of approximately 550° C. Though the slopes 
in figure 2 represent only average heat capacities 
over the finite temperature intervals of measurement, 
it is apparent that the heat capacity of each hydride 
sample reaches a maximum value within the range 
540° to 600° C, corresponding to this phase transi- 
tion. However, it will be observed that the partic- 
ularly high heat capacities are spread over a tem- 
perature range, a fact to be interpreted as due to 
considerable lack of phase equilibrium at these tem- 
peratures. This is not an unusual situation in the 
case of all-solid systems. The interpretation is sup- 
ported by the unusually poor reproducibility of du- 
plicate heat values at 560° and 575° C and the 
hysteresis found at 600° C, as may be noted from 
the last two columns of tables 4 to 8. 

As the temperature rises above 600° C, the heat 
capacities of samples 3 to 6 decrease rather abruptly 
to smaller values at temperatures that differ for the 
four samples. This fact also can be given a quali- 
tative interpretation in terms of the phase diagram, 
for above a certain temperature each hydride sample 
investigated should become all 8 phase, and there- 
fore at higher temperatures there would no longer 
be a contribution to the over-all, or observed, heat 
capacity by latent heat of phase reproportionation. 





TABLE 9. Mean observed heat capacities of zirconium metal 








= 
bserved (NBS N=6(ZrH, .. 
Deviations of smoott Valu f other - 
investigat ) N: Sample number used in tt : 
Mean heat poper (see Table Lathan 
Temperature capacity ‘ F li 
interval this Red Y 
invest Bul Cough- | Skinne mond Scott 2 
gation letir lin and ( nd {Ss} ’ 
7619 King [5 I N= 4(77- 
; cal | - g N= 3(ZrHo 556 
( deg ¢ y. 
0 to 100... 0. 0691 2 +] ; 19 2 \ 
100 to 200 O717 ; +4 2 1.5 3 
200 to 300 O772 6 +1 +9 s 2 j rs N ZrH aad 
300 to 400 O74 0 4 ( is 6 , 4 £ : 
100 to SOO 0793 +2 +3 l +Y 7 
500 to 600 OS809 15 3 1] a > 
600 to 700 ORS1 ( 2 L10 ) — 
700 to 800 OSH0) + 1 7 , 
SOO to Si 0932 ( 6 1.7 11 
S50 to 863 oyu ll lt +2 15 Zr 
Heat of trans ) 
formation v.38 - re tli rll 
is. 
So. tO WOO  oa7 18 4 " 
TABLE 10. Vean obse ed heat capaciti« of five amples of | 
irconium h jaride 
I irt Sal é ee | ] Sa t ~ Sa 
324 ‘ ’ 1 
7 caid 
( ( dé ( ( ( t 
- Figure 2. Heat content versus temperature, zirconiun 
0 to 100 0. 0738 0. O762 Ts r s RO3 
100 to 206 ot 4 OAS) ] s ‘ era ¢ ov / 
20 iM OSUuU 1028 ‘ l ] 
S00 t nM ws L087 Jt 1273 
400 to 500 1174 1244 1302 { 1464 
00 to 54 12 With due consideration to those few cases where 
iiaia p «1402 151 1808 phase equilibrium was obviously not reached during 
( f HAS | » : ° . 
the measurements of heat content, it is possible to 4 
e—— » = derive numerous quantitative energy and phase re- 
60 to 575.......|¢ . 408 M 13 lations in the zirconium-hydrogen system from the | 
M 6 378 20 precise thermal data reported in this paper. Such ’ 
iene - = - derivations are described in detail in a current pub- { 
Loot 2252 lication [2]. 
64 " } 1194 LISS 2176 - 
; 
” f 
wo 12 9 1K {0 
ree | si) 
The zirconium metal and its analysis were re- 
M 1165 12 : ceived from the U.S Department of the Interior, 7] 1 
0 to 900 0 sureau of Mines, Albany, Oregon, through the 


courtesy of M. L. Wright. Several members of the 
T 1] E , . oe National Bureau of Standards staff performed Val- Q] 
ABLI m ~quations of moothed heat conter ; j COMDOSiI- . . , i * 
ey ay wer $y f uable services during the work. J. K. Russell made 


lion of conium hudride pn to N°" ; 
a number of the heat-content measurements on two 
of the hydride samples. Martha Darr carried out 1 ] 
ties . aii ai all the spectrochemical analyses of the Sylvania hy- 
dride samples, and R. A. Paulson analyzed the ar- 
remper { B conlum metal and duplicates of the Sylvania hydride 
ature, samples for carbon and hydrogen. 
o 
100 637 114.0 
an 13 f et 
a0 2041.8 604.9 
sn) 2840 4 1039. 7 
‘ 3794. 2 1410. 0 
) 42 0 17 s 
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Unitorm Transient Error ' 


Edith L. R 


The economy of transient error is discussed. 


. Corliss 


Equations describing error in power- 


level measurements of transients can be used to compute the design of analyzers so as to 


distribute transient error in a way compatible 


In addi- 


with experimental requirements. 


tion, consideration of a limiting power-discrimination factor provides a measure of the 
largest number of bandpass filters that can be overlapped on adjacent channels to yield 


meaningful information about a rapidly changing signal. 


A scanning filter can be compared 


with a set of bandpass filters by calculating the effective number of overlapped bandpass 


filters providing the same resolution as the sez 
Restrictions on rapid signal analysis also are considered for the 


autocorrelation analyzers. 


inning filter. The results are applicable to 


special case in which power-level discrimination is limited by noise. 


1. Introduction 


Measuring systems are limited by natural laws, 
and several of these guarantee that errors are inevi- 
table. However, the distribution of errors can be 
controlled by the experimenter. One can choose 
experimental conditions that minimize errors in a 
region of particular interest; one can choose alterna- 
tively to spread the errors out as thinly as possible. 

The transient properties of a measuring system 
become dominant whenever there is a restriction in 
time. The constraint may arise either because a 
signal is to be analyzed by a rapid process or because 
a signal may be changing rapidly in time. 

In order to minimize transient error, experimental 
conditions must be chosen on a quantitative basis. 
Some prior knowledge of the signal to be analyzed 
is required. For the purpose of choosing the condi- 
tions of measurement, the problem of analysis can 
be restated as follows: Given an ensemble of signals 
characterized by a range of fluctuations in power 
over a range of intervals in time, one certain signal 
is to be distinguished from all other signals in the 
group. 

A signal is discrete in time. The fluctuations in 
power during the existence of the signal convey its 
identity. If the power-level fluctuations are to be 
indicated within a certain degree of error, the time 
required for the analyzer to build up to within a 
corresponding fraction of its steady-state value must 
not exceed the shortest time interval over which a 
significant fluctuation of power level occurs in the 
component of the signal being analyzed. Thus, the 
rate at which the indication can build up or decay 
in the analyzer must exceed a limiting value. It is 
useful to describe this rate by reference to the 
response time of the analyzer, i. e., the time required 
to reach substantially steady-state indication after 
an abrupt change in signal. 

Regardless of whether a frequency or a time- 
correlation analysis of the signal is undertaken, the 
results are observed as output power levels. Hence 
the ability to resolve a change in power level limits 
the discrimination that can be achieved for changes 
in frequency or intervals of time. The relation be- 


ee:-_-._ 

'This work was carried out as part of a program of pasic instrumentation 
research and development cooperatively supported by Office of Naval Research, 
Atomic Energy Commission, Air Force Office of Scientific Research, and the 
National Bureau of Standards 
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tween the various limitations on the performance of 
an analyzer can be illustrated by means of a three- 
dimensional system (see fig. 1) in which the Cartesian 
coordinates are frequency, time, and either power 
or a function of power. The transient error then 
becomes a least volume in this space, within which 
no information about the signal can be found. The 
size of this least volume is limited by the natural 
behavior of analyzers, but its shape can be distorted 
to suit the needs of the experimenter. 
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Figure 1. Three-dimensional signal space. 


Resolutions in frequency and time are implicit in 
the Q (figure of merit)? of a frequency analyzer or, 
correspondingly, in the integration time of a corre- 
lation analyzer. To minimize the errors in an anal- 
ysis, either Q or the integration time must be kept 
as high as possible consistent with the necessary 
resolution in power level. 

Consider the rate at which power builds up in a 
linear series-resonant system, in a simple filter or a 
galvanometer. A family of indicial response curves 
for such a system is shown in figure 2, in which the 
coordinates have been scaled to relative values. The 
amplitude response has been scaled with the steady- 
state value as a unit; unity on the time scale is the 
natural, undamped period of the system. Curves in 
the family are plotted for different values of Q. 


2 The ratio of the energy-storage capacity of a resonating system to its energy . 
dissipation rate is known as its “‘figure of merit’’, or Q. 
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Indicial response of a series-resonant system 
subjected to a step function. 


FIGURE 2. 


This shows the fraction of steady-state response attained as a function of the 
time following the impulse, scaled in units of the natural period of the system. 


The analyzer may be considered to have reached 
a sufficient fraction of steady-state indication when 
it is just capable of resolving all unlike signals, 
whereas all like signals are indicated as being the 
same. Attempts either to speed up the response or 
increase selectivity are obtained at the expense of 
decreased Q or increased response time, without pro- 
viding any additional information about the signal. 
Lengthening the response time, moreover, impairs 
the ability of the analyzer to indicate rapid power 
fluctuations, so that some features of the signal can 
be lost in this way. 


2. Transient Limitations of a Simple 
Analyzer 


The quantitative extent to which power-level res- 
olution affects resolution along the coordinates of 
frequency and time can be found from the natural 
properties of the analyzer. The discussion to follow 
will be based on the common type of system whose 
behavior is described by a second-degree differential 
equation with constant coefficients, but its applica- 
tion to other systems will be indicated. 

For example, consider the use of a linear series- 
resonant system as a tuned filter to indicate how a 
particular component of a complex signal fluctuates 
with time. If z represents amplitude, the transient 
response of the filter in terms of its center frequency, 
fo, and figure of merit Q, is 


r itd CF ie + i? 4Q2)t Les fray Si (5 /1Q? t). 1) 

Let the observation interval, Ar, be sufficiently 
long compared with 75, the undamped natural period 
of the filter, so that only the envelope of its response 
need be considered. 

The analyzer should be able to change its indication 
by at least the smallest detectable power increment 
during a specified time after onset or removal of the 
input signal. The power built up in the filter is 


proportional to |x|’. If Wo is the initial power in the 
filter, at any time after cessation of the signal th 
power W, is given by 


2ef, 22 ¢t 


—-—t 
W,=We ° 


Woe @ To 


Let AW=(i—e-*)W be the smallest discernibj. 
change of power level in the filter; i. e., the power 
must decay by a factor e~* for a change to be ob. 


served. If signal cessation is to be discerned withip 
time 
22 Ar 
W.—Wi,=Wi(1—e  °% 79 >(1—e-2)W, 
so that 


: (= =) 
e QT, <e¢, 
and the requirement is 


Q<28At_2nfo, P 
aT, a 

if the filter is to indicate the occurrence of a change in 
signal during the interval Ar. 

Equation (2) has been derived upon the basis of 
the finest resolution in power levels that can be 
achieved with the analyzer in the time interval Ar, 
Given that Ar is the briefest duration of any signal 
component to be observed, eq (2) also shows how the 
Q of a filter must be adjusted as a function of its 
tuning in order to distribute the transient error uni- 
formly over the frequency range. In addition, be- 
cause eq (2) deals with limits of resolution, it pro- 
vides a measure of the largest number of bandpass 
filters that can be overlapped on adjacent channels 
to yield meaningful information about a rapidly 
changing signal. 

Over a range of frequencies Af=f,—/; surrounding 
f,, the frequency to which the filter is tuned, the 
‘change in filter response does not exceed the smallest 
detectable power change. Consider the relative 
transmission characteristic of a filter having its Q 
given by eq (2) (see fig. 3). The relative power re 
sponse of the filter is given by 


» W, 
VW f ca1Po | ro Tw 
© LIA ECU—F fo)? 
The precise result of solving for the discrimination 
limits in terms of f, and f/; comes out as 


f —f; a , . a? 


fi. —xfAr V“ l6x?f2Ar? 


IV 


but if AW is small compared to W, (i.e., 
discrimination can be made) then 


fine-power 


9 


a 
fu-fi=Af2z5— 
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Figure 3. Relative power response of a series-resonant system 


as a function of steady-state driving frequency. 


3. Analyzers Made Up of Overlapped 
Channels 


For a set of filters overlapped at the limits of the 
power resolution, the limiting number of filters 
yielding meaningful information between f, and f,, 
the high- and low-frequency range limits, respectively, 
is given by 


f1—fy 2xAr 
nsf at info 


(The limits yield n,—1 rather than n,; because one-half 
of a filter extends beyond each end of the frequency 
range as defined here.) 

In practice, it is usual to overlap adjacent filters 
at their half-power points, so that their integrated 
response has a flat characteristic. For this condition, 
Af=f,/Q, and the number of filters becomes 


n mt Si—fo_ @Fi—Io) _ 


: _24Ar(fi—fo) 
. Af Jo a 


The ratio of this number to the maximum usable 
number of filters (overlapped at their limits of resolu- 
tion) turns out to depend only on power resolution, 
Le., 

No—1 


14 
ge 
ne l 


If we are willing to take a “yes-no”’ type of analysis 
from the filter, then for a set of filters overlapped at 
their half-power points we can choose as a com- 
patible condition that a=In 2, i.e., the power declines 
at least to one-half its initial value during the obser- 
vation interval Ar, and we get as the number of 
filters for this case: 


2rAr( fi—ho) 


Np— | 
. In 2 
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Under these conditions, n, invites comparison with 
the number of “‘bits’”’ available in a given bandwidth 
during the time interval Ar.’ 


4. Scanning Analyzers 


The process of scanning makes it possible to study 
the spectrum of a signal by means of single channel 
in place of the numerous channels needed when fixed- 
bandpass filters are operated in parallel. When one 
scans by sweeping the tuning of a filter over the 
frequency range to be covered, the filter shows a 
rather complicated type of response.t However, one 
can make a quantitative study of the effects of 
scanning by use of simple theory, if one recognizes 
that another multiplicative factor will be necessary 
to take accurate account of the spreading of the 
filter’s transient response due to the effects of scan- 
ning. The spreading occurs because the incoming 
signal beats with the ringing of the filter. 

If one sweeps through a band of frequencies in- 
cluded between the limits f, (high-frequency) and f, 
(low-frequency) in the time interval 7, the length of 
time available to build up the power response of the 
filter at any instantaneous tuning is generally much 
smaller than the whole sweep interval. Call the 
observation interval for each effective single filter Ar, 
as for discussion of the fixed-bandpass filter. This 
quantity is directly related to the length of time that 
it takes the center tuning of the filter to change 
through a range of frequencies equal to the band- 
width between its half-power points. In this simpli- 
fied discussion we shall take this relationship as 
unity, remembering that it is an approximation. 

Various types of scanning functions can be used, 
depending upon the nature of the selectivity charac- 
teristic desired. If one wishes to have a constant Q 
(constant percentage bandwidth), the condition that 
the transient error is to be uniform throughout the 
range of frequencies swept over can be satisfied by 
scanning at a rate proportional to the square of the 
instantaneous center tuning of the filter. When one 
wishes to scan through the frequency range with a 
filter whose bandwidth is a constant number of cycles 
per second, so that the Q is directly proportional to 
the center tuning of the filter, the rate of sweep 
should be constant independent of frequency. One 
type of scanning filter sweeps logarithmically over 
its frequency range; when it is adjusted to uniform 
transient response, its instantaneous Q is propor- 
tional to the square root of its instantaneous tuning. 
Thus it is intermediate in its properties between 
constant-frequency bandwidth and _ constant-per- 
centage bandwidth. 





3R..V. L. Hartley, Transmission of information, Bell System Tech. J. 7, 535 
1928). 
: 4F. M. Lewis, Vibration during acceleration through a critical speed, Trans. 
Am. Soe. Mech. Engr. 54, 253 (1932); G. Hok, Response of linear resonant system 
to excitation of a frequency varying linearly with time, J. Appl. Phy. 19, 242 
(1948); N. F. Barber and F. Urseil, The response of a resonant system to a gliding 
tone, Phil. Mag. 39, 345 (1948); J. Marique, The response of RLC resonant circuits 
to EMF of sawtooth varying frequency, Proc. Inst. Radio Engrs. 40, 945 (1952); 
M. C. Herrero, Resonance phenomena in time-varying circuits, Technical Report 
No. 69, (Electronics Research Laboratory, Stanford University, 1953); H. W. 
Batten, R. A. Jorgensen, A. B. MacNee, and W. W. Peterson, The response of & 
Panoramic Analyzer to CW and pulse signals, Proc. Inst. Radio Engrs. 42, 948 
(1954); S. S. L. Chang, On the filter problem of the power-spectrum analyzer, 
Proc. Inst. Radio Engrs. 42, 1278 (1954), 








The relations between @ and the sweep rate follow 


from the restrictions of eq (2) on minimum dwell 


time: 


aQ 


3 (3) 
2rfo. 


Ar 


Describing the instantaneous scanning rate by the 


expression 


h( fo \dfy- dt . 


an approximation to Ar is obtained by integrating 
the scanning function between the half-power points 
of the filter. (To an approximation good to 12% 
percent, even for Q as low as 1, the limits are f 
+ f,/2Q.) Thus, 


f ~f po 
Ar ~~ Afydfo- 
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Expanding this integral into Taylor’s series: 
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So that, to a close approximation: 
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Looking back at the original integral, evidently 


h(fo) 
and, substituting into the series 
Ar fo fo) +:( ; Jo )h "| fo) higher terms. 
Q 2 
But, from eq (3) 
it ace lf fo \? , 
Ae atl . htfo)+ ( : yn’ f., ete. 
2rfn Q * 3\2Q 
If Q is of moderate size or higher, then, very nearly 
. a 
h( fo) De fo? (4) 


The precise way in which Q depends on the total 
scan time, the scanning function, and the total 
bandwidth covered can be obtained by substituting 
for eq (4) in the expression 

h(fo) 


df,=dt. 
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Assigning the desired variation of Q with frequency, 
and integrating for the total scan time 7 between the 
band limits f, “and ti, We get 


“Se a Cd 
— 1 Or for J (9 


The magnitude of Q can be found by solving this 
equation, 

Usually + is the reciprocal of the repetition rate 
of the scanning filter. The value of @ is a measure 
of the fineness of power-level discrimination that cap 
be achieved. (Notice that eq (5) provides a means 
for assigning to a any desired dependence upon fre. 
quency, so that transient power discrimination can 
be distributed at will over the frequency range.) 

The functioning of a scanning filter can be com. 
pared with a parallel channel set of bandpass filters 
by calculating the effective number of bandpass 
filters that will provide the same resolution as the 


scanning filter. The filters may be taken to be 
overlapped at their half-power points, and are 
calculated to cover the frequency range scanned 


same power-level resolution and figure of 
those available from the scanning filter, 

narrow frequency differential, df, in the 
the “concentration” of 


with the 
merit 

For a 
tuning of the scanning filter, 
filters in the range is given by 


as 


F 
dn d, 2 
Al 


where Af is the width between the half-power points 
of the scanning filter at the instantaneous center 
frequency, fo. The figure of merit, Q, of course 
fo/Af. Sine e one-half of a filter exte nds beyond each 
limit of the frequency range defined, integrating 
over the range gives one less than the equivalent 
number of parallel-channel filters that have their 
center tuning extending over the range, or 


Is, 


5. Uniform Distribution of Transient Error 


Table 1 summarizes the results of applying the 
foregoing discussion to calculate Q, the figure of 
merit, &, the scanning-rate parameter, and Mm», the 
equivalent number of filters overlapped at their half- 
power points, in terms of the scanning time, 7, the 
scanning bandwidth limits f, and f,, and the atten- 
uation parameter, These calculations are based 
on uniform distribution of transient power-level dis- 
crimination, 1. e., w@ is taken as a constant. 

To facilitate adjustment of the response time of 
an analyzer to the least duration of any signal com- 
ponent to be observed, the tabulation also includes 
Ar the observation time interval for any single filter. 

Although the spreading due to ringing in the scan- 
ning filter is neglected, the relations shown in table 1 
indicate quantitatively how the transient error and 
scanning characteristics affect the resolution of an 
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Tapie l. Re lationships among variables describing the performance of scanning and fired bandpass filters in which error in wransient 
power-level discrimination is uniformly disiributed 


©, Figure of merit 


Seanning rate for selectivity 
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; ; > 
Quadratic scan, df/dt=kf¢ < = 
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analvzer. The spreading due to ringing can be 
taken eare of through the power-level resolution 
factor @. 

As a numerical illustration of how the factors oper- 
ate. consider what sort of analysis can be accom- 
plished over the frequency range from 100 to 10,000 
oa in 0.1 sec, if we are satisfied with a 10- percent 
power-level error. 

From a bank of parallel bandpass filters, over- 
lapped at their half-power points, you can get mean- 
ingful results from 6.26 10* filters, with a Q ranging 
from 6.28 & 10? at 100 eps to 6.28 10* at 10,000 eps. 

If vou scan uniformly at 10 sweeps/sec, the effec- 
tive number of filters is markedly decreased because, 
of course, much less response time is available for 
each filter; in fact each filter must now respond with 
90 percent of its final indication in about 4 msec. 
This sort of scan is equivalent to 252 filters, over- 
lapped at their half-power points, and the available 
a ranges from a (Y of 2.52 at 100 cps to 252 

t 10,000 eps. 

A logarithmic scan, which allows somewhat longer 
dwell time at low frequencies, is equivalent to only 
somewhat fewer filters, 211 being the number under 
the conditions of the example given above. How- 
ever, the Y at the low-frequency limit is 12.2, so that 
there has been some improvement in low-frequency 
resolution. This has been obtained at some loss in 
high-frequency selectivity; at 10,000 eps the @ is 122. 

Where constant percentage bandwidth is desired, 
the sean rate is proportional to the square of the fre- 
quency. A correspondingly larger part of the time 
of observation is devoted to the low-frequency region. 
The limits converge on a Q of 25.2 independeni of 
frequency, covering the range with 116 filters. 


6. Limitation by Noise 


The noise level associated with the signal provides 
a practical limit to the power-level resolution that 
can be obtained. Any further increase in resolution 
merely results in observing the noise, so that like 
signals may appear to differ. 

An especially interesting expression results from 
considering the least power change observable to be 
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limited by the noise present with the signal. If S is 
the signal power and N is the noise power, we take 
s a limit for detection 


AS>N., 
Considered as a fractional power change, 


AS _W—W,_,_-45_N 


S+N W, =—S+N 
When S/N>1, a—>N/S and 


afar>>(%) (6) 


and, correspondingly, 


V< ice ( a ) 2rfyAr ( < (7) 


A higher Q can thus be used when there is a good 
signal-to-noise ratio. Equation (6) has a provoc aiive 
resemblance to the Heisenberg uncertainty principle. 
The resemblance is real: both are derived by related 
mathematies, and both are based on the idea of 
least discernible power change. 

A special case of practical importance is the detec- 
tion of a sinusoidal signal in the presence of ‘white 
noise.”’ (Noise energy per unit bandwidth in cycles 
per second is constant.) For the filter calibrated at 
iis peak response by means of a sinusoidal signal of 
power W,, the same indication will be obtained with 
a white noise of /, units per cycle per second, given by 


QW. 
rf 
This expression is obtained by integrating over the 


filter response.® 


5D. Bierens de Haan, Nouvelles Tables D'Intégrales Définies, eq. 6, p. 47 
(Edition 1867 reprinted 1939). 
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Detection of a sinusoidal signal of power Sp (at the 
center frequency to which the filter is tuned) in the 
presence of a noise distribution of n; units per cycle 
per second proves to require a least time interval 
Aro, Which is independent of Q. 


any 
Aty> 255. 


Narrowing the filter does decrease the amount of 
noise passed, but it lengthens the response time by a 
proportionate amount. 


7. Transfer of Results to Correlation 
Techniques 


The foregoing discussion of selectivity and response 
time can be transferred directly to the design of 
correlation analyzers. In a paper on synchronous 
detectors ° Selgin derives a relationship between the 
Research 


* Paul Selgin, Harmonic of the synchronous rectifier, J. 


NBS 47, 427 (1951) R P2267. 


output 





integration time of a correlation analyzer 
of a filter exhibiting equal selectivity, 


and the Q 


Q= 7 


where 7 is the period of the synchronous detector 
and + is the R-C integration time constant. This 
relation permits calculation of the transient error 
distribution in correlation analyzers from the results 
in this paper. 


The author expresses appreciation to Chester ]. 
Page, who provided much-needed advice on this 
work, and to whom the generalized mathematical 
treatment deriving the scanning function from an 
arbitrary Q is due. 


WASHINGTON, July 26, 1956. 
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Optical I—Bench Method of Measuring Longitudinal 


Spherical Aberration 


Francis E. Washer 


A method of measuring longitudinal spherical and chromatic aberration by visual 


means is described. 


equipped with nodal slide and an angle-measuring telescope. 
the method is presented, together with a brief description of the apparatus used. 


The method employs an especially constructed optical T-bench 


The underlying theory of 
The 


results of measurement on a few typical lenses and an analysis of the sources of error are 


included, 
1. Introduction 


The determination of longitudinal spherical aberra- 
tion is frequently required during the course of the 
measurement of the optical characteristics of photo- 
eraphic objectives. While many methods have 
been developed that are capable of yielding this 
information, it, nonetheless, seems worth while to 
describe a method that has been in use at the 

' National Bureau of Standards since 1943. This is 
an indirect method, in that the lens is used to colli- 
mate the light from a point source located in the 
focal plane, and the angular deviations of adjacent 
beams from parallelism are measured. The de- 
termination of the longitudinal spherical aberration 
is based on the analysis of these angular deviations. 
At the time of inception it was necessary to develop 
asimple, rapid method of evaluating the longitudinal 
spherical aberration of reflector sight lenses.‘ In 
connection with the development of this rapid 
method, the present more accurate procedure was 
developed for the purpose of checking the results. 

Early work was done with relatively simple appara- 
tus that was rather unstable. The apparent worth of 
the method was, however, sufficiently high as to 
warrant the construction of new equipment in order 
that a higher degree of accuracy could be obtained. 
The resulting equipment consisting of the optical 
T-bench and a new nodal slide is shown in figures 


sand 4, 
2. Theory of the Method ! 


When a lens is affected with longitudinal spherical 
aberration, the effect can be described most simply 
by considering it as change in positions of the focal 
point for the various annular zones of the lens. 
For an ideal lens the focus is a point on the axis in the 
image space through which pass all rays of light 
that are initially parallel to the axis in the object 
space prior to incidence on the front surface of the 
lens. If the lens is reversed and a point source of 
light placed at the focus, the light incident on the 
back surface of the lens will emerge as a parallel 
beam of light concentric with the lens axis. This 


condition is shown schematically in figure 1, where 


—_ --——- 
F. E. Washer, An instrument for measuring longitudinal spherical aberration 
flenses, J. Research NBS 43, 137 (1949) R P2015 
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Figure 1. Schematic drawing of a lens showing light incident 
on the lens from a point source on the azis at O. 


For O located at the focus, the emergent beam is parallel to the lens axis. When 
O moves by amount, Af, to O’, the emergent beam is convergent, and a ray 
located at zone height, h, is deviated through the angle, «. 


light from the point source, O, is incident on the 
lens at height A and emerges parallel to the axis OX. 
Let the point source, O, be moved by amount Af, to 
O’. If O’ is farther from the lens than O, the 
emergent beam is convergent. For a ray incident 
at height h, the emergent ray will deviate from paral- 
lelism with the axis by the angle e. It can be shown 
that 

Afh 


“Fray 
where f is the equivalent focal length of the lens. 


For known values of ¢, the value of Af can be de- 
termined to a close approximation from the relation 


f(,. Af 7 
apa (14 ;) ( 


For Af, very small compared with f, the above 
relation may be written 


Af=_ 


(1) 


(3) 


It is clear that only a 1-percent error is introduced 
into the determination of Af for a value of Af=0.01f. 
This can be further reduced by inserting the approxi- 
mate value of Af obtained from eq (3) into eq (2) and 
redetermining Af. However, for most practical pur- 
poses the value of Af computed from eq (3) is satis- 
factory. 

If the lens shown in figure 1 is affected with longi- 
tudinal spherical aberration, rays of light from the 
point source at-0 incident on the lens at varying zone 
heights, A, will deviate from parallelism with the axis 








by varying amounts, «. If the various values of ¢ are 
plotted as a function of h, a curve such as No. 1 in 
figure 2, A, will result. For this curve, the variation 
of e with A, arises from a combination of two effects: 
(1) Out-of-focus effect, and (2) angular spherical 
aberration or parallax resulting from longitudinal 
spherical aberration of the lens. The two effects can 
be separated quite easily, as out-of-focus effect pro- 
duces a variation of € proportional to h, which plots 
as a straight line of constant slope passing through 
the zero point of the graph. This is shown by curve 
3 in figure 2, A. The nonlinear portion of the curve 
arises from angular spherical aberration. On sub- 
tracting the abscissas of curve 3 from those of curve 
| for corresponding values of h, the resulting values 
of « are those produced by the angular spherical 
aberration of the lens; this is shown as curve 2 in 





placed from the zero line of abscissas; this indicates 
that the target did not lie in the plane of best foeys 
for any zone of the lens. When the adjusted values 
of e (plotted as curve 2 in fig. 2, A) are used in deter. 
mining the values of Af, the result is as shown in 
curve 2 of figure 2, B; for this condition of focus, the 
target would lie in the focal plane of the zone fo, 
h=40 mm. 


3. Method of Measurement 


In order to measure the longitudinal spheric} 
aberration of a lens by measuring the variations jp 
parallax in the emergent beam, it is necessary to use 
equipment capable of making accurate measure. 
ments of parallax angle and zone radius h. The Op- 
tical T-bench, shown in figure 3, was developed aj 
this Bureau for use in making such measurements 
It consists of two optical benches mounted at right 














angles to one another in the form of a T on a strong 
metal table. The lens under test is mounted on the 
nodal slide, as shown, with the target and illumi. 
nating system mounted on a separate bench slide as 
indicated. The viewing telescope is mounted on g 
movable slide located on the T-bench. The tele 
scope is so arranged that its line of sight is approxi- 
mately normal to the bench ways upon which its 
base slide can move. The pointing of the telescope 
is done with micrometers near the ocular end, whieh 
produce a small rotation about effective horizontal 
and vertical axes located in the flexure plate mount- 
ing near the objective end of the telescope. The 
change in pointing caused by a movement of the 
horizontal micrometer of one scale division is approxi 
mately 5.8 see of are. (The accurate calibration 
factor is 0.0283 milliradian per scale division.) In 
reading the micrometer, the value of the setting to 
the nearest one-tenth of a scale division can readily 
be made. , 
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Figure 2. Parallax angle ¢ versus h, and Af versus h, for a 
typical lens 
4. Curve 1 shows the variation of « with h when the lens has longitudinal 
spherical aberration and the target is not at the position of best focu Curve 3 
shows the variation of e arising from out-of-focus effect alone of Af—0.67 mm, 
Curve 2 shows the angular spherical aberration alone (the difference between 
land 3 
B, Curve 1 shows values of Af computed before making the correction for out- 


of-focus effect. Curve 2 shows values of Af after making correction for an out-of 
focus effect of 0.67 mm. Curve 2 is also the longitudinal spherical aberration 
plotted with respect to the plane of best focus for zone height, A= 40mm 


When the value of f is known and the values of « 
as a function of A have been measured for the condi- 
tion of point source located on the lens axis and in the 
vicinity of the focal point, the values of Af can be 
computed with the aid of eq (3). When the values 
of Af are plotted as a function of /, the resulting curve 
shows the longitudinal spherical aberration. If the 
point source is not at the true focal point, the entire 
curve is displaced along the axis by an amount that 
is readily obtained from the graph or from the 
computed values of Af. 

In figure 2, B, curve 1 shows the values of Af ob- 
tained using the same values of e shown in curve | of 
figure 2, A. It is clear that the entire curve is dis- 


Arrangement of apparatus on optical T-bench for 
determination of longitudinal spherical aberration. 


FIGURE 3. 


lhe lens under test is mounted in the nodal slide in the center of the picture 
rhe target reticle and illuminating system are shown on the movable slide ™ 
the right of the nodal slide. The viewing telescope, equipped with diaphragm 
over the objective, is shown mounted on the transverse nodal slide in the left 
portion of the photograph, 
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In preparing to measure the longitudinal spherical 
aberration, the lens is so positioned on the nodal slide 
that its rear nodal point coincides with the axis of 
rotation of the nodal slide. The axis of the lens 
under test is brought into parallelism with the bench 
wavs, and the illuminated target 1S placed at the rear 
focal point of the lens. The positioning of the target 
at the focal point of the lens is done with the aid of 
the viewing telescope, which is itself set at infinity 
focus. The level of lens axis, target, and center of 
the viewing telescope object are adjusted so that all 
lie in approximately the same plane, which is parallel 
to the plane defined by the surfaces of the two bench 
ways. 

The telescope objective is now covered with a dia- 
phragm having a small central aperture, which is 
centered with respect to the telescope objective. 
When an observer looks through the telescope, he 
sees the image of an illuminated target (either a cross 
or pinhole pattern) superimposed on the cross hairs 
in the ocular of the viewing telescope. When the 
telescope is moved transversely by moving the slide 
upon which it is mounted, the image of the reticle 
will either stay on the cross hairs of the ocular or 
move toward one side of the visual field or the other. 
In the event of movement of the image, it can be 
brought again into coincidence with the cross hairs 
by appropriate movement of the transverse microm- 
eter. The angular movement thus noted is a measure 
of the change in direction of a small circular beam 
of light emergent from a given area of the lens under 
test from the direction of a similar beam from another 
given area of the lens and lying along the same 
diameter. 

lt isa simple matter to select a given size of dia- 
phragm opening, small in comparison with the area of 
the lens under test, and to select a series of steps so 
that measurements of the deviation of these succes- 
sive small emergent beams can be made for an entire 
diameter ot the lens If the lens has been properly 
focused, and the telescope so moves that it always 
remains parallel to itself, the angular deviations noted 
arise from parallax (or angular spherical aberration 
in the lens under test. The bench ways are so made 
that the telescope does maintain this condition of 
remaining parallel to itself at the various positions 
30 that possible errors from this source can be re- 
garded as negligible 

The scale on the bench carrving the telescope is 
sed in making the successive settings at selected 
intervals These settings are recorded, and the dif- 
lerence of successive bench scale readings from that 
one in the central position is recorded as h, the zone 
height. The micrometer readings for the successive 
pontings are taken for each value of h, and their 
differences from that one in the central position are 
recorded as e's. It is, of course, necessary to make 
icareful determination of the central position. This 
s usually done at the start of the experiment by 
hoting the total range of transverse movement, 
throughout which the image of the illuminated tar- 
get remains visible in the telescope ocular to the ob- 
server. These end points are noted, and a point 
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midway between them is accepted as the center of 
coordinates. 

Although it is believed that the possible deviations 
of the line of sight with movement of the telescope on 
its slide are negligible, it is presently the custom to 
avoid all possibility of error from this source by using 
an arrangement such as that shown in figure 4. 
Here the telescope slide remains fixed, and the scan- 
ning across the aperture of the lens under test is 
accomplished by moving a slide carrying an appro- 
priately mounted pentaprism. The deviation of the 
beam of light by the pentaprism is not changed by 
small rotations of the pentaprism about its vertical 
axis, Which might be caused by small rotations of 
the carrying slide as it moves along the bench ways. 

It is clear, then, that with this arrangement of 
apparatus it is possible to measure the values of € as 
a function of A in the emergent beam of collimated 
light from a lens. These values can be used to de- 
termine the values of the longitudinal spherical aber- 
ravion with the aid of eq (3). 





FiGuRE 4. 


Arrangement of apparatus on optical T-bench for 
determination of longitudinal spherical aberration. 


Chis arrangement differs from that shown in figure 3 in that a pentaprism has 
been placed on a movable slide on the transverse bench to eliminate transverse 
motion of the telescope 


4. Results of Measurement 


This method of measuring longitudinal spherical 
aberration has been successfully used for several 
vears at the Bureau for lenses having a wide range 
of focal lengths. In the succeeding sections, the 
results for a few typical lenses are given. 


4.1. Photographic Objective Having a Focal Length 
of 150 Millimeters 


a. Effect of Target Displacement Along Axis 


Although there is no occasion to doubt the validity 
of the relations expressed in eq (1) and (8), it is 
nonetheless of interest to report the results of an 
experiment that confirms these relations and simul- 
taneously puts to rest any lingering doubts concern- 
ing the proper sign of Af. In addition, considerable 








information is gained on the probable accuracy and 
limitations of this method of determining longitudinal 
spherical aberration. 

The lens selected for measurement was a photo- 
graphic objective, designated lens 2, having a nom- 
inal focal length of 6 in. and a maximum aperture of 
4/6.3. The measured value of the equivalent focal 
length was 152.4 mm, which means that the maxi- 
mum value of & (the half-diameter of the aperture) 
was 12.1 mm. The arrangement of apparatus was 
as shown in figure 4. The diameter of the diaphragm 
opening placed over the telescope was 2.0 mm. 

The viewing telescope was first focused for infinity 
for the condition of no diaphragm over its objective. 
The lens and target reticle were then accurately 
alined, so that the axis of the lens was parallel to 
the bench ways, and the center of the reticle was on 
the optical axis of the lens. The slide carrying the 
target system was moved along the bench until the 
reticle was in the plane of best focus of the lens under 
test as viewed through the telescope. The 2.0-mm 
diaphragm was then placed over the telescope and 
the center of coordinates determined. The results of 
measurement are shown in table 1 for five settings of 
the target. In this table, the target settings with 
respect to the lens bench scale are shown at the top; 
the highest numbers indicate the positions of the 
target nearest to the lens. The nominal values of 
Af (hereafter designated 6f to distinguish from the 
measured values of Af), with respect to the plane 
of best focus, are shown on the next line of the table; 
these values are approximate as the bench scale is 
read with a vernier to the nearest tenth millimeter, 
and no attempt was made to make these settings to 
TABLE 1. Variation of parallax angle, ¢, with transverse 
selling of telescope for five positions of the target with respect 
to lens 2, which has a focal length of 152.4 mm 


The column headed 6f=0.0 corresponds to the position of best focus; negative 
values of 6/ indicate the target is farther from the lens 
Target settings on lens bench, mm 
1 wf 91.1 1. 72.1 
Tele Scop" 
positions on 
bench, mm Nominal values of éf, mm 
—1.0 0.5 0.0 0.5 1.0 
Measured values of ¢, in scale divisions 
Left 
® 25.3 
28. 2 4.3 7.3 417.8 16.9 3.4 
30. 2 82. 2 73.3 62.9 “2 1} 
32.2 83.34 75.3 H6H.8 59.2 52 4 
34. 2 78.2 71.1 64.4 a8. 4 2 
6.2 72. ¢ 67.3 fi2 8.7 t 
8.2 67.8 ( l 62. 2 9.9 &. 0 
) 65.2 69 P 
42.2 63. 5 4.2 (4.7 65. 4 6. 6 
4.2 58.7 61.9 4. 1 tit 68.7 
4.2 63. § 57.9 62.6 66. 7 70.1 
418.2 47.4 3.8 0.9 66. 7 72.2 
KM). 2 $5.8 4.3 62.3 68.8 76. 2 
52. 2 2. 1 60.8 68.9 6. ¢ 84 ] 
b 55. 1 
Right 


* First image 
b Last image 
¢ Center setting 





better than +0.1 mm. The left column shows the 
successive positions of the telescope (or pentaprism) 
on the transverse bench. The words, Left and 
Right, or L and R, are included in the tables to show 
that the measurements are initiated at the left-hand 
side of the lens as viewed by the observer and are 
made at the designated intervals with the penta- 
prism being moved toward the right throughout the 
course of the experiment. The positions where light 
is first seen to enter the telescope is designated first 
image, and the position where it just disappears as 
the slide moves along the bench is designated las; 
image. No readings of the deviation € are possible 
at these positions, but these points serve to locate 
the center of coordinates, which in this instance 
corresponds to the 40.2-mm setting. Observations 
of the deviation, e, are made at 2-mm intervals along 
the transverse bench to a distance of 12 mm on each 
side of the central position. The values of the 
observed deviation are given in scale divisions for 
each of the five target settings. Each scale division 
is equivalent to 0.0283 milliradian deviation. ‘Table 
2 shows the values of A and e€ with respect to the 
center of coordinates; this table is derived from 
table 1. 


TABLE 2. Variation of parallax angle, €, as a function of zone 
height, h, for five positions of the target with respect to th 
¢ ‘ Je / ‘ : q / 
position of best focus of lens 2 
This table is derived from table 1. 
Nominal values of 6/, mm 
Zom 
height, h, 1.0 0.5 0.0 0.5 1.0 
mm 
Deviation ¢ from central position in scale divisions 
L 
14.9 
12 0.9 7.3 15.1 16. 1 23. 8 
10 17.0 8.7 0.0 SS 15.5 
s 1 10.7 $9 is 8 
f 13.0 6 l 4.6 14 
4 7.4 2.7 7 1 7.7 
2 2 6 0 7 3.1 4.2 
ti | { tv it iv 
2 1.7 4 1S 2.4 4.4 
t t 2.7 1.2 $ t 
6 11.7 6.7 0.3 3.7 7.9 
S 17.8 10.8 2.0 3.7 10.0 
10 14.4 10.3 0.6 Sg 14.0 
12 13. 1 3.5 6.0 13.6 21.9 
614.9 
R 


Last image 

The variation of ¢ as a function of A and Af for 
the conditions described in the preceding paragraph 
is shown graphically in figure 5. It is clear that 
curve 3 for Af=0 is a typical curve of angular spheri- 
cal aberration, whereas the other four curves show 
the contribution of varying amounts of out-of-focus 
effect. Values of Af for each value of A and « for 
each target setting listed in table 2 were computed 
with the aid of eq (3), and the results are shown in 
table 3. These results are shown graphically in 
figure 6, A. It is evident that the successive Af- 
versus-h curves are spaced along the axis of abscis- 
with separations closely approximating the 
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€,SCALE DIVISIONS 


Curves of parallax angle, ¢, versus zone height, h, 
for lens 2 showing the effect of successive target displacements 
of approximately 0.5 mm along the axis. 


Tp 


lhe results for a target position 1.0 mm farther from the lens than the plane 


best focus are shown in curve 1. Curve 5 shows the results for a target posi- 
1.0 mm nearer to the lens than the plane of best focus. Curves 2, 3, and 4 


the results for intermediate position (see tables 1 and 2). 


TABLE 3. ompute d values of Af as a function of zone he ight, 
h, for five positions of the target with respect to the position of 
hest focus for lens ? (see tables 1 and 2) 

e value of focal length, f, is 152.4mm. In converting « in scale divisions to e 
radians, | scale division =0.0283 milliradian. 
Ni nal values of 6f, mm 
Zone 
height, 1.0 0 0.0 0.5 1.0 
h, mn 
Computed values of Af, mm 
I 
12 0.05 0. 40 0.83 0. 88 1.30 
10 1.12 57 Oo 5S 1. 02 
s 1. 49 SS 32 31 0.81 
t 1, 42 71 16 50 1.03 
i 1. 22 14 12 71 127 
2 0). St lt 23 1. 02 1. 38 
0 
2 OM 13 ) 0. 74 1 5 
i 1. O07 +4 ‘*) 58 1. 07 
' 1. 28 ; 03 41 0. 87 
) 1. 4 ) 16 30 R2 
lf 1. 28 f O4 38 Y 
12 72 33 74 1. 20 
R 


iominal intervals separating successive positions 
M the target. The large negative values of Af are 
issoclated with negative values of the nominal 


he displacement with respect to the lens. 


gives the correct sign of 
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Figure 6. Af versus h curves for a series of target positions. 


A. Values of the displacement from best focus, Af, as a function of zone height, 
h, for the conditions shown in figure 5 (see table 3). 
B, Averaged values of Af versus A, derived from part A (see table 4). 


TABLE 4. Average values of Af as a function of zone height, 


h, for five positions of the target with respect to the position of 


best focus for lens 2. 


These results are derived from table 3. 


Nominal values of 6f, mm 


Zone | 
height, 1.0 —0.5 0.0 0.5 | 1.0 
h | in 
Average measured value of Af, mm 
mim 
2 0.71 0.14 0. 41 0. 90 1, 42 
i 1.14 44 .16 4 1.17 
6 1.35 72 —. 10 . 46 0.95 
S 1. 48 SS —, 24 . 30 . 82 
10 1. 20 —, §2 —, 02 . 48 07 
12 0.34 10 +, 58 81 1, 25 


The curves shown in figure 6, A are not completely 
symmetrical about the axis but can be made so by 
averaging the values of Af at corresponding values of 

h and —h. The resulting values of Af given in 
table 4 are shown graphically in figure 6, B. The 
curve for nominal 6f=0 is the normal longitudinal 
spherical aberration curve for the lens with respect 
to the plane of best focus for a zone height, A= 10 mm, 








b. Sources of Error 


When the averaged values of Af given in table 4 
are compared with the initial values shown in 
table 3, it is clear that the departures of the individ- 
ual values from the average value tend to vary some- 
what with A. This is shown in table 5, where the 
departures of the initial values of Af from the average 
are listed for each value of A and nominal 6f. The 
average departures for the five values of nominal 
5f are also given in the table. It is noteworthy 
that the average variation for zones 4, 6, 8, and 
10 are of the order of + 0.05 mm, whereas the varia- 
tion is +0.10 mm for zone 2 and + 0.21 for zone 12. 
A small error in the location of the zero position of 
h on the transverse bench scale will produce appre- 
ciable asymmetries in the initial values of « and 
Af; the values will tend to be too high on one side 
and too low to about the same degree on the other. 
However, this type of error does not affect the 
final average value to an appreciable extent 
Errors in Af arising from errors in / are proportional 
to the fractional error in A; thus an error of +0.1 mm 
produces a 5-percent error in Af for h=2 mm, but 
only a 1-percent error in Af forh=10mm._ Similarly 
errors in € are more serious for small values of h. 
For the case under consideration, an error in ¢ of 

0.1 scale division produces an error of + 0.03 mm in 
Af for h=2 mm and an error of +0.01 mm in A/ 
for h=6 mm. The large errors at h=12 mm are 
caused by asymmetrical use of the 2-mm viewing 
aperture at the edge of the lens 
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It is of interest to compare the values of Af ob- 
tained for each zone with each other and with the 
nominal vaues of 6f through which the target is dis 
placed. This comparison is given in table 6, where 
the difference in Af for each zone from the value for 
Af=0 1s given. It is clear that the successive values 
compare quite favorably with each other for zones 2 
through 10. It is, therefore, clear that quite reliable 
values of Af are obtained by this procedure, and il 
appears probable that the values of the longitudinal 
spherical aberration so obtained are reliable to 

0.04 mm for this particular lens 





TABLE 6 Computed values of the separations of th target 
. “ get 
from the position of best focus for lens 2 
rhese results are derived from table 4 
N nal value f I 
Z height. / Lf , 
Average measured value ¢ 4 
d in 0.49 1.0] 
1 60) M 1s , at 
1.2 2 iM ‘ L¢ 
‘ 1.24 4 ”) j 1 0 
1 18 ‘) ow ) 0.99 
] ow $s 67 


4.2. Longitudinal Chromatic Aberration of a Tele. 
scope Objective Having a Focal Length of 120 
Millimeters 


The axial longitudinal spherical aberration of 
lens usually varies with the wavelength of the image- 
forming light. Thus if the longitudinal spherical 
aberration Is measured for two different 
light, the characteristic curves of Af versus A will by 
displaced along the axis with respect: to each other 
The method of measurement described in section 4] 
provides a simple method of determining not only 
the longitudinal spherical aberration curve for 4 
given color but also the displacement zone by zon 
for differing illustrate this, results of 
measurement are given for a lens, designated 3. 
having a focal length of 120 mm and a clear apertur 
of 40 mm 

The same procedure was followed as outlined in 
section 4.1] except that three successive runs were 
made 
between the leht and the target In the 
instance Wratten monochromatic filters 74 
ereen), 75 (blue and 70 (red) were used 
and the results of measurement of parallax angle, 
height, A, are shown im 


colors of 


colors. To 


SOUuTCce 
present 


oreenD 


function of zone 
table 7 The variation of « with A is shown graphi- 
cally in figure 7, A for each of the three colors. The 
values of A} computed for each value of h are like- 
wise shown inh table : ana the averaged values of 
Af for paired values of h are shown graphically in 
figure 7, B. The position of the target was left un- 
changed during the entire experiment, consequently 
the Af-versus-A curves in figure 7, B, show both the 
longitudinal spherical aberration and the longitudinal 
chromatic aberration. It is noteworthy that the 
chromatic difference in between the light 
passed by filters 74 and 75 is negligibly small, where- 
as the difference for 74 and 70 is appreciable. Some 
out-of-focus effect is present, which shows that the 


é, as a 


focus 


target was approximately 0.27 mm farther from the 
lens than the position of best focus for filter 74. It 
can be inferred from the average values of Af, shown 
in table 7, that the chromatic difference in focal 
length between the light passed by filters 74 and 7) 
is 0.03 mm and between filters 74 and 70 is 0.10 mm 


For each run, a different filter was placed , 
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TaBLE 7. Varvation of parallax angle, e, and Af as a function 
of zone height, h, for three colors (74 green, 75 blue green, and 
70 red) 


The results are for lens 3, a telescope objective having a focal length of 120 mm 


Filter number 


Zome 74 75 70 
height, h ~ 
e, scale \/ e, Scak At e, scale Af 
liv mn divisions divisions 
mm mm mm mm 
2 4 0.31 12.5 0. 26 15.9 0.32 
It 11.9 30 9.7 25 9.0 23 
12 7.6 3) 6.6 22 44 —.15 
s 4.5 24 1.1 I 2.4 12 
4 26 2.3 23 1.4 14 
2 1.2 24 5 22 0.6 12 
) 0 0 0 
2 —L3 2t 13 2 5 10 
} 2 2t 2.5 26 —.9 09 
8 1. 25 1.9 2 2.4 12 
12 2t 7.1 24 1.8 16 
If 12.0 $1 11.0 28 &.7 22 
1) 16.4 ; 13.8 —. 28 13.1 27 
Average. ). 27 0). 24 0.17 


Objective Having a Focal Length of 
167 Millimeters 


4.3. Telescope 


In the course of making measurements on lenses 
of various types and focal lengths, one occasionally 
finds effects that are of unusual interest. The re- 
sults of measurement of longitudinal spherical 
aberration on this lens are here included because of 
the unusual form of the longitudinal spherical aber- 
ration curve. The lens, designated 4, is a telescope 
objective having a focal length of 167 mm and a 
clear aperture of approximately 46 mm. 

The test procedure was essentially the same as that 
described in section 4.2, except that only two filters, 
74 and 70, were used. The data for filter 75 are 
omitted, as the results obtained with it are nearly 
identical with those for filter 74. The results of 
measurement of parallax angle, e, as a function of 
zone height, A, are shown in table 8 and figure 8, A. 
The values of Af were computed for each value of h 
and are shown in table 8. The averaged values of 
Af for paired values of h are shown graphically in 
figure 8, B. The curves in figure 8, B show the 
longitudinal spherical aberration for each of two 
colors. It is noteworthy that the two curves are 
similar but are separated from each other by approx- 
imately 0.2 mm. 

It is interesting to note the shape of the Af versus 
h curves in figure 8, B. In the vicinity of the axis, 
the curves have a bulge toward the direction of 
shorter focal length. This anomaly does not occur 
in all lenses measured, but it does occur in enough 
to create a belief that this effect does not arise from 
experimental error. It is possible that it may result 
from a variation in the radius of curvature of one of 
the components of the lens under test. 
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Fiagure 7. Parallax angle e versus h, and Af versus h for three 
colors. 

A. The values of parallax angle e, as a function of zone height A, are given in 
curves 1, 2, and 3, for Wratten filters 74 (green), 75 (blue green), and 70 (red), 
respectively 

B. Values of longitudinal chromatic and spherical aberration Af derived from 
the measurements shown in part A (and table 7). 


TABLE 8. Variation of parallax angle, «, and Af as a function 
of zone height, h, for two colors (74 green and 70 red) 


rhe results are for lens 4, a telescope objective having a focal length of 166.7 mm, 


Filter number 


Zone ; 
height, 74 70 
h — aa wines 
¢, scale Aj ¢, scale Af 
divisions divisions 
mm mm mm 
23 12.2 —(). 42 10.1 —0. 34 
—21.5 8.1 —. 30 .8 —, 21 
20 5.¢ —. 22 2.7 —.11 
—15 2.8 15 —1.4 +.07 
—10 1.6 .13 —2.2 17 
7.5 1.3 .14 —1.7 18 
—5 1.5 —. 24 —0.9 .14 
—2.5 0.8 25 —.4 .13 
0 0 m S | «mnose 
2.5 7 —. 22 om 06 
5 1.3 —. 20 — ll | 
7.5 1.2 —. 13 1.5 16 
10 1.3 10 as | 2 
15 2.1 11 1.9 . 2 
20 —4.3 —.17 33 | —.09 | 
21.5 6.5 . 24 —4.8 | —.18 | 
23 9.4 32 —7.2 +. 25 
—EEE | 
Average —(), 22 +0. 01 
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Fiagure 8. Curves for a lens showing unusual characteristics 
of the longitudinal sphe rical abe rration oblaine d for iwo 
colors. 


curve 2 shows the 


in an unusual 


Curve 1 shows the results for Wratten No, 74 filter, and 
results for Wratten No. 70 filter. The curves appear to 
fashion away from the lens in the axial region (see table 8). 


5. Sources of Error 


In section 2, it is noted that an error in the location 
of the target reticle with respect to the focal point 
of the lens would produce an error in the measured 
values of the parallax angle, e, and that this error 
would be propagated into the computation of those 
values of Af from which the values of the longitudinal 
spherical aberration are derived. It is, however, a 
simple matter to so position the reticle that variations 
in the value of Af from this cause are less than 1 
percent and can be regarded as negligible. It 
worthwhile to consider the magnitude of the errors 
induced in Af by errors in the measured values of 
the zone height, 4, and the parallax angle, «. The 
effect of errors in h and ¢ are discussed briefly in 
section 4.1, b in connection with the measurements 
on a particular lens, but some additional comments 
are given in the following sections. 


Is 


5.1. Effect of Error in the Values of A 


There are two potential sources of error in the 
determination of the zone height, A. If the values 
of A are measured accurately but with respect to an 
incorrect zero position, the curve of ¢« versus fA is not 
symmetric about the origin and creates the appear- 
ance of large error. This is particularly noticeable 
in the derived curve of Af versus h, as it produces a 
very pronounced asymmetry. This error may be 
minimized or eliminated by averaging the values of 
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Af for corresponding values of +A and —h, JE x) 
usually, however, more aesthetically satisfying » 
make a more accurate determination of the correes 
zero position and to make the computations wij 


respect to this corrected position. I 
When the values of A are measured with respeg ! 
to the correct zero position, then the values of Af gy 


affected by a given fractional error in h by approx. 
mately the same fractional error. Thus for a ¢op. 
stant error in h, the fractional error increases wit! 
decreasing h and so places a lower limit on the yaly 
of A that may be profitably used. Assuming that} 
can be measured to the nearest +0.05 mm, th 
lowest value of A that can be used is + 1.0 mm if oy 
does not wish the error in Af from this source 
exceed 5 percent. For the same constant error in | 
the error in Af at h 10 mm 0.5 percent 
which may usually be regarded as negligible. 


Is 


5.2. Effect of Errors in « 


The magnitude of the observed parallax angle, , 
is usually quite small, so that errors in € provide th 
primary limitation in the determination of Af. Iti 
difficult to determine ¢ with an accuracy better tha 
+ 0.5 see because of limitations placed by pointing 
accuracy.2 It is probable that determinations of , 


can be made by using reasonable care, so that tly 
error of the mean does not exceed 1.0 sec. Fo . 
purposes of illustration, computations of the magn.’ * 
tude of the error in Af caused by an error of +2.0s . 
were made for a series of values of equivalent focal 
length, f, and the zone height, 4, and the results an 
shown in table 9, A. In these computations, th “ 
equation . 
a is . 
sie ‘ 
was used, with value of ¢ maintained constant at +)! : 
sec. Consideration of the values in part A of table A 
makes it clear that the magnitude of the error in 4 
increases steadily with f for a fixed value of h, and . 
inversely with A for a fixed value of f. For value 
of f above 100 mm, the errors in Af become appre- 
ciable for the selected value of the error in ¢«. Iti, 
therefore, interesting to consider these magnitudes of : 
the error in Af with the theoretical depth of focus, 1 
d,, of a lens in the image space. Approximate values 
of d, in millimeters are computed on the basis of the 
equation SI 
4h (| 
aa iy "6 
n 
where 6 is the f-number and /?,, is the maximum 0 


theoretical resolving power for light having a wave 
length of 575 mu.* The value of /?,, in lines per mill- 
meter can be computed with the aid of the formula 


, 1426 
hi, b 


?F. E. Washer, Effect of magnification on the precision of indoor teleseop 
pointing, J. Research N BS 39, 163 (1947) R P1820 
E. Washer, Region of usable imagery in airplane-camera lenses, J. Resea® 
NBS 34, 175 (1945) RP1636 
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TABLE 9. Effect of errors on the determination of longitudinal 
spherical aberration 

Comparison of the error produced in the determination of the longitudinal 

spherical aberration, 4/, by an error of +2 sec in the measured value of the paral- 

+ ungle,«, with the theoretical depth of focus in the image space, dy, that yields 


maximum theore tical resolving power, Ra, for light having a wavelength of 
575 Mg. 
Values are given in millimeters for a series of values of equivalent focal length 


f,and zone height, & 


4. Errors produced in the longitudinal spherical aberration, 


Zone \/, by an error of +2 seconds for various values of / 
height, - — = 
h 
100 200 400 800 
mm 
9 5 0. O10 +-() 039 +-() 155 +0). 62 2.49 
F OOS 019 O78 31 1. 24 
10 oo2 O10 039 16 0. 62 
”) OO5 019 O8 32 
40) .010 04 .16 
B. \ f theoretical total depth of focus, dy, yielding a 


resolving power F,, for the following values of / 
100 200 100 SOO 
2 (0). 280 1. 120 4.48 17. 92 71.7 
f 0) 280) 1.12 4.48 17.9 
10 O18 070 0. 28 1.12 4.5 
a 18 07 0). 28 L.1 
02 07 0.3 


It must be mentioned that values of d, obtained with 
eq (4) are for an ideal lens; for an actual lens having 
aberrations the value of d, may range from approxi- 
mately 2.86//?,, 1b/ Pn. 

To facilitate the computation of d,, the values of 
the f-number, 6, and the maximum theoretical re- 
solving power are shown in table 10 for each value 
of the zone height, h, and equivalent focal length, f, 
shown in table 9, A. The values of d, are computed 
with the aid of eq (4), and the resulting values are 
shown in part B of table 9 for ready comparison 
with the errors in Af shown in part A of the same 
table. It is readily seen that the values of errors in 
Af (or 6Af) are quite small in comparison with d, for 
small values of h. The relation can be shown to be 


to 


bAf 0.014hd,, (6) 
indicating that the relative magnitude of 6Af with re- 
spect to d, is independent of the focal length, f, but 
does increase with increasing 4. From these con- 
siderations, it is evident that although the error in 
Af is large at small values of h, these large errors are 
small in comparison with the usable depth of focus, 
d,. In addition, 6Af becomes small with large values 
of h, thus compensating for the increase in relative 
magnitude of 6Af with respect to d, for large values 


of h. 
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TABLE 10. Values of the f-number, b, and maximum theoretical 
values of the resolving power, Rp», used in computing the 
values of dy shown in part B of table 9 





A. Values of /-number, 5, corresponding to successive zones 
for the following values of f 
Zone . ae ae ee Te = 
height 
h 50 100 200 400 800 
mm 
2.5 10 } 20 40 80 160 
5 5 10 20 40 80 
10 2.5 5 10 20 40 
20 2.5 5 10 20 
40 2. 5 5 10 
B. Values of maximum theoretical resolving power, R», in 
lines per millimeter for \=575 my, for the following values of f 
50 100 200 400 800 
2.5 143 72 36 is | 9 
5 286 143 72 36 18 
10 572 286 143 72 36 
20 572 286 143 72 | 
40 a . 572 286 143 | 
6. Discussion 


This method of measuring longitudinal spherical 
aberration and longitudinal chromatic aberration has 
been used at the Bureau for some time. It has been 
used for lenses having focal lengths ranging from 
2 in. to 15 ft. and having apertures ranging from 0.5 
in. to 12 in. It has an accuracy comparable to the 
Hartmann photographic test method. It has the 
advantage of being a visual test method. It can, 
moreover, provide reasonably accurate values quickly 
and easily. 


The author expresses his appreciation to other 
members of the staff of the National Bureau of 
Standards for assistance during this work. Some of 
the measurements were made by Walter R. Darling. 


‘J, Hartmann, Z, Instrumentenk. 2, 1, 33, 97 (1904), 


WasHINGTON, February 26, 1958. 
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Mass Spectrometric Study of the Rate of Thermal Decompo- 
sition of Hydrazoic Acid 


J. L. Franklin,? Vernon H. Dibeler, and Preston P. Morris, Jr.’ 


Hydrazoie acid vapor at low pressures (0.15 mm Hg or less), has been decomposed 


thermally and the rate of reaction studied mass spectrometrically at 265° to 325° C. 


Nitro- 


gen and ammonia are the only reaction products, and no intermediates were observed. 
The disappearance of HN; is first order and the reaction occurs on the walls of the glass 


reaction vessel, 


The experimental activation energy is 31 


The 


kilocalories per mole. 


activation energy for the reaction HN;—NH-+- Ny, must be at least 39 kilocalories per mole. 


1. Introduction 


Recent studies [1,2] * of hydrazoic acid have shown 
that the heat of the reaction 
HN, — HN+N, (R1) 
is 6 to 9 keal/mole. However Franklin et al. [2] 
have pointed out, from considerations of spin-orbital 
correlations that this reaction cannot proceed to the 
>- ground state of NH and instead would have to 
vo to the lowest singlet (‘A) state of the radical. 
Thev estimated that this would be some 28 keal/mole 
above the ground state, so that the activation energy 
for the decompositions would be about 37 kcal/mole. 
For so weak a bond, one would expect the molecule 
to decompose readily and a measure of the activation 
energy for the thermal decomposition of hydrazoic 
acid should ordinarily correspond to the activation 
energy of reaction. However, Ramsperger [3] and 
Meyer and Schumacher [4] have shown that the de- 
composition of hydrazoic acid in glass occurs on the 
surface of the vessels and presumably is catalyzed by 
the glass wall. Consequently, in such a system the 
activation energy would at best set a lower limit to 
the activation energy for reaction 1. In spite of this 
limitation it was deemed worthwhile to study the 
kinetics of decomposition of HN, and this paper 
presents the results of this investigation. 


y- 


2. Apparatus and Procedure 


A gas analysis mass spectrometer, Consolidated 
Electrodynamics Corporation type 21-102, was used 
to observe the rate of decomposition of hydrazoic 
acid and the appearance of decomposition products. 
The all-glass sample introduction system and 2-liter 
sample reservoir was similar to that described by 
O'Neal and Wier [5] and used by them for obtaining 
mass spectra of high molecular weight hydrocarbons. 
The sample reservoir volume was enclosed in a well- 
insulated heating jacket and provided with thermo- 
couples at several points. Isolation from the sample 
troduction manifold and the evacuating pumps was 
accomplished with a gallium-sealed, magnetically 
,. This research was performed in part under the National Bureau of Standards 
free Radicals Research Program, supported by the Department of the Army. 

“iuest scientist on leave from the Humble Oil Co., Baytown, Tex 


Present address, Conv iir, San Diego, Calif 
‘Figures in brackets indicate the literature references at the end of this paper. 
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operated valve. The gallium valve, connecting 
tube, and the tube containing the effusive leak were 
wrapped with individually controlled heating tapes 
and with thermal lagging. Chromel-alumel thermo- 
couples were provided at the valve, the leak, and at 
several intermediate points for temperature meas- 
urement. The maximum operating temperature 
was about 350° C. At 300° C, temperatures could 
be read and controlled with an uncertainty of about 
+4° C. 

Conventional mass spectrometric techniques were 
used for both chemical analysis and observation of 
the kinetics of decomposition. The sample reservoir 
also served as a reaction vessel. Prior to each run, 
it was thoroughly outgassed, usually overnight, at 
temperatures somewhat above the anticipated tem- 
perature range. At the start of a run and after a 
stable operating temperature was attained, neon was 
introduced into the reservoir to a nominal pressure 
of about 0.05 mm of mercury as measured by gas 
pipet. The mass range of 16 to 46 was scanned 
quickly to serve as a background. Then hydrazoic 
acid was added to a partial pressure of about 0.1 mm 
of mercury and the time of introduction was noted. 
In order to make initial observations as frequently 
as possible, only selected portions of the mass scale 
were scanned : e. g., m/e=16 to 20, 28 to 32, and 40 to 
46. Each range uniformly required about 1 min to 
scan and the total lapsed time was noted at the start 
of each scan. Measurements were usually continued 
until the molecule ion of HN; (m/e=43) had de- 
creased to about 1 percent of the original abundance. 
The electron energy was 70 v. 

Hydrazoic acid was prepared from sodium azide 
and 60 to 80 percent valle acid. Preparations were 
carried out in a vacuum and the HN, was purified 
by simple bulb-to-bulb distillation and by pumping 
on the product held at liquid nitrogen temperature. 
A mass spectrometric analysis of the product showed 
no impurities of higher molecular weight or of most 
simple gases condensible at liquid nitrogen tempera- 
ture, the exception being occasional small amounts of 
ammonia. For reasons of safety the quantity of 
stored hydrazoic acid was limited to that giving a 
pressure less than 50 mm of mercury at room tem- 
perature in a 2-liter bulb. As HN, was observed to 
decompose slowly in storage under these conditions, 
periodic analyses were necessary to assure material 
of suitable purity. 








3. Results 


Table 1 summarizes the experimental data. 
Hydrazoic acid and neon concentrations after various 
reaction times are given in all runs. The nitrogen 
and ammonia concentrations are given for purposes 
of illustration only for run 6. The only products of 
the decomposition of hydrazoic acid that were ob- 
served were ammonia and nitrogen. Although both 
hydrogen and hydrazine might have been expected 
none was ever found. Data from a typical experi- 
ment are plotted in figure 1. Here peak heights are 
corrected for differences in sensitivity and the mass 
28 peak is corrected for the contribution from HNs. 
The curves for ammonia and nitrogen are character- 
istic of those for primary products of reaction formed 
simultaneously from the initial reactant. 

It is shown later that the decomposition of hydra- 
zoic acid is kinetically first order. For this to occur 
with ammonia as the final product, at least one 
transient intermediate required. This might 
include NH or electronically excited HN;, as the 
first active species and various secondary intermedi- 
ates such as N,H», NH», ete., can be postulated. 
Most mechanisms for the decomposition of HN 
have postulated NH as the initial product and we 
thought it worthwhile to try to establish this by 
identifying the products of its reaction with other 
materials. Specifically we made experiments in 
which large excesses of benzene and toluene were 
added to the reaction mixture and the products 
analyzed for postulated products of reaction. Thus 
if NH were present, one might expect either biphenyl 
or aniline from the reaction in the presence of ben- 
zene, but none was identified, although the reaction in 
solution gives aniline. Similarly, in the presence of 
toluene one would expect bibenzyl, benzyl amine or 
dibenzyl amine to be produced but none was ever 
found. Rice and Freamo [6] report that cyanic acid 
was formed when CO was added to HN, at high 
temperatures but the identity of masses of HNCO 
and HN, rendered this method of identification un- 
suitable to our studies which depended upon mass 
spectra for identification and analysis of the reaction 
products. Similar considerations prevented our 
using ethylene which might have been expected to 
add NH giving ethylene imine. 

Each mass spectrum was also examined carefully 
for evidence of transient intermediates but none was 
ever found. Of course NH would be obscured by 
the NH* peak from HN,. However, NoH, and N.H, 
if present would have been strong evidence for the 
presence of NH and NHb, respectively. In two or 
three runs small peaks at m/e=30 and 32 were found. 
However they always occurred together and 
peak at m/é 31 was ever observed. Therefore i: 
was concluded that the m/e=32 peak was due to O, 
introduced as an impurity and the m/e=30 peak was 
NO formed by reaction O, with either N. or HNsg. 
Thus it was impossible to obtain any direct evidence 
as to the nature of the active intermediates involved 
in the reaction. 

Although the reaction appeared to be kinetically 
first order, it might conceivably have involved a 
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FIGURE 1. Variation of concentrations with time for run 6. 


self-regenerating chain process. Accordingly a few 
percent of NO was added in one of the runs to ascer- 


tain whether the rate of reaction was markedly 
reduced thereby. It is evident from table 1 that 
NO caused only a small reduction in rate. This 


reduction may be enough to indicate that free radi- 
cals are involved in the reaction and are partly 
removed by the NO. This very small reduction in 
rate, however, clearly implies that a chain reaction 
is not being interrupted. 


- 4. Kinetics 


The gases in the reaction vessel passed into the 
mass spectrometei by molecular flow, the rate for 
various components of the gas mixture being inversely 
proportional to the square root of the mass. In all 
runs except one, neon was mixed with the hydrazoie 
acid in the reactor and the rate of disappearance of 
the neon, as indicated by the peak at m/e=20, was 
employed to calibrate the rate of flow of other com- 
ponents into the mass spectrometer. The rate ol 
flow of neon from the reactor should be proportional 
to the neon concentration so that 


[Nev] 


n [Ne] ket 


| 


where k, is the rate constant in reciprocal seconds. 
Data from a typical run plotted im figure 2 bear this 
out. 
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TABLE 1. 


(Peak heights (PH) 


Reactor packed with 


Summary of experimental data 


are given in arbitrary units.) 


3 4 
300 275 
0. 047 0. O85 
0. 150 0. 060 
HIN Ni HIN Ne 
PH t PH t PH t 
mi min min 
1740 0 2700 
1140 2 2085 ; 1954 l 
45 17 1800 19 1593 17 
iH5 1 1545 4 1290 52 
189 47 1320 14 493 47 
a) (2 1110 6 729 62 
1! 7 W440) 79 531 77 
21 v2 790 WY ISO} 92 
s 113 620 108 21 106 
143 430 124 138 121 
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Hydrazoic acid will disappear from the reaction 
vessel by reaction and by flow into the ionizatiorr 
chamber. The rate of disappearance is 


d{HN3] 


Fr k.[HNs], 2 


k{HN3]’ 


are respectively, the rate constant 
for reaction and for flow out of the reactor and 7 is 
the order of the reaction. From figure 2 it appears 
that [HN,] is an exponential function of time and 
that the reaction is first order i. e.,n=1. However, 
since the over-all disappearance is influenced by the 
flow rate, which is first order, it seemed desirable to 
apply a more sensitive test of reaction order. For 
this the differential rate equation is more suitable 
than the integrated expression. Equation (2) can 
be rewritten 


where Rp and k. 


dt 
By calculating the slopes of tangents to the curve 
of [HN,] versus time, such as that in figure 1, 
din{HN,}/dt is readily calculated for various points 
along the curve. If this parameter varies with the 
concentration of HN, the order can be obtained 
graphically since k, is known. If dln{HN,]/dt is 
independent of the HN; concentration the reaction is 
first order. A typical set of results is given in table 
2 and shows that within experimental accuracy dln 
[HN] is independent of the HN, concentration and 
thus the reaction is indeed first order. Since n=1, 
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(2) becomes —d[HN3]/dt=(k,+-k.)[HNs] and (k, +k) 
can be determined from the slope of a plot of 
InfHN,] against ¢. Also k,=kyw M, M,. where M is 
molecular weight, so k, is determined. The rate 
constants for all runs were determined in this way 
and the results are given in table 1 . 

Ammonia and nitrogen are formed simultaneously 
by reaction and are withdrawn by effusion into the 


mass spectrometer. Since the total reaction jg 
3HN,—-NH,;+4N, we may write 
cUN al _ (HN,|—ks [NH ‘ 
cMal_* (HN,|—k, (NJ 


where k; and ky are, respectively, the rate constants 
for effusion of ammonia and nitrogen from the reae- 
tor. Integration leads to the following relations 
between concentrations and time: 


(HN;] 5. aa i 
[NH,} 3(k Lke—k exp kt) —exp k,—ke)t 
(NH;] exp kst) (6 
», 1.88% TEN, 
INI="y “ ‘ exp (—k,t)—exp (—k,—k,)t) 
[N.] exp kt). (7 
Figure 1 gives typical curves of |NHj,] and [N 


against time obtained experimentally and by calcula. 
tion from eq (6) and (7). The agreement leaves 
something to be desired but is close enough to be 
considered as confirmation of the method. 

diflerentia ale equation 


TABLE 2 Test of reaction orde hy 


lor run oo 


f HN jJHN 1000d In HN 
\s 
m +, 890 70 18. 43 
200) 4 18.73 
4) vn) 2 »”).2 
) 1 10.0 18.6 
t 1, S00 2 18.1 
ri 1, 500 26.3 I 
si 27 1.3 
| 1, OF 18.7 17.4 
loo 400) 15.0 lf 
11 75 13.8 x4 


5. Activation Energy 


An Arrhenius-type plot of reaction rate constant 
versus temperature is given in figure 3. 
although somewhat define a reasonably 
good straight line from which the activation energs 
Eft, is found to be 31 keal/mole. The rate constant 


is thus 


scattered, 


k---4.5« 108 exp (—3 1000 RT) see. is 
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FicgurRE 3. Arrhenius plot of activation energy. 
Rice and Freamo [7] decomposed HN, thermally at 
1,000° C and allowed the gases to impinge on a cold 
surface, thereby obtaining a characteristic blue 
solid. By varying the distance from the furnace 
to the cold surface they obtained a half-life of about 
10-* seconds for the active species which formed the 
blue solid and which they concluded was NH. 
Their results might also be explained as giving the 
half-life HN, which decomposes to an active 
intermediate (such as NH) of very short life. If 
we calculate the half-life of HN, at 1,000° C from 
eq (8) we obtain a value of 3 107* see which, in 
view of the fact that the two sets of measurements 
were made under very widely differing conditions, 
is a reasonable check on the results of Rice and 
Freamo [7]. 

Glasstone, Laidler, and Eyring [8] show that the 
rate constant may be expressed in terms of the 
experimental activation energy 


of 


_ _ Et 
oh RT) 


St 
exp (2+ PR (9) 
‘ 


For reasons that will be given later (see eq 10, 11, 
and 12) the right-hand side of eq (9) must be multi- 
plied by 3 to give the over-all rate constant. Taking 
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Et=31 kcal/mole and employing eq (8) and (9) we 
calculate AST to be —26 e. u. This is much too 
large a negative entropy for a gas phase, first order 
reaction. The value is best explained as resulting 
from a reaction occurring on a solid surface. Indeed, 
such an entropy value corresponds approximately 
to the entropy of condensation of a gas. 

In order to corroborate this conclusion a run was 
made using a reactor packed loosely with glass wool. 
The results are not as satisfactory as those obtained 
with an unpacked reactor, probably because diffusion 
was impeded by the packing. However, they do 
show that the added surface greatly increased the 
reaction rate. We may conclude, therefore, that 
the reaction occurs predominately on the glass 
surface of the reaction vessel. 

Since the reaction is kinetically first order and the 
surface is involved, two kinds of mechanism can be 
adduced which would be in accord with our observa- 
tions. 

1. HN; is adsorbed reversibly and the surface is 
sparsely covered. The adsorbed material undergoes 
a series of reactions resulting eventually in the pro- 
duction of nitrogen and ammonia which are Preset st 
The following mechanism is one of several possibilities 
of this kind 


’ 3 f ‘5 
HN;(g) s HN, (ads) 
b 
HN;(ads) — HN (ads) +N,(g) \ 
(R2) 


‘ 
HN (ads) + HN;(ads)—>N,H, (ads) + N,(g) 





d 
N.H.(ads) +HN,(ads) — NH;(g) +2N,(g) J 


If one assumes that reaction b is the slow, rate- 
controlling step and employs a steady-state treat- 
ment for all transient species the following rate 
expression is obtained: 


—d{HN;(g)]_ 


7 = 3k, Kua [HN;(g)] 


(10) 


where A is the equilibrium constant for the adsorp- 
tion of hydrazoic acid on glass. 

2. HN, reacts in the gas phase to give an active 
species which adsorbs rapidly and then undergoes a 
series of reactions with gaseous HN;. It is assumed 
that the surface is almost completely covered by the 
active species. As an example 


HN;(g) — NH(g)+N,(g) (Rl) 
HN(g) <> NH (ads) 
HN (ads) +HN;(g) — N,H,(ads) + N2(g) (R3) 


N.H,(ads) +HN;(g) — NH;(g)+N,(g) 


Again, if one assumes the first reaction to be the 
slow, rate controlling step and applies a steady-state 
treatment, the reaction rate is found to be 








= 


—d{HNs( 


g)]__ a, lle 
7 3k, [HN;(g)] 


(11) 


Both kinds of mechanisms give first order kinetics 
for the decomposition of HN,;, but the over-all rate 
constant for mechanisms of type 1 includes the 
adsorption equilibrium constant whereas those of 
type 2 do not. 

Mechanisms of type 2 usually prevail for the re- 
combination of atoms or simple radicals where a 
third body is required to remove the energy of reac- 
tion. It seems unlikely that this should be the case 
here since various reactions of the probable interme- 
diates other than simple recombination can readily 
be visualized and these would provide a means of 
removing excess energy without the need of a third 
body or wall. Further, the rate of disappearance of 
HN, by this mechanism is given simply by 3/; and so 
is independent of surface effects. This would not 
account for the negative entropy of activation and 
so mechanisms of this kind can be discarded. 

Since the rate constant for mechanism 1 includes 
the equilibrium constant for adsorption of HN, on 
glass, we may write 


' ap w.okL (4S? AH* AHgas , ASuss 

ker = Sho K=—, exp ( R-RT~ RTT R) 
3kT —— (AS*+- AS ya —AH'— AH 
exp (Re) exp (Re ) 


(12 


ASt for a unimolecular reaction is usually small 
and, as a first approximation may be taken as zero. 
The entropy of adsorption is negative and is often 
comparable in magnitude to that obtained for AS' 
in these studies. 

In deriving eq (10) and (12) we assumed that the 
adsorption of HN, was reversible and that equilib- 
rium existed at all times between HN, in the gas 
and on the surface. At equilibrium, AF=O and so 

AH TASyas 


ids 


The heat of adsorption then, at temperatures near 


300° C would be about 8 keal/mole. Since 
Fi. = 4H + Abas 
AH; =39 keal/mole. 
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The activation energy, AH}, is, of course, the ac. 
tivation energy of reaction b, the slow Step in the 
In the example given this is shown 8 4 
bond breaking reaction but it might also be a simple 
elevation of the molecule to a higher electronic state 
Since reaction b occurs on the surface, AH? is po 
necessarily the same as it would be for the gas phage 
reaction. Further, it is probably somewhat smaller 
that it would be for the gas phase process, gine, 
otherwise, in view of the small entropy of activatio, 
of a first order reaction in the gas phase, the reaction 
would occur in the gas and not on the surface, Algo 
if the rate controlling step is some reaction other 
than the simple bond breaking process shown in the 
example then it follows that the activation energy 
for the gas phase reaction 


process. 


HN, — NH-+N, (Ri 
is greater than AH? since otherwise the latter rege. 
tion would occur at a greater rate and would be the 
process observed. It follows therefore that the meas. 
ured AH? is equal to or less than the heat of activa- 
tion of reaction R1; i. e., AH* for R1 is at Jeast 39 
keal/mole. This is surprisingly close to the value of 
37 keal/mole estimated by Franklin et al.2 It would 
appear from these data that the true value oi the 
activation energy of R1 is somewhat greater than 
this, but the amount would be difficult to estimate 
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Ta High-Speed 
ae Erosion Damage to Solids Caused by High-Spe 
Simple 
7 . *. . . * 
C State 
c stat Collision With Rain 
S hot 
S Phase Olive G. Engel 
Smaller 
5, Since The stresses that are produced when a liquid drop coliides with the planar surface of a 
iVatioy solid are discussed. These are a result of the impact pressure and the consequent radial 
eacti : flow of the drop. It is concluded that a rain-erosion resistant material may be either soft 
— ind rubbery or hard and rigid. If it is soft and rubbery it mitigates the applied stresses but 
Also it must be able to withstand the mitigated stresses. If it is hard and rigid it does not miti- 
n Other gate the applied stresses and it must be able to withstand the unmitigated stresses 
1 1n the , ie : . 
energy 1. Introduction feet per second (fps). The diameter of this damage 
mark, as measured directly on the aluminum plate 
The problem of the erosion damage that is pro- with a steel rule and magnifying glass, is about %e in., 
R} duced on objects that fly at high speed into rain has which is considerably less than the diameter of the 
become increasingly more serious as flight velocities steel sphere. Figure 1, B shows an example of a 
er reac. have been increased. damage mark that was made by a 0.22-in.-diam de- 
L be th During the past 10 vears, a program of testing forming lead pellet on \-1n. 1100 H14 aluminum 
€ meas- has been carried out to find materials that will resist Sheet. This particular damage mark was obtained 
activa. | the stresses that are brought to bear in high-speed ® 4 relative impingement velocity of 609 fps. 
least 39 collision with rain. The work has been done mainly 
value of at the Cornell Aeronautical Laboratory. Speci- 
t would mens of materials have been fastened to the ends of 
> Or the a propeller and have been rotated at velocities up 
er thar to approximately 600 miles per hour (mph) through 
‘timate artificial rain. Tests of a different kind have been 


’ carried out at higher velasialal by the Convair Divi- 
sion of General Dynamics Corporation. 
During the past 5 years, research has been in 
[A] 235, progress at the National Bureau of Standards to 
determine the mechanism of the erosion process. 





and M . . » 

™ A previous publication ' has presented three fac- 

+ (1929 tors that are important in analyzing the damage " 

hem. {A that is produced as a result of high-speed wate rdrop a lelmm 


collision with poly(methyl methacrylate): (a) The 
waterdrop behaves like a hard sphere or pellet in 
75, 548 that it produces a cup-shaped cavity in the plastic; 
(b) although materials do not usually fail under uni- 


hem. 23, ; 


73, 552 form compression, a local compression can produce 
e Theor tensile stresses that result in failure; and (c) in the 
Co. th case of projectiles that flow as a result of the colli- 


sion, the solid plastic can be damaged by the stresses 
exerted by the flowing projectile against surface ir- 
regularities that are restrained by the underlying 
layers of material. 


2. Collisions of Waterdrops with Aluminum 


Further evidence of the behavior of water under 
impact conditions was obtained in a study of the 
damage marks left by deforming lead pellets, steel 
spheres, and waterdrops as a result of high speed 
collision with soft 1100 aluminum sheet. — 
1, A shows an example of a damage mark made by : 

-in.~iam steel sphere on -in. 1100-H14 shenab- 





num sheet. This particular damage mark was ob- 
tained at a relative impingement velocity of 228 c — Imm ke 
Figure 1. Damage pits in 1100 aluminum formed as a result 
i The work dk scribed in this paper os sponsored by the Weta Air Develop- of collisions with (a) a steel sphere at a velocity of 228 fps, 
men enter le material presented here has been taken largely from reports 
submitted to them. Their assistance and interest are gratefully acknow —— (b) a lead pellet ata velocily of 609 fps, and (c) a wate rdrop 
' Olive Engel, J. Research NBS 54, 51 (1955) RP2565 at a velocity of 2,500 fps. 
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At a collision velocity of 880 fps (600 mph), 
waterdrop damage marks that have been observed 
in 1100—O aluminum resemble the mark left in 
1100-H14 aluminum by impingement of a de- 
forming lead pellet at very low velocity. When a 
waterdrop collides with 1100 aluminum at about 
2,000 fps, an impression is made which resembles 
that made by a low-velocity steel sphere more nearly 
than that made by a relatively low-velocity lead 
pellet (see fig. 1). These observations lead to the 
conclusion that in the range of collision velocity up 
to about 1,000 fps a waterdrop behaves like a pro- 
jectile that flows as a result of the collision, but at 
collision velocities of the order of 2,000 fps a water- 
drop behaves very much like a projectile that does 
not flow as a result of the collision. 


3. Destructive Action of the Radial Flow of a 
Waterdrop 


When a waterdrop collides with the planar surface 
of a solid, or when the planar surface of a solid runs 
into a stationary waterdrop, the impact pressure that 
results reaches a high value in a very short time.’ 
This high pressure drives the liquid that is close to 
the solid surface radially outward around a central 
stagnation point. The radially flowing liquid exerts 
a shear stress on the surface of the solid over which 
it is running. There is a shear stress between the 
separate layers of the flowing liquid itself, and it 
can be expected that a shear stress exists at the inter- 


face between the liquid and a solid surface over, 


which it is moving. The shear stress, r, between 
layers of liquid in laminar flow is given by the prod- 
uct of the viscosity, uw, and the velocity gradient 
through the moving sheet of liquid perpendicular to 
its direction of flow. That is, r—u(Or/Oz) where v 
is the velocity of the moving sheet of liquid and 
is the direction through the thickness of the liquid 
sheet. The layer of liquid molecules in direct con- 
tact with the solid has zero velocity, but the velocity 
gradient is not zero and the shear stress is applied to 
the solid. 

When a surface protrusion exists, the liquid from 
a drop that is flowing radially over it exerts forces 
against the protrusion acting in Opposition to the 
cohesion of the protrusion to the underlying layers of 
material. Pressure, o, of figure 2, is exerted against 
the protrusion by the flowing liquid. The pressure 
that is exerted by the liquid tends to move the pro- 
trusion along the planar surface of the solid and re- 
sults in a shear stress, 7’ of figure 2, at the base of 
the protrusion. The flow of the liquid results in the 
shear stress, 7, that was discussed above. The pres- 
sure exerted by the liquid also results in a turning 
moment that tends to bend the protrusion over 
The turning moment is the integrated cross product 
of the compressive force exerted by the liquid and 
the elevation of the protrusion above the planar sur- 
face of the solid at the point where the force is 
applied. As the protrusion bends, a tensile stress, 
o, of figure 2, appears on the side of the protrusion 


2 Olive Eng J. Research N BS 54, 281 (1955) RP2591 
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flow of an impinging wate rdrop 


against which the compressive force is applied and q 
compressive stress, oj of figure 2, appears on the 
opposite side of the protrusion. If the force exerted 
by the rapidly running liquid is sufficiently great, or 
if the protrusion has a sufficient elevation above the 
planar surface, failure may occur. Whether the pro- 
trusion is simply bent over, whether it is broken off. 
or whether part of the solid material below the sur- 
face is torn out with it depends on the strength of the 
material and on where failure occurs first. Likely 
points of failure are marked with notches in figure 2, 

If a waterdrop impinges against the planar surface 
of a rubberlike material the effect of the compressive 
stress, o, of figure 3, which is the direct result of the 
collision, is to compress the solid material at the point 
of impingement. The local compression results in 
a dimple or cup-shaped cavity. Tensile stresses, ¢, 
of figure 3, appear in the sides of the dimple and 
especially in the knee in the solid material at the rim 
of the dimple. If these tensile stresses are sufficiently 
great, tears may form in the surface layers of the 
solid material. The material that formerly occupied 
the volume that now forms the hollow of the dimple 
is displaced. It is moved radially outward and up- 
ward around the dimple and this results in a shear 
’ of figure 3. The radial flow of the liquid of 
the drop imposes a radial tensile stress, o; of figure 3, 
and a shear stress, 7, which was discussed above. In 
the case of very thin rubbery coatings it is possible 
that a portion of the shear stress, +r, may be trans- 
mitted through the coating to the primer bonds, 


stress, 7 
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which hold the coating to the base metal, and that it 
may cause failure of the primer-to-rubber or of the 
primer-to-metal bond. However, the removal of a 
circle of rubber coating around the central point of 
the collision (stagnation point of the flow), as oc- 
curred in the case of a Neoprene coating (see fig. 6), 
i; not necessarily the result of transmission of this 
shear stress through the coating. It is possible that 
a deep circular cut may form as a result of the tensile 
stress, , of figure 3, in the material at the rim of the 
dimple or cup-shaped cavity. If this cut extends 
through the coating the radially flowing liquid may 
lift the coating with a wedge action and peel it off. 


3.1. Waterdrop Traces in Fine Moist Sand 


The effect of the forces that come into play when 
the liquid of an impinging drop flows radially can be 
seen in the marks that are left when waterdrops fall 
into fine moist sand. To obtain the picture of such 
traces shown in figure 4, very fine mesh sand that 
was just moist was gently pressed into a large petri 
dish. The surface of the sand was leveled and was 


then lightly covered with carbon black to as uniform 
a degree as possible. Waterdrops were allowed to 
fall from a flat-nosed pipet through a distance of 
approximately 40 ft to the surface of the carbon- 
covered sand. The waterdrop traces in figure 4 
show clearly the existence of the stagnation point at 
the center of the flow where the radial flow velocity 
is zero; at this point the carbon has not been washed 
from the surface of the sand. The circular trench 
about the stagnation point is the result of the forces 
that act as a consequence of the radial flow of the 
liquid. The relative collision velocity was about 
25 fps. 


3.2. Waterdrop Traces Produced in a Neoprene 
Coating 


Traces similar to those that formed when water- 
drops were allowed to fall into the petri dish of 
carbon-covered sand were also found on a Neoprene * 
coated specimen that was rotated on the Cornell 


This was not one of the standard Neoprene coatings that fulfilled the require 
ments of Specification MIL-—C-7439. 
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CO lide d u ith 


sand that was dusted with carbon black. 








Aeronautical Laboratory rain erosion test apparatus 
at a velocity of 600 mph through 1l-inch per hour 
(iph) artificial rain for 1.5 min (see fig. 5). The 
average coating thickness was 8 to 10 mils. Because 
the surface is not completely covered with such 
damage marks it seems likely that traces of this kind 
are not produced by the first waterdrop collision 
with the unfatigued Neoprene coating. They are 
probably formed during a single collision but only 
after the coating has suffered deterioration of some 
kind (possibly permanent set) as a result of numerous 
radial-stretch-and-recover cycles imposed by the 


radial flow of waterdrops that collided earlier but 
produce 


that were not themselves able to the 


observable damage. 





—| Imm — 
Fiaure 5. Circular damage marks produced on an experi- 
mental Ne oprene coating that collided with waterdrops at a 


velocity of 880 Ips. 
3.3. Waterdrop Traces in Polymeric Materials 


Figure 6 shows damage marks that were produced 
on four different polymeric materials by collision of 
the specimen material with waterdrops at a velocity 
of 2,250 fps. The experiment was carried out at 
Convair. The Neoprene coating was one of the 
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standard Neoprene coatings that fulfilled the require. 
ments of specification MIL—C-7439. Arrows jn the 
picture point to the circular w aterdrop traces. The 
undamaged centers, which mark the stagnation point 
of the radial flow. are very clear in the case of the 
two traces at the left of the picture. Circular dam. 
age traces were also produced in Teflon, Selectron 
5003 polyester resin, and methyl methacrylate 
plastic. These pictures show that the forces exerted 
by the radial flow of a waterdrop as a result of 
collision with the planar surface of a solid at 
velocity of 2,250 fps are effective damage tools, 


4. Response of Structural Materials to High. 
Speed Waterdrop Impingement 


The damaging properties of an impinging liquid 
drop that have been considered in the preceding 
sections are similar for drops of all liquids but vary 
in intensity depending on the density of the liquid, 
on the relative collision velocity, and on the extent 
to which the solid material vields under the blow. 
If the characteristic properties of all structural mate- 
rials were the same, the appearance of damage marks 
produced DN the impingement of drops of a riven 
liquid against all of them at some arbitrary velocity 
would also be the same: there would be only one 
type of damage and only one mechanism by which 
the damage is produced. Because the character- 
istic properties of structural materials are different. 
there are as many different damage processes as 
there are broad groups of material properties, 
A specific question that was posed at the time that 
the research program on the mechanism of high- 
speed rain erosion was initiated was whether a soft 
rubbery material or a hard rigid material should be 
sought as a solution to the rain-erosion problem. 
These two extreme are in the 
following sections. 


ol 


causes considered 


4.1. Rubberlike Yield 

A material that behaves like rubber under the 
blow that results from collision with a liquid drop 
has the advantage of being a pressure reducer for 
the blow. When the impact pressure is reduced, 
the velocity of the radial flow of the liquid of the 
drop, which is driven by the impact pressure, is also 
reduced. Because the radial flow velocity is reduced, 
the shear stress that it exerts, which decreases as 
the radial flow velocity decreases, is likewise reduced. 
The material is, however, depressed at the point 
where the liquid drop. struck and radial tensile 
stresses appear in the sides of the cup-shaped clepres- 
sion and around and especially over the rim of it 
it forms. The surface layers of the material 
around the restricted area of contact between the 
drop and the surface of the solid are also stretched 
as a result of the rapid flow of liquid from the drop. 
A material of this kind need only have strength 
properties sufficient withstand the mitigated 
stresses in order to remain undamaged on colliding 
with a liquid drop. 
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TEFLON SELECTRON 5003 


NEOPRENE METHYL METHACRYLATE 


Circular damage marks that were produced by collision with waterdrops at a ve locity of 2,250 fps. 
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Increase in the indentation hardness of a rubber is 
commonly accompanied by a decrease in its elonga- 
tion under a given stress and, therefore, by a decrease 
in its ability to mitigate the stresses that are imposed 
on it by high-speed collisions with a waterdrop. The 
indentation hardness associated with maximum re- 
sistance against waterdrop-impingement damage at 
an arbitrary collision velocity is that hardness above 
which the gain in strength that often accompanies 
further hardening is unable to offset the loss in 
ability to mitigate the stresses that are imposed by 
collision with a waterdrop at the velocity in question. 
In terms of the rupture energy per unit volume, the 
optimum indentation hardness is that hardness 
above which the area under the stress-strain curve 
begins to decrease. 


4.2. Resistance of Hard Materials 


All materials yield to some extent under a collision 
blow. However, materials that do not display a high 
elongation under a given stress as a result of collision 
with a liquid drop do not mitigate to a notable degree 
the stresses that the colliding drop exerts. To be 
erosion resistant, materials of this class must be able 
to withstand the unmitigated stresses. Whether or 
not they can withstand the unmitigated stresses de- 
pends on whether they can absorb the energy of the 
collision before their yield or fracture strength iS eXx- 
ceeded. If the relative collision velocity between the 
solid material and the liquid drop is increased without 
limit all the known hard materials will be found to fail 


at the point where their yield or fracture strength is~ 


exceeded ; they will either yield with plastic flow or 
they will shatter. 

Poly(methyl methacrylate) is a practical, rigid- 
material solution to the rain-erosion problem for flight 
velocities up to 300 mph. When specimens of 
poly(methyl methacrylate) are rotated through 1-iph 
simulated rain at a velocity of 300 mph, initial pitting 
of the plastic occurs after 45 min. If, however, the 
velocity is increased to 400 mph, initial pitting occurs 
after 2 min. Test results show that polystyrene and 
the glass-cloth-reinforced polyester resins are less 
resistant to rain erosion than poly(methyl metha- 
crylate). The combined properties of rigidity and 
low strength eliminate these polymers as practical 
rain-erosion resistant materials at ordinary airplane 
flight velocities. 

Metals undergo permanent plastic flow when their 
yield strength is exceeded. In the case of 1100-O 
aluminum, the depression in the surface caused by a 
single collision with a waterdrop at a velocity of 600 
mph is barely perceptible. The pummeling action of 
repeated collisions produces, by plastic flow, a general 
unevenness of the surface. This unevenness appears 
to play a role in enhancing a delayed, drastic failure 
of the material. Moderately thin specimen plates 
of the soft and medium hard metals undergo con- 
siderable plastic flow during a single collision with a 
waterdrop at velocities in the range of about 1,000 to 
2,000 fps; the characteristic damage mark is a pit 
that resembles an empty spherical segment. See 
figure 1, C 
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The hard metals are known to afford a practicg} 
solution to the high speed rain-erosion problem aj 
moderately low impingement velocities. Incipient 
failure shows up only after hours of waterdrop im. 
pingement in the erosion test at a velocity of 500 mph 
in a 1-iph rain density. The failure of the hard metals 
in this range of impingement velocity cannot be ex. 
plained in terms of the impact pressure produced by 
collision with a waterdrop. It appears to involve 
pressure-raiser defects, such as microscopic pits and 
pores on the surface of the specimen, or weak spots jn 
the material at or close to the surface, or both. 

Of the materials that have been tested for resist. 
ance against collision with waterdrops at high relative 
impingement velocities, the ceramics have been found 
to be among the most resistant. However, if the 
impingement velocity is increased without limit, these 
materials also shatter at the point where their frac. 
ture strength is exceeded. A practical tensile 
strength of 30,000 psi can be expected from poly- 
crystalline alumina compared to lower values for the 
soft metals and plastics. 


5. Requirements for Rain-Erosion Resistance 


The question that was posed several years ago as 
to whether the solution to the high-speed rain- 
erosion problem should be sought among the soft 
rubbery materials or among the hard rigid materials 
can be answered in the following way. The material 
may be hard and rigid, or it may be soft and rubbery 
and still have a high degree of rain-erosion resistance, 
If it is hard and rigid, it must have strength proper- 
ties that will enable it to withstand, without fracture 
and without plastic flow, the maximum unmitigated 
stresses imposed by collision with a waterdrop in 
the velocity range for which the material is being 
tested; these stresses increase in magnitude as the 
impingement velocity is increased. If it is soft and 
rubbery, it need only have strength properties suffi- 
cient to withstand the mitigated stresses that are 
imposed by collision with the waterdrop. The 
velocity ceiling for the rubbery materials is the point 
at which they are no longer able to withstand the 
mitigated stresses. 

These remarks apply to the ability of a material 
to withstand a single waterdrop blow. Under actual 
flight conditions successive blows on the same spot 
have a certain probability of being very closely 
spaced in time. Therefore, under flight conditions 
a rain-erosion resistant rubbery material must not 
only have strength properties adequate to withstand 
the mitigated but must also recover fast 
enough to be able to mifigate the stresses of an addi- 
tional blow to the extent that they do not exceed its 
strength properties. To be a practical rain-erosion 
resistant material the rubber must, furthermore, not 
lose its ability to mitigate the stresses through de- 
terioration of any kind during its expected service 


life. 


stresses 


WaAsHINGTON, April 9, 1957. 
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Emission Spectrum of Carbon Monoxide From 2.3 to 2.5 
Microns’ 
Earle K. Plyler, Harry C. Allen, Jr., and E. D. Tidwell 


The 2-0, 3-1, 4-2, 5-3, and 6 


in the region from 4000 to 4360 em~ with a grating spectrometer. 
second order resolved lines separated by 0.08 em. 
measurements were made on the 2-0 band and the molecular constants were calculated. 
0.000025 Do 


grating used in the 


constants in em are Bo=—1.922511 


tf emission bands of carbon monoxide have been measured 


A 10,000 line per inch 
Several sets of 
The 
Using these 


6.13+ 0.02 X 10-°. 


values and the By from 1-0 band in conjunction with the microwave constant a new deter- 


mination of the speed of light is obtained, ¢ 
1. Introduction 


The vibrational-rotational spectrum of carbon 
monoxide in emission was first measured in the infra- 
red region in 1951 |1].2. The observations were made 
on the harmonic band in the region from 4000 to 
4360 cm~'. Due to the high temperature of the 
flame the 2-0), 3-1, 4-2, and 5-3 transitions were ob- 
served. Recently Goldberg and Muiler |2] have 
measured the absorption spectrum of carbon monox- 
ide in the solar spectrum and have observed rota- 
tional lines with J values up to 75. 

Other bands of carbon monoxide have also been 
measured. Herzberg and Rao [3] have measured 
the 4—0 band of CO in absorption in the near infrared, 
and the fundamental band at 4.67 uw has been meas- 
ured by Plyler, Blaine, and Tidwell [4]. More re- 
cently 33 lines of the 2-0 band of CO have been 
measured by Rank and his co-workers |5)}. From 
each set of measurements the various workers have 
calculated a set of molecular constants for the CO 
molecule and the results of different authors [2, 4, 5] 
are in close agreement. 

Since the earlier measurements were made on the 
harmonic band of CO, improvements in detectors 
and infrared spectrometers permit more accurate 
measurements to be made. The bands in the 4000 
to 4360 cm™ region have been remeasured and this 
paper gives the results which have been obtained. 


2. Experimental Method and Results 
The emission from carbon monoxide in the 2.3-u 
region was produced in two ways. First, a tube con- 
taining CO was placed in the coil of a hf oscillator. 
(The use of a hf oscillator for the production of 
molecular spectra has been previously studied by 
Wilkinson, Ford, and Price [6].) The pressure of 
the CO was 6 em and a small flow of fresh gas was 
maintained. A spectrum of high intensity was ob- 
served showing transitions up to 11-9. The inten- 
sity of the rotational lines, however, varied with the 
pressure and the current through the oscillator coil 
and it was felt that more precise measurements of the 
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99,794 


3 kilometers per second. 


lines could be made of the spectrum from a flame. 

Second, the CO spectrum was observed in emission 
from an oxyacetylene flame. The gases were pre- 
mixed and were slightly rich in fuel. The burning 
took place in a welding torch and a section of the 
flame about 2 in. above the inner cone was focused 
on the front slit of the spectrometer. The radiation 
falling on the slit was increased by reflections of 
mirrors around the flame [7]. A grating spectrom- 
eter with a collimating mirror of one-meter focal 
length and a grating of 10,000 lines/in. was used for 
resolving the bands. The spectrum was observed 
in the first and second orders, the second order giving 
higher resolution. Emission lines separated by 0.08 


cm! were partially resolved. 

The observed spectrum (4000 to 4360 cm) 
includes the transitions 2-0, 3-1, 4-2, 5-3, 6-4, and 
7-5. Because of the ove rlapping of the lines from 


the different transitions, precise measurements were 
made on the bands arising from the three lowest 
transitions. In figure 1 are shown the band heads 
and some of the rotational lines for the 2-0, 3-1, 4 
5-3, and 6-4 transitions. The intensity of the 
rotational lines with the same J values of the 2—0 
and the 3-1 is about equal. In the other transitions 
the intensity is lower. In order to obtain sizeable 
deflections it was necessary to increase the amplifier 
gain for the other transitions; the noise level was the 
greatest for the 6—4 band. 

The small lines appearing between the strong lines 
from R30 to R45 in each band arise from the band 
folding back on itself; these lines are part of the R 
branch extending up to about J=68. The high J 
values are almost superimposed on the lower J’s 
in the 2—0 transition, but in the 3-1, 4-2, and 5-3 
transitions the high J lines are well separated from 
the stronger lines in the region R30 to R45. 

Some of the low-intensity lines observed in the 3-1 
region are the R3 to R6 lines of the 2-0 band. Their 
low inte nsity shows that very few of the molecules 


in the flame exist in the lower rotational states. 
On figure 1 are given the identifications of the 
different lines. Precise measurements were made 


on the 2-0, 3-1, and 4-2 bands. When the lines 
of the different series overlapped they were not 
measured. In order to have a more complete set 
of measurements of the 2—0 band the lines from R20 





to P20 were measured in absorption. The caleu- 
lated and observed wave numbers of the 2-0, 3-1, 
and 4—2 bands are given in table 1. 


3. Discussion of Results 


The molecular constants of CO have been caleu- 
lated from measurements of the 4—0. 2-0. and 1-0 
bands by several observers [1, 2, 3, +, 5). In most 
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cases the constants determined by 
Muller; Plyler, Blaine, and Tidwell: and 
agree within the probable error given by the 

The line positions measured in the 2 
been used to calculate Bo, Dy and 
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authors, 
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The constants were calculated both with and without The first value was obtained from the least Squares > = Ji 
the m‘ term. Only very small differences occurred treatment of 13 A,F’’ values from the measurements 
in %, Bo, and Dy when the fourth power term was’ and the second value was the result of using a ¢g). 


omitted. Since D,—D, is about 210~* em™', the culated value for the constant J) and including 
contribution of this term for J=20 is approximately only the 8 largest A,F’’ values in the least squares 


0.001 em~' which is about one-fourth the experi- treatment. Rank et al. [5] have calculated the 
mental error in the measurements of the lines. velocity of light using the second By given aboye 
Sixty-five lines of the 2-0 band were used in the and from By=57,635.965 +0.005 Me determined | 
calculation. They extended from R46 to P24 (see from the rotational spectrum. The result is ¢- 
table 1). The constants obtained, in cm™', were 299,793.7 +0.7 km/sec. On using the same value of 
By from microwave measurements and the average 
vy = 4260.069 +0.001 value from our infrared measurements, the speed of 
B.—1.922511 -+0.000028 light is determined as ¢ 299,794 13 km sec in 
excellent agreement with previous determinations 
D,=6.13 +0.022 107°. by this method. The central value is somewhat 
higher than the best value given by Bearden and 
The constants By and D, had been previously calcu- Thomsen [10], 299,792.8 and also that given by 
lated from measurements of the 1-0 band [8] and Cohen, DuMond et al. [11], 299,793.0; however 
found to be both of their values are within the probable error 
of this determination. 
Bo=1.922523 +0.000025 
D,=6.26 +0.17 107°. 


In conclusion, the authors thank Carroll Danne- 
miller for doing the least squares calculations on 


The standard deviations were first reported but the 
the Bureau’s electronic computer. 


uncertainty is now given as probable error. The 
values of B, from the two bands fall within the 
error of measurements. A better value of B, can be 4. References 
obtained by combining [8] the two 4,’s determined 


. ° ° > - » r : ene . “are - silwe an 
by analyses of the 1-0 and 2-0 bands, which gives — |"! Earle K. Plyler, W. S. Benedict, and 5. Silverman, J. j 


Chem. Phys. 20, 175 (1952). 


the result, [2] Leo Goldberg and Edith A. Muller, Astrophys. J. 118 
397 (1953). 
B, 1.922517 -O.000019, [3] G. Herzberg and Kk. \ Rao, J. Chem. Phys 17, 1099 
1949). 
the uncertainty again being the probable error. [4] Earle IK. Plyle r, Lamdin R. Blaine, and Eugene D 
nn ee > + : . ridwell, J. Research NBS 55, 183 (1955) RP2617 
The central value of Vy as determined from an [5] D. H. Rank, A. H. Guenther, G. D. Saksena, J. N 
analysis of the 2—0 band is in better agreement with Shearer, and T. A. Wiggins, J. Opt. Soc. Am. 47, 686 
the recent value of Do=6.13110~° obtained by . 1957). 
Gordy from the microwave spectrum [9], and the [6] Wilkinson, Ford, and Price, Molecular spectroscop; 
: ; : oa ae . (Institute of Petroleum, London, 1955 
value of Dyo=6.117X10 from a& measurement of [7] Earl K. Plyler and H. J. Kostkowski, J. Opt. Soc. Am 
the 2-0 by Rank et al. [5] than was the earlier de- 42, 360 (1952). 
termination from the ] () band. [8] Earle Kk Plyler, L. R Blaine, and W S. Connor, J, 
pe . ‘ . _ . Loe 4 ) 4; 
There are now two determinations for By which Opt. Soe. Am. 45, 102 (1955 
ee chs f 1.999517-+-0.000019. Thi [9] W. Gordy (private communication). 
eve an = rage value of 1.922. hi J WOLD. MS {10} J \. Bearden and J. 8S. Thomsen, A survey of atomi 
value of B, Is 1 excellent agreement with the \ ilues constants (Johns Hopkins University, 1956). 
obtained bv Rank et al. [5], ill] > R. Cohen, J. W.°M DuMond, et al., tev. Mod 


Phys. 27, 363 (1955). 
By=1.922519 + 0.000013 


and 
Bo= 1.922521 +0.0000035 WasHINGTON, February 5, 1958 
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Further Studies of the Influence of a Ridge on the 
Low-Frequency Ground Wave 


James R. Wait and Anabeth Murphy 


Computations are presented in graphical form for the perturbation of a plane wave by 
a semicylindrical boss on an otherwise flat ground plane of perfect conductivity. The 
height of the ridge is comparable to the wavelength. This is an extension of earlier work 
on the semielliptical boss 


A suggested model to study the effect of a nonsmooth ground for low frequencies is a 
semielliptical or semicylindrical boss on an otherwise flat ground plane of perfect conductivity. 
Using this model of a ridge, calculations were presented previously ' for a base width of two- 
thirds of a wavelength and various ellipticity ratios. Recently, the computations for the semi- 
evlindrical case have been extended for other diameters. These results are presented here. 

In the special case of the semicircular boss, the elliptic wave functions reduce to Bessel 
functions. The revelant solution is given in outline in the following. 











FiGuRE 1. Amplitude and phase of the field in the vicinity of a semicilindrical boss. 


Semicylindrical boss on flat ground plane and incident plane wave 


The perfectly conducting semicylindrical boss of radius @ is illustrated in figure 1, along 
with the cy lindrical coordinate system (p, 9, 2). The incident wave is incident from the right, 
and the magnetic vector has only a z component H/?* given by ? 


iT H,et*e 28 ¢, (1) 
where #7, is the amplitude of the magnetic field and k=2z/wavelength. The associated electric 


field is vertical and has an amplitude /. Using a well-known addition theorem in Bessel 
functions,’ the incident field can be rewritten in the form 


}7=-— FH, > > ce, "* *? J, (kp) cos md, (2) 
m v0 
where 
@=1, é&a=2(m=—1,2.3.. .). 


The secondary, or scattered, field #7. is now chosen to satisfy the boundary condition 
that the normal derivative of the total field is zero at p=a for r>@>0 and at @=0, 7 and 
p>a. Furthermore, the secondary field must contain the Hankel function H,? (kp) to comply 


9797 


tJ. R. Wait and A. Murphy, Influence of a ridge on the low-frequency ground wave, J. Research NBS 88, 1 (1957) RP2727 
> The time factor implied is exp(iwt 
G. N. Watson, Theory of Bessel functions, Cambridge University Press, 2d ed. (1944) 
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with the radiation condition at infinity. It then follows that 


, H ~y (kp) J,.(x) 
se _ — imr /2 m r 
HH: H, Do Ent 2") 


COS 171.0, (3) 


where 
r=ka=2ra/d, 


and where the prime indicates a derivative of the Bessel function with respect to its argument. 
The “‘relative amplitude”’ of the field is then defined by the quantity 
Hy +H? 
H, 
The “phase lag’ of the field is defined by 
] , Hy +H? 
cp cos @—arg . 
| H, 
In other words, the relative amplitude is the magnitude of the total magnetic field for an 
incident wave of unit amplitude. The phase lag, as defined above, is the argument of the 
incident field less the argument of the total field. 
| 


The “relative amplitude” and “phase lag’”’ for the magnetic field on the boss (p=a, 180 
@>0°) are plotted in figures 2,a and b, respectively, for various values of z(—ka). The 
abscissa in both cases is the angular positon ¢@ on the boss. It is immediately apparent that the 
amplitude of the field on the “illuminated” side (¢<.90°) of the boss is approaching 2.0 as the 
cylinder diameter is increasing. The amplitude of the field on the shadow side (od 9O~) of the 
evlinder is strongly reduced for the larger cylinders. It is also interesting to note that the 
“phase lag’ on the illuminated side is negative corresponding to a phase lead relative to the 
incident field while it tends to rise rather sharply on the shadow side. Geometrical optics 
would predict an amplitude of 2.0 and a zero phase for 0<¢<90° and a vanishing field for 
90°<@<180°. ‘This is the asymptotic behavior as x tends to infinity. 


Ficure 2 in plitude and phase of the field in the — 
vicinily of a semicylindrical boss. 


: Amplitude of field on bo I Phase lag of field on be 


The “relative amplitude” and ‘“‘phase”’ of the magnetic field on the ground plane in front 
of the boss (p>a,@=0°) are shown plotted in figures 3,a and b, respectively, for a range of z. 
The abscissa is the distance measured from the center of the boss and expressed in radians. 
Here it is apparent that the secondary, or scattered, field is interfering with the incident wave 
to form a standing wave pattern. As would be expected, the period of the oscillation is ap- 
proximately w radians or one-half wavelength. 
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FIGURE 3. Amplitude and phase of the field in the vicinity of a semicylindrical boss. 
I / J ' A ! 


Amplitude of field in front of boss (¢=0°). (b) Phase lag of field in front of boss (¢=0°). 
The same quantities for the ground plane to the rear of the boss (p>a,¢ 


180°) are shown 
in figures 4,a and b. 


It is seen that the relative amplitude of the field is approaching its 
unperturbed value of 1.0 as the observer recedes from the ridge and the phase lag is being 
restored to its unperturbed value of zero. 
more pronounced influence. 

The above model is admittedly highly idealized. The restriction, however, that the 
conductivity of the ridge is infinite can be easily removed in the formulation. 
condition on the evlinder is then to be replaced by 


As would be expected, the larger cylinders have a 


The boundary 


ER, —ZH., (4) 


where Z is the surface impedance.* 


‘M.A. Leontovich, Approximate boundary conditions for the electromagnetic field on the surface of a good conductor, Bul. Acad. Sci., U. 8. 
rie physique, in Russian, 9, 16, (1944). 


Z for a good conductor, expressed as a ratio of the 
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FIGURE 4 ly plitude and phase of the field in the 


mcinily of a semiculind ical boss 


intrinsic impedance of free space, is given by 
Z/1 202 ™ (14-1) (€w/2e)!. 5 


€) 1S the dielectric constant of free space, and o is the conductivity of the evlindrical boss or 


ridge. Using this boundary, condition, the secondary field becomes 


J, (2) +G,,(a 
ITs 7. > .40"" —, . ITS (kp) cos me 6 
> € He (x) GHe a) p) co 7 
where 
1Z 
1207 eo 
When |G< <1 corresponding to a highly conducting evlinder, the above formula for //= 


becomes indistinguishable from its perfectly conducting counterpart given by eq (3). 

The assumption of a perfectly conducting ground is not as simple to remove. However, 
if attention is confined to within several wavelengths of the ridge, the incident ground wave 
would behave locally as a plane wave at low-radio frequencies for a moderately or a well- 
conducting ground. 

The restriction to normal incidence of the plane wave on the perfectly conducting cylin- 
drical boss is easily removed. The only change necessary is to replace k by k sin 6 everywhere 
(@ is the angle subtended by the incident wave normal and the axis of the cylinder). The 
parameter z in the figures is then equal to ka sin @ and kp is replaced by kp sin @. 

Ground-wave propagation over ridges or obstacles of finite length could be treated in the 
same fashion, using semispherical or semispheroidal bosses, although the computations would 
become excessively tedious. The results presented here for the semicylindrical boss and the 
earlier ones for the semielliptical boss are believed to be satisfactory to \ ield an initial estimate 
of the magnitude of the effect. Subsequent computations for more elaborate models would be 
justified if careful experimental measurements could be carried out under controlled conditions. 


Bovu.per, Coxo., March 13, 1958. 
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Infrared Studies on Polymorphs of Silicon Dioxide and 
Germanium Dioxide 
Ellis R. Lippincott,’ Alvin Van Valkenburg, Charles E. Weir, and Elmer N. Bunting 


The infrared spectra of coesite, low-temperature tridymite, low-temperature cristobalite- 
low-temperature quartz, vitreous silica, hexagonal GeQOs, tetragonal GeO», and vitreous ger- 
Wherever possible an assignment of frequencies 


mania are reported from 4,000 to 300 em~'!, 


has been made on the basis of the selection rules for the crystal symmetry. 


Three character- 


istie group frequencies near 1,100, 800, and 480 cm~! are common to all the polymorphs of 


S10, 


These frequencies respectively correspond to a stretching mode involving displace- 


ments associated primarily with the oxygen atoms, a stretching mode involving displacements 


associated primarily with the silicon atoms, and a Si-O bending mode. 


The presence of 


these group frequencies in coesite indicates that the coordination of silicon in coesite is tetra- 
hedral and that its high density is associated with the packing of tetrahedral units at an 


angle approximating 120 degrees. 


The tetragonal and hexagonal GeO, polymorphs show a 


marked difference in spectra due in part to the change from sixfold to fourfold coordination. 
The assignment of observed frequencies in hexagonal GeO, is consistent with that made for 
low-temperature quartz if allowance is made for the beavier mass of the Ge atom. 


1. Introduction 


Several investigators have presented a detailed 
analysis of the vibrational spectrum of quartz [1]. 
Infrared spectra of cristobalite and fused silica have 


been reported but no specific assignments of their 


spectra appear to have been viven 11] Infrared 


spectra are presented of the SiO, polymorphs of 


coesite, cristobalite, tridymite, and vitreous silica; 
also the GeO, polymorphs of hexagonal GeQkg, 
tetragonal GeQO,, and vitreous GeQ,, in the region 
from 4,000 em™ to 300 em A partial interpreta- 
tion of the spectra in terms of respective structures 
and selection rules is given. Particular emphasis has 
been placed on locating characteristic frequencies and 
on classifving the observed frequencies in terms of 
types of vibrations. The assignment of frequencies 
is admittedly tentative in part because either the true 
erystal symmetry is not known, the spectral data are 
incomplete, or the spectra are too complex. Because 
coordination of silicon in coesite }2 has not been es- 
tablished, one of the objectives of this work was to 
furnish evidence for either the fourfold coordination 
of the other polymorphs of SiO, or for a different co- 
ordination. In the tetragonal form of GeQ, the 
germanium atom has a sixfold coordination with oxy- 
gen, Whereas in the quartz-type hexagonal GeQ, the 
germanium atom has the fourfold coordination of 
silicon in quartz [3, 4 A similar tetragonal form of 
silica is not known but has been predicted by analogy 
with germania. 


2. Experimental Procedure 


2.1. Samples 


The cristobalite was prepared from ethyl ortho- 
silicate by the usual hydrolysis technique. The SiO, 
thus obtained was heated to 1.600°C for 5 hr. Its 


Consultant to the Miner I jucts D on, National Bureau of Standards 
ve Department Chemistry, University of Maryland, College 


Figures in bracket licate the literature references at the end of this paper 
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X-ray powder pattern was identical with that found 
for natural cristobalite. A quartz sample of high 
purity was obtained from the Geophysical Labora- 
tory (standard sample). Crystals of coesite were 
prepared from precipitated silicic acid which was sub- 
jected to 30,000 atm at 600°C for one-half hour. 
The X-ray powder pattern and indices of refraction 
of this material agree with those reported by Coes 
(2|. Tridymite was prepared by C. N. Fenner of the 
Geophysical Laboratory. The preparation was ac- 
complished by heating pure quartz and Na,WQ, at 
1,000° C for 71 hr, followed by careful purification of 
the tridymite. The fused silica was obtained by 
fusing a sample of pure quartz in the usual manner. 

The germanium oxide was obtained commercially. 
Microscopic examination showed extremely fine par- 
ticles less than 1 uw in size. X-ray powder patterns 
of the material coincided with those reported for the 
hexagonal form of GeO, [4]. The tetragonal form 
of GeO, was prepared hydrothermally. The hex- 
agonal GeO, was sealed with water in a platinum 
capsule and heated at 800° C at a pressure of 20,000 
psi for 3 days. Microscopic examination showed 
some crystal growth of tetragonal GeO,, but crystals 
were still too small (cirea 1 ») to permit more than an 


estimation of the index, approximately 1.99. X-ray 
patterns showed tetragonal GeO, with a small 
amount—estimated at less than 5 percent—of un- 


converted hexagonal material. Spectroscopic anal- 
vsis of the tetragonal material showed principal im- 
purities to be Pt 0.1 to 1.0 percent; Al, Ca, Fe, Pd, 
and Si 0.01 to 0.1 percent; and smaller quantities of 
Ag, Ba, Cu, Mg, Mn, and Pb, which were found as 
traces. Vitreous GeO, was formed by quenching 
fused hexagonal GeQ, in a platinum container. 


2.2. Infrared Studies 


Infrared spectra of the samples were obtained on 
a double-beam spectrophotometer with NaCl prism 
in the region 4,000 to 650 em™ and on a single-beam 
double-pass spectrophotometer with CsBr prism in 


the region 700 to 300 cm. The samples were 








ground in a boron carbide mortar, mixed with dry 
ground KBr and pressed into KBr pellets at 100,000 
psi by standard techniques [5]. The concentrations 
range from 2 to 3 mg of sample per gram of KBr. 
The spectra from 4,000 em to 300 cm™! are given 
in figures 1 to 15. The region above 1,500 em™! 
contained no fundamental frequencies for any of 
the polymorphs of SiO, or GeO,. Absorption peaks 
near 3,500 cm~! were due to water present in the 
KBr. The observed frequencies given in figures 
1 to 15 are tabulated in tables 2, 4, 5, 7, 8, and 9. 


3. Interpretations 


A SiO, group in a rigid framework has 3N-3 or 
six vibrational degrees of freedom. There is some 
ambiguity as to how the corresponding modes of 
vibration should be classified in terms of bond 
stretching, bond bending, and bond-distortion types. 
Somewhat arbitrarily a bond-stretching and bond- 
bending mode will be associated with each O atom. 
These O-stretching modes should correspond to the 
highest observed frequencies in the spectrum. Two 
remaining degrees of freedom need to be assigned 
to the Si atom. One of these must be a stretching 
mode corresponding appreximately to motions of 
the Si atom between the two O atoms. The one 
remaining mode must correspond to a low-frequency 
distortion or bending motion of the Si atom. 

For each SiO, group which is added to the first to 
make up the unit cell, one must add nine vibrational 
degrees of freedom. These nine modes may _ be 
classified approximately into the following types: 
Because each added SiQ, group corresponds to the 
addition of four bonds, four of the nine modes 
correspond to bond-stretching frequencies, two in- 
volving displacements associated with the O atoms 
and two for the Si atom. (Using another descrip- 
tion these would correspond to two antisymmetric 
stretching modes of the type <-Si O Si, and two 
symmetric stretching modes of the type —SiO Si 
Two more of these nine modes may 
as bending motions associated with the bending of 
Si-O-Si angle. The remaining three modes are 
associated with low-frequency bending or distortion 
frequencies. Table 1 summarizes this classification 
of the modes of vibration for various numbers of 
SiO, groups per unit cell. The case for N equal to 


be classified 


mode of vibration 


ini cell 


TABLE 1. 
for diffe rent numbers of SiO» q Oups per 


Descriptive classification of the 


510: groups per unit ce Sper t nh as- 
gned poly- 
mor} SiO» 

l 2 $ 1 t s 

Number of Si-O stretch- 2 4 6 8 12) le 1 20K 1 OO 
ing modes involving 
motion primarily 
sociated with the 
oxygen atoms. 
Number of stretching l 3 7 | l 82 Ht 
modes involving 
motion primarily 
associated with the 
Silicon atoms 
Si-O bending modes 2 ' 6 8 | 12 600 m 
Distortion mode l { 7\/10> 1 22 380 to 100 cn 
Total mode 6/15 24 33 51 | 68 
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3 corresponds to quartz and that for N equal to 8 to 
cristobalite. It should be emphasized that. this 
classification into types is only approximate and 
that in general a given mode will be a mixture of 
stretching, bending, and distortions types. The 
proper method of classifying these modes would 
be in terms of the symmetry elements of the crystal. 
For quartz this has been impressively illustrated 
by Saksena [1]. Because there is some doubt ag 
to the true crystal symmetry of cristobalite, and 
that of coesite [2] has not been described in detail, 
this method of classifying normal modes on the 
basis of symmetry is not possible at the present 
time for these latter two forms of silica. 

The feasibility of the method outlined above for 
the classifications of the normal modes of various 
forms of silica may be tested on the assignment given 
for quartz by Saksena [1]. From table 1 one would 
expect, in order of decreasing frequency, four dis. 
tinct spectral groupings consisting of six Si-O streteh- 
ing modes involving displacements associated pri- 
marily with the O atoms, five stretching modes in- 
volving displacements associated primarily with the 
silicon atom, six Si-O-Si bending modes and seven 
low-frequency distortion modes. These are readily 
associated with the frequencies observed in the fol- 
lowing regions of the spectrum: 1,200 to 1,000 em=. 
825 to 600 ecm™', 600 to 390 ecm~', and 380 to 106 
em~', respectively. The distinction between the lat- 
ter two groupings is somewhat arbitrary. If one 
allows for the fact that frequencies assigned to species 
E in table 2 must be counted twice because they are 
doubly degenerate, it is found that these groupings 
contain 6, 5, 6, and 7 degrees of freedom, respectively, 
in agreement with the numbers expected from a unit 
cell containing three SiO, groups. 


3.1. Cristobalite 


The observed spectra were taken on low-temper- 
ature a-cristobalite. Since little information is avail- 
able about its structure, the spectra are interpreted 
on the basis of that for high cristobalite for which 
more structural information is available [6, 7}. High 
cristobalite was studied by Wyckoff for which he 
proposed a cubic structure of space-group svmmetry 
Of=Fd3m with eight SiO, groups per unit cell {6}. 
Application of O, selection rules shows that only 5 
triply degenerate fundamental frequencies of species 
F. corresponding to 15 of the 69 modes of vibration 
for the unit cell should appear in the infrared spee- 
trum. Because the observed spectrum of low-tem- 
perature cristobalite contains at least seven frequen- 
cles plus additional low -frequency peaks below 300 
em~', it is not feasible to interpret the infrared spee- 
trum on the basis of Of—=Fd3m _ space-group sym- 
metry. Barth [7] reexamined the structure of high 
temperature 8-cristobalite and concluded that its 
structure corresponded to a lower symmetry and 
proposed a structure with space-group symmetry 
‘ P2,3. Application of T-selection rules to a unit 
cell with eight SiO, units predicts that seventeen 
triply degenerate fundamental frequencies should be 
observed in the infrared spectrum. At least twelve 


| 





of t 
cm 
inte 
basi 

T 
assis 
may 
whit 
as @ 
stru 
rule 
The 
met 
spec 
men 
prec 
abo 
clos 
In t 
stre! 
freq 
spec 
Si-O 
obse 
tivel 
disp 
the | 
at | 
benc 
obse 
acti 
200 
class 
T, : 
the 
that 
assig 

A 
crist 
spec 
Wav 
anal 

A 
crist 
of t 


to 8 to 
it this 
te and 
ture of 

The 
would 
‘Tystal, 
strated 
ubt as 
e, and 
detail, 
on the 
present 


ve for 
Various 
t given 
would 
ur dis- 
tretch- 
“d pri- 
des in- 
ith the 
seven 
readily 
he fol. 
) em 
to 100 
he lat- 
If one 
species 
ley are 
upings 
tively, 
a unit 


pm per- 
; avail- 
preted 
which 

High 
ich he 
metry 
ell [6). 
only 5 
species 
ration 
| spec- 
v-tem- 
pquen- 
w 300 
| spec- 
) svym- 
f high 
iat its 
v and 
metry 
a unit 
pnteen 
uld be 


twelve 


TABLE 2. 


Spectral activity Saksena * 


Kirshman * 


Spec 1eS 
207 207 
A Raman 356 258 
| 466 467 
1, O82 1, O82 
| 364 378 
with no « "I 00 480 
B Infrared with no dipole ’ Z ! oo 480 
1, 149 
f 128 127 
392 to 403 104 
695 696 
1, 159 1. 160 
! Ed Infrared with no dipole ” Z ¢ 265 267 
479 
795 to 807 704 
1, 063 1, 065 
1, 228 


« See reference in text. Figures in parentheses indicate the 
doubly degenerate and the observed frequencies correspond to two moc 


| of these should be observed in the region above 300 
1 and it is considered unlikely that a satisfactory 


cm- 
interpretation of the spectrum can be made on the 
basis of space-group symmetry T*=P2,3. 
ese considerations sugges iat a tentative 
T lerat uggest that tentat 


assignment and classification ‘of observed frequencies 
may be made on the basis of a different symmetry 
which has vibrational selection rules not as stringent 
as an QO, structure but more restrictive than a T 
structure. To accomplish this, the Tg selection 
rules are assumed to be of approximate value [8]. 
The vibrational species for a structure of Tg sym- 
metry are A,, A,, E, F; and F,. Of these the only 
species which is infrared active is F,. Ten funda- 
mental triply degenerate infrared frequenci les are 
predicted, of which eight should appear in the region 
above 300 cm This latter number corresponds 
closely with the number of observed frequencies. 
In table 3 are summarized the distribution of Si-O 
stretching, Si stretching, Si-O bending, and low- 
frequency distortion modes among the vibrational 
species for a structure of Tg symmetry. The three 
Si-O stretching modes may be assigned to the peaks 


- 


observed at 1,204, 1,160 and 1,104 em™', respec- 
tively. The two Si stretching modes involving 


displacements associated primarily with motion of 
the Si atoms may be assigned to the peaks observed 
at 798 and 620 em", respectively. The three 
bending modes have been assigned to the bands 
observed near 515 and 485 em™!. The two infrared 
active distortion modes are expected in the 150 to 


200 em=' region and would not be observed. The 
classificaton given in table 4, based on assumed 
I; selection rules and group frequencies, explains 


the fact that the observed spectrum is simpler than 
that of quartz and enables one to obtain a tentative 
assignment of all the observed frequencies. 


A more detailed and complete assignment for 
cristobalite must await the obtaining of a Raman 


spectrum and additional infrared studies at longer 
wavelengths on single crystals, followed by an 
analysis based on the true crystal symmetry. 

A comparison of the spectra of quartz and of 
cristobalite shows that the most intense component 
of the Si-O antisymmetric stretching frequencies 


use of Saksena’s assignment. 


les. 
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Classification of fundamental frequencies (em) in a-Quartz 


Barriol * Simons and | This work ° | Descriptive 


Me Mahon assignment 
207 (207) Distortion. 
356 (356) Distortion. 
466 (466) Si-O bending. 

1, O82 (1, 082) Si-O stretching. 
385 374 Distortion. | 
4188 513 Si-O bending. 

800 790 780 Si stretching. 

1, 190, 1, 227 1, 055 1, 150 Si-O stretching. 
127 (128) Distortion, 
304 397 Si-O bending. 
696 695 696 Si stretching. 

1, 163 1, 162 1, 176 Si-O stretching. 
268 (265) Distortion. 

462 to 475 Si-O bending, 
798 795 to 802 801 | Si stretching. 
1, 065 1, 097 Si-O stretching. 

1, 228 


© Direction of optic axis. 4 Species E is 


TABLE 3. Section rules and group frequency classification for 


a Tq structure with 8 SiO, groups per unit cell 


Total Si-O Si Si-O Low fre- 
Species * Spectral fre- stretch- | stretch- | bendings| quency 
activity quencies ings ings distor- 
tions 
Ay Raman... 5 2 1 2 dotosinindliti 
A» Inactive. 1 . ‘ 1 
E Raman. 6 1 1 1 3 
Fi Inactive 7 1 2 1 3 
F2. Infrared 10 3 2 3 2 
and 
Raman. 
“The E and F species are doubly and triply degenerate and the corresponding 


frequencies must be counted two and three times, respectively. 


Classification of observed infrared frequencies for 
cristobalite * 


TABLE 4. 


em~! Description 
485 Si-O bending. 
515 Si-O bending. 
620 Si stretching. 
798 Si tretching. 

1, 104 Si-O tretching. 

1, 160 Si-O stretching. 

1, 204 Si-O stretching. 


® F: is the only active infrared species 
for Ta symmetry. The Raman active 
species are Aj, Ey, and F9. 


differs by only 7 em™', with the peak for cristobalite 
being at the larger cm™! value. This result is some- 
what different from that found by Simmons and 
McMahon [1] with infrared reflection studies, where 
the Si-O stretching frequency was found at different 
cm”! values for quartz, cristobalite, and fused silica. 
Likewise the Si stretching frequencies near 800 cm™! 
do not differ greatly between the two forms. The 
only other frequencies that appear to be somewhat 
the same in both forms are those that correspond to 
the strong peaks occurring in the region 460 to 515 
cm”'. These frequencies have been assigned to a 
mixed Si-O-Si and O-Si-O bending mode. However, 
corresponding low-frequency peaks in quartz are 
definitely shifted to lower frequencies from their 
counterparts in cristobalite. This shift will be 








tentatively associated with the change in angle that 
the tetrahedral units make with each other in the 
respective structures. The fact that the Si and 
Si-O stretching frequencies occur at nearly the same 
positions, respectively, in both compounds indicates 
that the difference in structure between the two is 
associated with the orientation of the tetrahedral 
unit rather than any great difference between Si-O 
bond distances or Si-O bond-stretching force con- 
stants, resulting from a change of coordination. 


3.2. Coesite 


The infrared spectrum of coesite given in figure 3 
and tabulated in table 5 is richer and more complex 
than that for either quartz or cristobalite. This is 
expected since the cry stal svmmetry of coesite [2] is 
lower or C2/c) and 16 to 18 
SiO, groups per unit cell. Two factors make it 
difficult to obtain an assignment of frequencies for 
coesite. First, the exact crystal svmmetry is not 
known, nor can a pseudo structure of high symmetry 
be assumed to simplify the anaylsis. Second,the 
exact number of SiO, groups per unit cell is uncertain. 
However, the infrared spectrum of coesite shows the 
presence of group frequencies similar to those found 
The Si-O stretching fre- 


(space group Ce has 


in quartz and cristobalite 
quencies appear near 1,100 cm the Si 
stretching frequencies appear near 800 em This 
indicates that the tetrahedral unit in coesite Is still 
intact and that there has been no change 
nation of the Si atom or marked change in Si-O bond 
change in Si-O bond stretching- 
It is interesting that this is not the 


W hereas 


In) coordi- 


length or marked 
force constant. 
case for the hexagonal and tetragonal forms of GeO 
The difference in structure appears to be due to the 
packing of the tetrahedral units. Evidence for this 
is that the mixed Si-O-Si and O-Si-O bending fre 
quencies, which appear at 515 and 485 em™ in ertsto- 
balite and at 475 and 462 em“ in quartz, shift to 442 
and 430 in We interpret this shift to be 
associated with the change in angle between the tetra- 
hedral units of the three structures. Empirically one 
may estimate this angle in coesite by correlating this 
shift with the known angles in cristobalite and quartz, 


coesite. 


which are near 180° and 150 respectively This 
indicates that this angle in coesite is approximately 
1?0) 10 

TARLI ‘ ( f ( 


64 





3.3. Low-Temperature Tridymite 


Low-temperature tridymite forms orthorhombie 
crystals with 64 SiO, groups per unit cell [9]. The 
true crystal symmetry is not known, but it is believed 
to be a distorted form of the simpler high-tempera- 
ture 8-tridymite, which has a simple hexagonal 
structure of Ds Cé62ce or Ds C6/mme Space. 
groupsymmeiry. In order to obtain an approximate 
interpretation of the infrared spectrum of the low. 
temperature tridvymite it is assumed that the selee. 
tion rules for the high-temperature tridymite of 
D5, symmetry are approximately valid despite the 
fact that the true symmeiry is a distortion of the 
D. svimmetry. The infrared active species are A’ 
and E’ corresponding to nine fundamental fre. 
quencies (4.45 and 5K;). The D.», selection rules 
[8S] and a classification of Si-O stretching, Si stretch. 
ing, Si-O bending, and low-frequency distortioy 
modes are given in table 6. The two Si-O stretching 
modes in species A? have been assigned to the peaks 
observed at 1,109 and 1,175 em. The AS Si-O) 
bending mode has been assigned to the 478 em” 
peaks The distortion mode in AS would be ex. 
pected near 200 em and has not been observed 
The five frequencies of species K’ have been classi- 
fied as corresponding to an Si-O stretching, an § 
Si-O bendings, and a clistortion 


stretching. two 


mode. The Si-O stretching has been assigned at 
1.109 em the Si stretching at 792 em! and the 
two Si-O bendings to the 478 em peak and its 
nearby component. The low-frequency distortio 
would be expected in the 150 to 200 em recions 
and would be missed by the mstruments used in this 
work. 
| ABI a) S/ hior ‘ pid / neation o mode oO Drie 
na 81 group D4 —=C62e 
s I 

\ | 

\ 

‘ ' 

\ I 


This classification of frequencies based on selec- 
Cé2¢ space group enables one 
tentative assignment of all observed 
frequencies table 7 and Lo understand why the 
observed spectrum for tridymite is simpler than that 
found for both Quartz and cristobalite In particu- 
lar it explains the observation that only one fre- 
is observed in the 850 to 625 em=! region of 


iion rules for the DA 


LO obtain il 


quency 
the spectrum corresponding to a stretching mode 
associated with motion of the Si atoms 

Although the actual symmetry is lower this i 
dicates that the deviation from Do», selection rules 
is apparently not enough to cause a significant change 
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in the appearance of the infrared-spectrum., A more 
detailed analysis of the spectrum could be made if 
a Raman spectrum were available. _ Tt is interesting 
that the three observed frequencies in tridymite 
corresponding to an antisymmetric Si-O stretching, 
a Si stretching, and to a characteristic mixed O-Si-O 
and Si-O-Si bending, respectively, appear near the 
came em~! values in quariz, cristobalite, and coesite 
indicating that they may be classified as characier- 
istic group frequencies. The mixed-bending fre 
quency observed at 478 em7! corresponds to a tre- 
quency close to that observed for cristobalite. lhe 
slightly lower values observed in quartz and coesite 
are interpreted as being due to the change of angle 
of the tetrahedral units from 180 


TABLE 7. Classification of observed infrared frequen ves of 

tridymite on the basis of Ds, selection rules 
, Dype | 
vy 178 Si-O bending. | 
ng Si-O bending 

y2 Si stretching | 

n 1 1090 Si-O stretching 

| 7 | 
br 1 | Si-O stretching 

\ 

e E frequet ire doubly degenerate and each I 

to two modes of vibration 


3.4. Vitreous Silica 


The infrared spectrum of vitreous silica is given in 


figure 5. The spectrum is characterized by the 
presence of strong bands at 1,108 em 805 em, 
and 468 em! and may be assigned as the Si-O 


stretching, Si stretching, and Si-O-Si bending char- 
acteristic group frequencies, respectively. These 
frequencies are at essentially the same positions as 
in the other SiO, modifications, but they are some- 
what broader. The broad peak near 805 cm™! has 
its peak absorption at slightly higher em™ values 
than the other polymorphs of SiQ,. This would 
indicate that there again is no change in the coordina- 
tion of the tetrahedral units. The broadening of 
the bonds may be due to the random packing of the 
tetrahedral units. The infrared spectrum of vitreous 
silica resembles tridymite in appearance with the 
same number of observed frequencies. 


3.5. GeO, 


The infrared spectra of three different modifica- 
tions of GeQ, are given in figures 6, 7, 8, 14, and 15, 


corresponding io the tetragonal, hexagonal, and 
vitreous forms. The tetragonal form has a rutile 
structure with two GeQ, groups per unit cell, 


corresponding to sixfold octahedral coordination. 
The hexagonal form has a quartz-type structure with 
three GeO, groups per unit cell, corresponding to 
& tetrahedral coordination. This difference in 
coordination results in a large difference in frequency 
of the Ge-O-Ge antisymmetric stretching mode of 
vibration observed in both spectra. In tetragonal 
GeO, this frequency at 720 cm whereas in 
hexagonal quartz-like GeO, it is found at 885 em™, 


IS 
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a difference of 165 em™'. This difference is due to 
the markedly different Ge-O bond length and bond 
energy for the two forms. For tetragonal GeO,, 
the Ge-O bond length [3] is 1.86 A. The bond 
length for the hexagonal GeO, has not been given, 
but it is likely that it is shorter than 1.86 A, cor- 
responding to a stronger bond and larger bond- 
stretching force constant and vibrational frequency. 
The large difference in the Ge-O-Ge antisymmetric 
stretching frequency for those two forms of GeO, 
is in marked contrast to the positions of the Si-O-Si 
antisymmetric stretching frequency for the various 
forms of silica. There is only a small shift in the 
position of this frequency for the various forms of 
SiO, indicating that there is no change of coordina- 
tion and no, or only minor, changes in Si-O bond 
lengths, bond energies, and bond stretching force 
constants. The difference in structure of the 
various forms of silica must be associated with the 
packing of the tetrahedral units. 

The selection rules [8] for a tetragonal GeO, 
structure with Dj{=P4/mnm space-group symmetry 
and two GeQ, groups per unit cell predict the 
presence of six infrared active fundamental fre- 
quencies (3 A,, and 3 E,). The observed frequencies 
and assignment are given in table 8. The three 
frequencies of species A,, have been classified as 
antisymmetric Ge-O stretching, a Ge-O bending, 
and a low-frequency distortion mode. The first 
two modes have been assigned to the peaks observed 
at 720 and 513 em™'. The third mode of species 
A. is expected in the regions 125 to 200 cm™ and 
would not be observed with the instruments used 
in this research. The three frequencies of species 
EK, have been classified as a Ge-O stretching, a Ge 
stretching, and a Ge-O bending, and are assigned 
at 720, 525, and 425 cm™', respectively. No attempt 
will be made to explain the shoulders appearing on 
the 425 em™ band, 

TABLE 8. Classification of the observed infrared frequencies of 
tetragonal GeO» 


Species ® em! Type 

As Distortion. } 

Eu 425 | Ge-O bending. 

Agu... 513 Ge-O bending. 

Ee 515 Ge stretching. | 

Aou, Ew | 720 Ge-O stretching. | 
| | 


* Each E frequency is doubly degenerate and corre- 
sponds to two modes of vibration. 

A detailed interpretation of hexagonal GeO, 
would require the obtaining of an infrared spectrum 
run at lower frequencies as well as a Raman spectrum. 
However, the observed spectrum is consistent with 
that found for low-temperature quartz. The strong 
peak at 885 cm™ is clearly the analogue of the 
1,100 em Si-O-Si antisymmetric stretching fre- 
quencies of quartz in species B and E and must 
consist of two fundamental frequencies. The 801 
em! and 782 cm™! peaks of quartz correspond to the 
Si stretching modes and would be expected to shift 
markedly in hexagonal GeO, because of the heavier 
mass of the Ge atom. Tentatively these two peaks 








Classification of the observed infrared frequencies for 
hexagonal GeO, (cm) 


TABLE 9. 


| Species GeO, ype Quartz 
I yl 
| issignment 
' 
| A (Raman | 
| active | 
| Distortion 374 
| 345 Ge-O bending 513 
| 515 Ge stretching 780 
B (Infrared 885 Ge-O stretching 1, 150 
ictive 128 
| 307 
| 85 Ge stretching 696 
963 Ge-O stretching 1,176 
. 265 
F Infrared 332 Ge-O bending 462 to 475 
| active 540 Ge-stretching 801 
\ RSS Ge-O stretching 1, 097 


of quartz will be respectively associated with the peaks 
of hexagonal GeO, at 540 em™ and 515 cm. The 
first of these has been assigned to species E and the 
latter to species B of the quartz structure. The 696 
em~' band of quartz must then be shifted to 585 em™! 
in GeO, and assigned to species E. The two remain- 
ing frequencies observed in hexagonal GeO, at 345 
em! and 332 cm™ must be considered as the analogs 
of the mixed Si-O-Si, O-Si-O bending frequencies 
observed in quartz at 513 em™ and 462 cm, respec- 
tively. The first (845 cm™') must be assigned to 
species B, and the second (332 cm™') to species E of 





the quartz structure. This assignment is given jp 
table 9. 

The Ge-O-Ge antisymmetric stretching of vitreoys 
GeO, is found at 895 cm™ as a broad band and indi- 
cates its similarity with hexagonal GeO, in that both 
consist of tetrahedrally coordinated GeQ, units. No 
additional low-frequency bands were observed jp 
fused GeQy. 
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Infrared transmission spectrogram of quartz from Lisbon, Md. 
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